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ABSTRACT 


—  Motivated  by  a  lack  of  explosive  test  data  on  nuclear 
submarines,  the  Navy  has  sought  other  means  to  qualify  in¬ 
stalled  equipment  in  submarine  shock  environments.  The 
currently  used  method  for  non-shock  testable  items  is  the 
Dynamic  Design  Analysis  Method  (DDAM)  developed  by  the 
Naval  Research  Laboratory  in  the  early  1960's.  With  the 
advent  of  large-scale  computing  power,  newer  numerical 

methods  have  become  available  to  predict  equipment  re- 

. 

sponses.  This  i  nv^t-rg*44-on-  is  a  comparative  study  of 
DDAM  and  ELSHOl^;  a  new  generation  numerical  shock  response 
code.  The  limitations  and  strong  points  of  both  methods 
are  examined  using  illustrative  examples.  4 
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I .  INTRODUCTION 


A .  BACKGROUND 

Modern  combat  vessels  must  be  designed  with  the 
capability  to  survive  moderately  severe  shock  loadings 
induced  by  the  underwater  explosion  of  conventional  or 
nuclear  weapons.  In  the  scope  of  this  investigation, 
survival  includes  the  mission  survival  o-f  the  platform.  It 
must  maintain  its  ability  to  utilize  all  machinery  and 
weapons  systems  to  carry  out  its  primary  mission  following 
an  underwater  explosive  attack.  Although  the  vessel  may 
withstand  considerable  hull  damage  and  maintain  its 
structural  integrity,  failure  of  internal  equipment  can 
render  it  useless  as  a  weapon  and  thereby  eliminate  its 
value  in  time  of  crisis.  This  investigation  examines  the 
shock  loading  of  internal  equipment  in  submarines  and  some 
of  the  me thods. ava i 1 abl e  to  analyse  shock  response. 

In  view  of  the  fact  that  a  submarine  can  not  be 
expected  to  survive  a  direct  hit  by  any  modern, 
h i gh- i n tens i ty  weapon,  many  attempts  have  been  made  to 
specify  a  level  of  shock  loading  for  the  purpose  of  design 
evaluation.  The  current  specifications  for  building 
submarines  contain  the  shock  requirements  to  be  met  by  the 


builders  and  vendors  of  installed  equipment.  Generally, 
all  equipment  is  to  pass  a  series  of  shock  tests,  in  the 
installed  configuration  as  outlined  by  MI L-S-90 1 C  CRef.13, 
where  shock  testing  is  practical.  This  document  specifies 
standard  explosive  and  mechanical  test  fixtures  and 
procedures  for  conducting  the  tests.  The  suitability  of  an 
equipment  design  or  installation  is  evaluated  according  to 
its  ability  to  function  as  intended  during  and  after  each 
shock  impulse.  Equipment  testing  is  very  expensive  and 
requires  a  great  deal  of  preparation.  In  the  case  of 
equipment  design,  it  is  often  not  practical  to  evaluate  an 
installation  prior  to  the  final  design.  In  some  instances, 
the  direct  testing  of  installed  components  may  not  be 
possible  from  the  aspect  of  prohibitive  size  or  weight.  In 
these  cases,  a  design  analysis  is  specified. 

The.  design  analysis  for  shock  response  was  first 
proposed  in  the  United  States  during  World  War  II.  Hull 
damage  reports  indicated  the  degree  of  suceptibi! i ty  of 
installed  equipment  to  shock  loadings  from  bombs  and  depth 
charges.  Based  on  limited  testing  and  observations,  a 
shock  design  factor  was  specified.  A  series  of  curves  were 
used  to  depict  the  variation  of  shock  design  factor  with 
equipment  weight  to  yield  an  equivalent  static  design 
force.  This  force  was  used  to  specify  mounting  hardware 
and  main  structural  members  for  the  equipment  installation 


and  varied  with  motion  input  direction.  The  method  proved 
to  be  simple  in  application  but  reflected  a  lack  of  realism 
with  experimental  tests. 

In  the  late  1 950's,  nuclear  power  plants  were 
incorporated  into  submarine  designs  to  afford  them  greater 
independence  from  the  frequent  surfacing  requirements 
inherent  in  conventional  submarine  operations.  Nuclear- 
power  plant  technology  was  new  and  thus  had  no  previous 
history  of  shock  resistance.  The  acquisition  cost  of 
nuclear  components  is  high  and  the  new  dangers  associated 
with  their  damage  from  underwater  explosion  sparked  a 
renewed  interest  in  underwater  explosion  testing  arid  design 
analysis.  From  the  late  1950's  through  the  middle  1960's, 
most  of  the  currently  used  standards  for  submar i re¬ 
installed  equipment  shock  design  were  adopted. 

The  Dynamic  Design  Analysis  Method  (DDAM)  was  proposed 
by  the  Naval  Research  Laboratory  and  accepted  by  the  Navy 
as  a  design  evaluation  requirement  in  1963  [Ref  23.  The 
details  of  this  specification  will  be  discussed  more  fully 
later;  however,  the  basic  principle  is  stated  here.  DDAM  is 
a  simplified  modal  analysis  method  which  utilizes  shock 
inputs  which  were  empirically  derived  from  underwater- 
explosion  tests  of  realistic  ship  and  submarine  instal¬ 
lations.  It  is  assumed  that  the  equipment  and  its  foun 
dation  together  make  up  a  system  which  responds  as  a  linear 
elastic  structure  to  the  input  which  is  described  by  a 


design  shock  spectrum.  The  successful  uti i ization  of  DDAM 
to  evaluate  equipment  response  to  a  design  shock  input  is 
dependent  on  the  abil it/  of  the  design  shock  spectrum  to 
reflect  accurately  the  structural  environment  existing  in 
given  installation.  This  will  be  influenced  by  the  size 
and  weight  of  the  submarine  test  subject,  'the  structural 
system  employed  in  the  construction  of  the  submarine,  and 
the  interaction  between  the  equipment,  foundation  and  hull 
structure  used  to  derive  the  design  shock  spectrum. 

B.  THE  STATE  OF  CURRENT  SHOCK  DESIGN  ANALYSIS 

Since  its  adoption,  DDAM  has  been  used  without 
modification  to  qualify  ship  and  submarine  equipment 
installation  designs.  The  original  design  shock  spectrum 
[Ref.  23  remains  intact  without  revision.  Shortly  after 
the  adoption  of  DDAM  as  a  design  requirement,  underwater- 
explosive  testing  of  nuclear  submarines  was  banned  due  to 
the  inherent  risks  involved.  This  ban  was  relaxed  in  1983 
for  the  low  level  explosive  test  of  a  SSN  688  class 
su bmar i n e .  Me c h an i c a  1  and  explosive  tests  of  e q u i pme n t 
installations  on  shock  simulation  platforms  are  the  only 
regular  check  of  designs  meeting  the  standards  imposed  by 
DDAM.  A  danger  exists  with  this  kind  of  verification  in 
that  the  tests  may  be  designed  to  accurately  reproduce  the 
input  motions  provided  for  in  the  design  shock  spectrum 
utilized  with  DDAM . 


In  recent  years  the  increasing  availability  of  large 
computers  and  the  rapid  development  of  numerical  methods 
have  provided  the  engineer  with  new  tools  to  use  in  the 
analysis  of  submarine  shock  response  problems. 
Finite-Element/Finite-Difference  methods  allow  for  the 
analysis  of  structure  and  fluid  responses  and  the  work  of 
Geers  [Ref.  3]  provides  a  means  to  incorporate  fluid- 
structure  interaction  effects.  Several  computer  codes  have 
been  developed  to  analyze  installed  equipment  response  to 
shock  waves  based  on  these  formulations.  One  of  them  is  the 
Elastic  Shock  <ELSHGK>  code  developed  by  Weidlinger 
Associates  under  Defense  Nuclear  Agency  and  Office  of  Naval 
Research  funding  [Ref. 4].'  The  motivation  for  this  work  has 
been  to  investigate  modern  submarine  underwater  explosive 
shock  response  in  conjunction  with  a  testing  program  using 
small  to  large  scale  models  and  shaped,  explosive,  tapered 
charges.  Throughout  its  development,  ELSHOK  has  been 
validated  using  these  highly  controlled  tests.  In  the  low 
level  explosive  test  of  the  SSN  686  class  submarine,  ELSHOK 
was  used  to  predict  the  level  of  equipment  responses  pr i or 
to  the  tests  and  to  increase  the  level  of  confidence 
associated  with  the  tests.  This  code  has  now  been  made 
available  to  the  Naval  Postgraduate  School  by  the  Defense 
Nuclear  Agency  along  with  support  from  Weidlinger  Associates 


C.  PURPOSE  FOR  THIS  INVESTIGATION 

DDAM  is  a  design  tool.  Its  function  is  to  provide  a 
numerical  method  by  which  the  engineer  can  check  his  design 
as  to  its  adequacy  for  installation  in  a  submarine.  It  is 
intended  to  be  conservative  and  stra i gh tf orward  to  apply. 
When  DDAM  was  first  proposed  [Ref. 43,  its  authors  warned 
that  the  design  shock  factors  specified  [Ref. 2 3  were  not 
absolutely  determined  and  should  be  reviewed  as  explosive 
testing  progressed  and  data  was  accumulated.  Although 
testing  procedures  for  installed  equipment  components  have 
become  very  sophisticated,  these  simulators  may  not  accur¬ 
ately  reflect  the  response  of  installations  in  modern  day 
submarines.  Without  the  ability  to  carry  out  high- 
intensity  shock  tests  on  real  submarines,  verification  of 
the  design  shock  spectrum  used  in  DDAM  can  not  be 
accomplished  as  intended.  The  evolution  of  numerical  codes 
capable  of  simulating  the  response  of  installed  equipment 
to  a  specified  shock  loading  makes  available  a  means  to 
check  the  applicability  of  the  originally  proposed  design 
shock  factors  provided  in  DDAM  to  modern  submarines.  There 
is  no  flexibility  built  into  DDAM  to  allow  design  for  other 
than  the  originally  prescribed  shock  input  magnitude.  DDAM 
makes  no  separate  account  for  equipment/structure  inter¬ 


action  effects. 


The  purpose  of  this  investigation  is  to  examine  the 
response  of  several  equipment  models  installed  in  a  present 
day  scale,  6900  long  ton  <LT)  submarine  using  DDAM  and 
ELSHOK.  The  equipment  models  range  in  weight  from 
1,000-20,000  pounds.  Simple  models  are  used  to  make 
correlations  possible  between  the  two  methods  and  worst 
case  results  are  obtained.  Additionally,  equ i pmen t/hu  1  1 
interactions  are  investigated  for  their  possible  influence 
on  overall  equipment  response. 


I  I  .  DYNAMIC.  DESIGN  ANALYSIS  METHOD 


A.  DESCRIPTION  OF  DDAM 

As  men t  i  oned  in  the  introduction,  DDAM  is  a  simplified 
modal  analysis  method  which  utilizes  shock  inputs  which 
were  empirically  derived  -from  underwater  explosion  tests  of 
realistic  ship  and  submarine  installations.  Although  the 
analysis  is  identical  -for  both  types  o-f  vessels,  all  refer¬ 
ences  in  this  document  will  be  to  submarines.  In  order  to 
utilize  DDAM  to  evaluate  a  design,  several  assumptions  must 
be  sat i sf i ed . 

1.)  The  equipment  and  foundation  make  up  a  linearly 
elastic  system. 


2. )  The  structure  can  be  represented  reasonably 

well  using  a  lumped  parameter  model. 

3. )  The  structure  is  hard  mounted  to  the  rest  of 

the  submarine.  No  bottoming,  or  damped  mounts 
are  allowed. 

In  Navy  shock  requirements,  equipment  shock  classi¬ 
fications  are  based  on  a  graded  system.  These  grades  run 
from  A  for  mission  critical  items,  to  C  for  items  that 
could  indirectly  affect  the  ability  of  the  vessel  to 
function  as  intended.  For  instance,  a  locker  which  becomes 


adrift  and  could  injure  personnel  or  surrounding  equipment 
could  be  grade  C.  The  grade  of  the  item  is  specified  by 
contract  and  determines  the  degree  of  shock  qualification 
it  must  undergo  to  be  considered  safe  for  installation.  In 
this  study,  all  equipment  models  were  considered  grade  A. 

The  ideal  way  to  compute  the  elastic  shock  response  of 
an  installed  equipment  is  to  consider  the  entire  structure 
(submarine  and  equipment)  as  an  elastic  system.  The  normal 
modes  and  natural  frequencies  of  this  system  can  then  be 
determined  and  the  response  of  each  normal  mode  to  the 
water  applied  pressure  loadings  computed.  The  modal 
responses  can  then  be  superimposed  to  get  the  resultant 
equipment  response.  In  an  everyday  design  method  such  as 
DDAM,  this  approach  is  impractical.  To  simplify  the  above 
procedure,  DDAM  utilizes  the  concept  of  a  shock  spectrum. 
Any  point  on  the  complete  structure  can  be  designated  as  a 
reference  point  which  can  then  be  considered  to  be  a  fixed 
base  for  the  equipment  on  one  side  of  it.  The  dynamic 
response  of  this  portion  can  then  be  computed  by 
superposition  of  normal  mode  responses  to  the  motion  of  the 
fixed  base.  In  DDAM,  this  motion  is  specified  as  a  step 
input  to  each  normal  mode  and  is  calculated  by  empirical 
equations  based  on  underwater  explosion  tests  of  models  and 
submarines  [Ref.  53.  A  detailed  treatment  of  the 
mathematical  basis  for  DDAM  is  given  by  Butt  [Ref.  63. 


In  order  to  use  DDAM,  a  lumped  parameter  model  of  the 
equipment  must  be  -formulated  to  obtain  the  matrix 
expression: 

CM]  <X>  +  CK]  <X>  =  <0>  (1) 

which  has  the  associated  eigenvalue  problem 
CM]  CXja]  [  ]  =  [K]  CXja]  (2) 

where  CXja]  is.  the  matrix  of  eigenvectors  and  C  CO2  D  is  the 
diagonal  matrix  of  natural  frequencies.  The  advantage  of 
this  procedure  is  that  it  allows  uncoupling  of  the 
equations  of  motion  and  subsequent  solution  of  the  dynamic 
problem  in  terms  of  its  c omp on e n t  modes.  Careful 
formulation  of  the  equations  of  motion  to  capture  all  the 
dominant  characteristics  of  the  equipment  leads  to  good 
results  with  only  a  small  number  of  modes.  When  DDAM  was 
developed,  it  was  intended  to  be  used  primarily  in  hand 
computations.  The  major  difficulty  encountered  in  the  use 
of  DDAM  was  the  solution  of  the  eigenvalue  problem  which 
can  readily  be  handled  today  with  small  computers.  The 
normal  mode  theory  upon  which  DDAM  is  based  is  well 
established  in  shock  and  vibration  practice.  The  area  of 
concern  is  the  adequacy  of  the  adopted  design  inputs  to 
reflect  present  day  submarine  building  practices. 

B.  DESIGN  SHOCK  INPUTS 

The  specification  of  design  shock  inputs  for  use  in 
DDAM  calculations  determines  the  usefulness  of  this  method 


in  the  evaluation  of  equipment  installation  designs.  DDAM 
is  a  standard.  The  design  input  values  must  provide 
consistent  fixed-based  equipment  excitations  tor  the  broad 
base  of  submarine  platforms  in  use  today  to  simulate  a 
"standard"  underwater  explosion.  To  simulate  this  require¬ 
ment  here,  charge  weights  and  standoffs  were  selected  so  that 
each  case  corresponded  to  a  constant  value  of  energy  flux 
density;  i  .  e  .  , 

Ul/R2  =  Constant  <3> 

In  this  equation,  W  is  the  equivalent  weight  of  the  charge 
in  pounds  of  TNT  and  R  is  the  standoff  distance  from  the 
hull  in  feet.  The  submarine  platforms  in  use  today  range  in 
size  from  about  4,600-18,788  LT .  Each  submarine  hull  is 
characterized  by  its  own  modal  properties  (mass  distri¬ 
bution,  natural  frequencies,  mode  shapes).  The  same  piece 
of  equipment  as  is  installed  in  a  small  submarine  will 
demonstrate  a  different  response  when  installed  in  the  same 
configuration,  in  a  large  submarine.  The  effect  of  hull/ 
equipment  interaction  becomes  an  important  consideration  in 
large  pieces  of  equipment  tuned  to  one  or  more  natural  fre¬ 
quencies  of  the  hull  structure.  The  design  shock  input  used 
with  DDAM  must  incorporate  these  factors  so  that  when  it  is 
said  that  a  piece  of  equipment  meets  the  specification  for 
qualification,  the  qualification  level  is  the  same  for  all 
submarine  classes. 


In  the  publications  regarding  DDAM,  little  is  said 
about  the  origins  of  the  particular  shock  input  spectrum 
utilized  to  evaluate  equipment  response.  It  has  been 
empirically  derived  to  provide  values  consistent  with  the 
data  upon  which  it  is  based.  In  the  event  that  it  is 
decided  to  design  submarines  to  resist  underwater  explosion 
of  a  different  shock  intensity  than  the  one  chosen  for  the 
present  method,  a  major  effort  would  be  required  to 
construct  a  data  base  upon  which  to  base  the  new  input 
values. 

C.  USING  DDAM 

Since  the  development  of  DDAM,  new  tools  have  replaced 
the  slide  rules  of  engineers.  Among  them  are  the  readily 
available  desktop  computers.  In  the  course  of  this  work,  i 
was  necessary  to  use  DDAM  to  analyse  simple  structures  of 
fewer  than  56  degrees  of  freedom.  To  this  end,  a  computer- 
program  has  been  written  in  the  BASIC  computer  language  to 
carry  out  the  required  computations.  A  listing  and  short 
users'  manual  can  be  found  in  Appendix  A.  The  program  was 
verified  by  comparing  results  with  hand  calculations  and 
published  sample  problems. 


III.  EXPLOSIVE  SHOCK  RESPONSE  ANALYSIS  USING  ELSHOK 

A.  GENERAL  PRINCIPLES  OF  OPERATION 

The  ELSHOK  computer  code  [Ref.  43  calculates  the 
transient  response  of  a  submerged,  ring-stiffened  shell  of 
revolution,  with  or  without  internal  structure,  to  an 
underwater  shock  wave  emanating  from  an  explosive  source 
placed  at  an  arbitrary  location  away  from  the  shell.  The 
structure  is  considered  to  be  linearly  elastic,  and  the 
surrounding  fluid  is  treated  as  an  infinite  acoustic 
medium.  Modal  structural  analysis  is  used  in  all  phases  of 
the  calculations.  Internal  equipment  response  is  treated 
by  coupling  the  free-free  modes  of  the  empty  r i ng-st i f f ened 
shell  and  the  fixed-base  modes  of  the  internal  equipment 
through  use  of  dynamic  boundary  conditions  [Refs.  7,83. 

The  structure-fluid  interaction  is  approximated  by  means  of 
the  Doubly  Asymptotic  Approximation  (DAA>  due  to  Geers 
[Ref.  33.  The  form  of  the  DAA  used  is  that  obtained  when 
the  normal  fluid  displacement  of  the  structure-fluid 
interface  is  expanded  using  surface  expansion  functions 
which  are  orthogonal  over  the  wet  surface  of  the  submarine. 
At  frequencies  of  zero  and  infinity,  the  pressure-ve 1 oc i ty 
relations  are  exact  so  that  in  transient  analysis,  the  DAA 


yields  exact  results  at  early  and  late  times,  and  by  the 


nature  of  its  formulation  provides  a  smooth  transition 
between  these  two  limits.  In  effect,  the  DAA  allows  for 
uncoupling  of  the  fluid  field  problem  from  the  structural 
field  problem. 

The  structural  problem  solved  by  ELSHOK  is  separated 
into  two  parts.  A  modal  substructuring  procedure  is  used  to 
solve  the  dynamic  response  problem  for  internal  equipment. 
The  advantage  of  this  is  to  eliminate  the  need  to  handle  the 
modes  and  natural  frequencies  of  the  combined  structural 
problem,  as  well  as  the  requirement  for  a  combined  system 
stiffness  matrix.  Interaction  forces  and  moments,  and 
compatibility  of  deformation  at  the  she  11 -substr uc tur e 
attachment  points  are  used  to  solve  for  the  dynamic  response 
of  the  component  parts. 

Referring  to  the  ideal  case  of  DDAM  formulated  without 
the  use  of  a  shock  input  value;  ELSHOK  is  a  numerical  means 
to  arrive  at  the  input  to  the  substructure  without  depend¬ 
ence  on  an  explosive  testing  database  and  with  the  added 
advantage  that  interaction  effects  between  the  hull  and 
substructure  are  taken  into  account.  ELSHOK  performs  a 
transient  analysis  whereas  DDAM  utilizes  a  simplified 
"front  end"  to  arrive  at  the  maximum  forces  and  deflections 
in  a  given  response  problem  for  a  single  magnitude  of 


B.  ORGANIZATION  AND  IMPLEMENTATION  OF  ELSHQK 

The  ELSHOK  code  is  implemented  as  a  series  of  seven 
programs.  The  major  components  are: 

1.  B0S0R4  -  structural  analyzer  for  shell  [Ref.  93 

2.  ACESNID  -virtual  mass  processor 

3.  PI  FLASH  -  shell -fluid  processor 

4.  SAPIV  -  structural  analyzer  for  substructure 

[Ref.  103 

5.  PICRUST  -  substructure  processor 

6.  USLOB  -  time  integration  processor 

7.  PUSLOB  -  plotting  processor 

Figure  1  depicts  the  general  relationship  among  the 
code  modules.  The  numbers  above  correspond  to  the 
chronological  execution  order  to  be  followed  in  a  given 
analysis.  In  order  to  carry  out  an  analysis  using  ELSHOK, 
a  lumped  parameter  model  of  the  submarine  being  analyzed 
must  be  formulated  or  available  as  well  as  the  equipment 
model.  More  will  be  said  about  models  later.  The  compu¬ 
tations  are  carried  out  in  four  phases. 

B . 1 .  Phase  I  -  Shell  and  Fluid  Analysis 

B0S0R4  is  a  finite  difference  code  used  to 


capture  the  effects  of  the  submarine  hull  structure  on  the 
overall  response.  Two  models  are  actually  used:  a  full 
model  containing  information  about  the  major  configuration 


Figure  1  -  Organization  of  the  ELSHOK  Computer  Code  [Ref. 43 


of  the  submarine  and  its  general  dynamic  properties,  and  a 
compartment  model  representing  the  localized  area  of 
interest  on  the  hull.  The  -full  model  is  used  to  capture 
the  gross  effects  of  the  response  such  as  rigid  body 
translation  and  whipping  <N  =  1  modes),  and  rolling  and 
twisting  <N  =  0  torsional  modes)  of  the  entire  structure. 
The  compartment  model  is  used  to  capture  local  details  of 
the  response  in  the  area  of  interest  (N  =  0  breathing  and 
N>2  modes).  N  is  the  integral  number  of  circumferential 
wau e  s  or  harmon i c  s  in  the  c i r  c  umf  e  r  e  n  t i a  1  d i s  t  r i bu  t i on s  of 
the  surface  expansion  functions.  For  each  N,  a  separate 
B0S0R4  calculation  is  required.  The  BGS0F:4  code  provides 
the  in-vacuo  free-free  modes  and  natural  frequencies  of  the 
submar ine  hull. 

The  second  part  of  phase  I  is  the  use’  of  ACESNID 
to  compute  the  virtual  mass  array,  which  is  used  in  the 
late  time  contribution  of  the  DAA,  The  virtual  mass  array 
is  determined  from  the  solution,  based  on  simple  sources, 
of  a  steady-state  low  frequency  problem  in  which  normal 
displacements  corresponding  to  the  surface  expansion 
functions  are  appl ied  to  the  surface  of  a  cavity  of 
revolution  in  an  infinite  fluid  having  the  same  size  and 
shape  as  the  wet  surface  of  the  submarine.  Only  one  run  of 
this  module  is  required  for-  all  values  of  N  considered. 


PIFLASH  to  compile  the  data  -from  all  B0S0R4  runs  and 
ACESNID  into  a  she  1 1 -and-f 1 u i d  file  for  later  use.  The 
analysis  done  for  the  work  in  this  thesis  considered  only 
the  hull  response  of  one  submarine  so  only  one 
she  1 1 -and-f 1 u i d  file  was  needed  for  all  subsequent 
calculations. 

B . 2  Phase  II  -  Substructure  Analysis 

SAP IV  is  a  commonly  available  finite  element  code 
used  in  the  ELSHOK  code  to  perform  the  modal  analysis  of 
each  substructure  (equipment  model).  The  equipment 
installation  is  first  cast  into  an  appropriate  finite 
element  representation  and  then  SAPIV  is  used  to  solve  for 
the  desired  mode  shapes  and  frequencies.  This  information 
is  placed  in  a  substructure  mode  file  along  with  the 
geometry,  stiffness  and  mass  information.  The  mass,  mode, 
and  frequency  information  were  also  used  in  the  DDAM 
analysis  for  comipar  at  i  ve  purposes. 

The  final  step  in  phase  II  is  the  execution  of  the 
PICRUST  code.  This  code  prepares  the  input  for  phase  III 
calculations  by  placing  the  above  modal  information  on  a 
substructure  file.  Additionally,  it  calculates  influence 
coefficients  corresponding  to  the  base  or  support  motions  of 
the  substructure  and  accounts  for  the  user-specified 


con nee  t i v i ty . 


Following  the  completion  of  Phases  I  and  II,  the 
she  1 1 -and-f 1 u i d  file  is  merged  with  the  substructure  file 
to  form  the  input  file  for  USLOB  which  is  ELSHOK's  time 
.integration  processor.  USLOB  allows  specification  of  the 
charge  weight  and  geometry.  For  the  problems  of  concern 
here,  an  exponentially  decaying,  empirical  model  was  used 
to  represent  the  explosive  shock  source.  It  is  of  the 
form : 

P1(R)  =  Ki<Wl'®/R)K2  exp  <-t/60)  (4> 

where  ; 

PI<R)  =  the  incident  pressure  at  a  radius  R  from 
the  explosion  <on  the  hull) 

=  a  multiplicative  constant  for  incident 
pressure 

k2  =  a  spatial  decay  constant  for  incident 
pressure 

t  =  time  after  arrival  of  shock  wave  at  point 

of  interest 

W  =  weight  of  spherical  charge 

So  =  Ks  W1/s  <  W1/s  /  R  >K4 

K3  =  multiplicative  constant  for  time  constant 
of  exponential  decay 

K*  =  spatial  decay  constant  for  time  constant 
of  exponential  decay 

In  all  the  cases,  the  geometry  was  such  as  to  maintain  con¬ 
stant  energy  flux  according  to  eq.<3>.  The  integration  in 

time  is  done  using  a  modified  version  of  the  Runge-Kutta 


numerical  integration  scheme.  During  the  integration, 
USLOB  outputs  velocity  histories  -for  user-spec i f i ed  points 
on  the  hull  and  substructure. 


B . 4 .  PhaselV  -  Plotting 

Phase  IY  capabilities  were  not  utilized  in  this  inve 
ti gat  ion.  However,  ELSHQk  does  have  the  capability  to  gen¬ 
erate  plot  -files  containing  velocity  history  i  reformation  -for- 
processing  on  a  graphics  terminal. 

B .  5 .  De  -f  1  e  c  t  i  on  Calculations 

As  mentioned  previously,  ELSHOK  produces  velocity 
histories  at  the  points  of  interest  on  the  substructure. 

DDAM  produces  maximum  deflections  or  forces  at  the  same 
points  for  similar  models.  In  order  to  make  comparisons 
between  the  two  methods,  it  was  decided  to  convert  the 
velocity  histories  from  ELSHOK  to  deflection  histories  using 
an  integrator  based  on  Simpson's  1/3  Rule.  For  each 
equipment  case,  a  separate  program  was  written  in  BASIC  to 
accomodate  the  variations  in  configuration  between  each 
model  and  the  differences  in  output  from  the  USLOB  code  due 
to  point  of  interest  specification.  These  programs  are 
listed  in  Appendix  B. 


I U .  MODELS  USED  IN  THE  ANALYSIS 

A.  FINITE  ELEMENT/FINITE  DIFFERENCE  MODELS 

The  idea  of  representing  a  given  domain  as  a  col  lectio 
of  discrete  elements  or  values  was  not  first  applied  in 
finite  element/finite  difference  methods.  Early  mathemati 
cians  estimated  the  value  of  It  to  nearly  48  decimal  place 
by  representing  the  circle  with  a  polygon  of  a  finitely 
large  number  of  sides.  The  polygon  used  by  the  mathemati¬ 
cians  is  a  model  of  the  circle  they  wished  to  represent. 

In  modern  numerical  methods,  it  is  realized  that  very  few 
problems  in  the  real  world  can  be  anal yzed  wi thou t  the  use 
of  simplifying  models.  BQSQR4,  the  finite  difference  code 
in  ELSHOK  used  to  describe  the  dynamic  behavior  of  the  sub 
mar ine  hull,  requires  a  greatly  si mp 1  i f i ed  mode  1  of  its 
real  world  counterpart.  Similarly,  both  DDAM  and  ELSHOK 
use  a  substructure  model  to  capture  the  fixed-base 
properties  of  installed  equipment. 

B.  SUBMARINE  HULL  MODELING  USING  B0S0R4 

Actual  submarine  hulls  are  very  complex  structures. 
However,  they  do  possess  certain  characteristics  that  can 
be  used  to  simplify  their  modeling.  From  the  aspect  of 
submarine  hull  shock  response,  the  hull  is  a  free-free  r , n 
stiffened  cylinder  with  suitable  end  closures.  The 


internal  weights  of  machinery  and  other  out-fitting  must  be 
accounted  -for  but  these  items  can  be  smeared  into  their 
adjacent  hull  structure.  The  model  thus  becomes  a  cylinder 
whose  cross  sections  reflect  the  density  o-f  the  submarine 
at  a  given  f  ramie.  Since  a  cylinder  is  a  sur  -face  o-f 
revolution,  the  properties  of  each  cross  section  need  only 
be  specified  at  one  point.  The  reader  will  notice  that  any 
hard  structure  serving  to  stiffen  the  hull  in  a  non- 
ax  i  symmetric  sense,  such  as  non-circular  frames  or  massive 
machinery  foundations,  can  not  be  represented  in  this 
fashion.  If  these  structures  are  in  the  area  of  interest 
in  the  hull,  they  should  be  modeled  as  internal  structures 
in  phase  II  of  ELSHGK. 

The  B0SQR4  code  allows  the  submarine  hull  model  to  be 
constructed  by  specifying  the  properties  of  the  submarine 
hull  along  a  line  of  revolution.  A  numerical  comp  1  i c at i on 
enters  the  problem  when  dealing  with  models  of  submarines 
because  the  structure  often  repeats  itself  regularly  along 
its  axis  of  revolution.  With  such  a  configuration,  there 
are  many  modes  in  which  the  motion  of  the  stiffeners  is  of 
small  amplitude  compared  to  that  of  the  shell.  The  bays 
between  the  rings  vibrate  at  frequencies  which  may 
approximate  those  correspond! ng  to  a  s i mp 1 y-suppor ted 
cylinder  of  the  same  geometry  as  the  bay.  Multiple  or 
closely-spaced  eigenvalues  correspond  to  modes  in  which 
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one  or  more  of  the  bays  is  vibrating  while  others  are 
unaffected.  To  avoid  eigenvalue  difficulties,  the  prudent 
measure  is  to  analyse  as  small  a  segment  of  shell  as 
possible  in  order  to  avoid  numerical  difficulties  with 
multiple,  or  closely-spaced  eigenvalues  as  indicated  by 
Bushnell  [Ref.  ?] .  In  submarine  shock  analysis,  the 
interest  is  in  the  gross  response  contributions  of  the 
entire  hull  to  a  given  section  of  interest.  The  authors  of 
ELSHOK  have  allowed  for  this  through  the  use  of  two  models. 

A  full  model  is  used  to  obtain  the  gross  contributions 
of  the  entire  hull  structure  in  the  N  =  0(torsional>  and 
N  =  1  < tr ansi  at i onal  and  whipping)  modes.  This  will  include 
such  things  as  bulkheads,  heavy  ring  stiffeners,  and  end 
closures.  A  compartment  model  is  then  used  to  obtain  more 
detailed  information  about  the  particular  section  of 
interest.  In  this  way,  responses  away  from  the  section  of 
interest  do  not  have  to  be  carried  through  the  remaining 
calculations  since  they  do  not  affect  the  response  of  this 
area  anyway. 

Due  to  security  restrictions  placed  in  information  re¬ 
garding  the  construction  techniques  and  arrangements  of  U.S. 
nuclear  submarines,  it  was  not  possible  to  model  an  actual 
submarine  for  use  in  the  analysis  work  contained  in  this 
thesis.  With  the  aid  of  Weidlinger  Associates  of  New  York, 
however,  a  model  resembling  a  general  submarine  hull  struc¬ 
ture  of  <S?00  long  tons  displacement  was  obtained.  This 


30 


model  is  available  through  the  U.S.  Naval  Postgraduate 
School  but  is  not  included  here  so  that  restrictions  do  not 
have  to  be  placed  on  the  circulation  of  this  work. 

Figure  2  is  a  schematic  depiction  o-f  the  full  subma¬ 
rine.  The  submarine  being  modeled  has  been  cut  into  149 
slices  along  the  longitudinal  axis.  In  the  model,  infor¬ 
mation  regarding  Young's  modulus,  Poisson's  ratio,  mass 
density,  thermal  expansion  coefficient,  plating  half¬ 
thickness  and  whether  or  not  stiffeners  are  smeared  is  pro¬ 
vided  for  each  one  of  these  points.  Twenty-six  discrete 
rings  are  represented  along  the  axis  of  the  model,  twenty 
of  which  are  in  the  local  area  of  interest.  These  rings 
represent  frames  in  the  real  submarine  where  T-stiffeners 
give  increased  stiffness  to  the  hull.  In  actuality,  these 
stiffeners  are  spaced  fairly  evenly  along  the  entire  model, 
but  smearing  them  outside  the  area  of  interest  will  not 
appreciably  affect  the  vibrational  characteristics  in  the 
area  of  interest.  Using  the  B0S0R4  code,  modal  analysis 
using  finite  difference  techniques  is  carried  out.  Table 
1  is  a  listing  of  the  natural  frequencies  and  wave  numbers 
for  each  mode  retained  in  the  analysis.  Seven  N  =  0 
torsional  modes,  and  thirty-eight  N  =  1  rigid  body 
displacement  and  whipping  modes  have  been  retained  for  the 
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Submar i ne 


MODE  N  FREQ  MODE  N  FREQ 

<  Hz  )  ' Hz  > 


1 

1 

3.597568E-02 

3 

1 

2. 377235E+08 

5 

1 

1 . 29881 6E  +  @ 1 

7 

1 

1 .801 155E+0 1 

? 

1 

2 . 305725E+0 1 

1  1 

1 

2 . 992269E+0 1 

1  3 

1 

3.41 1510E+01 

15 

1 

3 . 563254E+0 1 

17 

1 

3.71 6095E+0 1 

19 

1 

3 . 847244E+0 1 

21 

1 

4 . 027563E+0 1 

23 

1 

4 . 295875E+0 1 

25 

1 

4 . 492591 E+0 1 

27 

1 

5.261 150E+01 

29 

1 

5 . 544796E+0 1 

31 

1 

5 . 994869E+0 1 

33 

1 

6 . 31 21 91 E+0 1 

35 

1 

6 . 5741 09E+0 1 

37 

1 

6.761 096E+0 1 

39 

0 

4.01 3475E-05 

41 

0 

1 . 465699E+0 1 

43 

0 

3.021 633E+0 1 

45 

0 

4 . 579604E+0 1 

2 

1 

4.71 6345E-02 

4 

1 

7.8901 12E+00 

6 

1 

1 .793631 E+0 1 

8 

1 

1 . 897726E+0 1 

1  0 

1 

2 . 7209 1 8E+0 1 

12 

1 

3 . 228725E+0 1 

14 

1 

3 . 498256E+0 1 

1  6 

1 

3 . 590094E+0 1 

18 

1 

3 . 306356E+  0 1 

20 

1 

3 . 962965E+0 1 

22 

1 

4 . 0591 42E+0 1 

24 

1 

4 . 390227E+0 1 

26 

1 

4.9161 27E+0 1 

28 

1 

5 . 289662E+0 1 

30 

1 

5 . 63470 1 E+0 1 

32 

1 

6.131 432E+0 1 

34 

1 

6 . 521 966E+0 1 

36 

1 

6.682896E+01 

38 

1 

6 . 345927E+0 1 

40 

0 

7 . 584749E+  00 

42 

0 

2 . 342790E+0 1 

44 

0 

3 . 855953E+0 1 

TABLE  I 


Wave  Numbers  and  Frequencies  of  Full  Model 


Figure  3  is  a  schematic  representation  ot  the  compart¬ 
ment  (local  area  ot  interest)  model  used  in  this  analysis. 
In  this  model,  all  -frames  are  modeled  as  discrete  ring 
st i ffeners.  Since  compatibility  must  be  maintained  between 
the  two  models,  (compartment  and  -full),  there  is  a  one  to 
one  correspondence  in  the  area  of  interest  between  points 
in  each  model.  The  comp ar tmen t  model  is  used  to  obtain  the 
modal  con tr i but i ons  corresponding  to  N  >  2  (normal,  merid¬ 
ional,  and  c  i  r  cum-feren  t  i  al  displacements).  As  mentioned 
previously,  deformations  or  motions  occur ing  beyond  the 
region  of  interest  should  be  unimportant  for  the  N  =  0 
breathing  modes  and  modes  corresponding  to  N  >  2.  Each 
bounding  ring  or  bulkhead  is  very  stiff  in  its  own  plane 
and  moves  primarily  in  that  plane  for  localized  loading. 
Table  2  is  a  listing  of  the  natural  frequencies  and 
c i rcumf er en t i a  1  wave  numbers  for  each  mode  retained  in  the 
analysis.  For  the  compartment  model,  25  modes  corres¬ 
ponding  to  N  =  2,  24  modes  corresponding  to  N  =  3,  24  modes 
corresponding  to  N  »  4  and  25  modes  corresponding  to  N  =  0 
have  been  retained.  A  total  0+  254  modes  were  found  by 
B0S0R4  for  the  full  and  compartment  models.  Of  these 
modes,  143  were  retained  for  further  calculations,  in  phase 
II  and  III  of  ELSHOK.  The  111  modes  that  were  dropped 
affect  responses  primarily  outside  of  the  area  of  interest. 
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Figure*  4  is  a  schematic  of  the  general  shell  structure 
mode  1  and  notation  conventions  used  in  the  remainder  of  th 
ELSHOK  analysis.  Two  coordinate  systems  are  used.  X  ,Y  an 
Z  refer  to  the  submarine  global  coordinate  -frame  and  x.y 
and  z  refer  to  the  substructure  local  coordinate  frame.  I 
both  systems,  the  >:-axes  run  from  aft  to  forward  and  the 
z-axis  orientation  is  related  to  the  Z  axis  in  the  global 
system  by  th*  angle  cc(j  where  0  refers  to  the  substructure. 
The  variable  s  is  used  to  locate  points  along  the  meridian 
of  the  full  model ,  and  u,v,  and  w  are  used  for  local  shell 
displacements  and  motions  located  c i rcumf er en t i al 1 y  by  the 
angle  3.  This  concludes  the  discussion  of  the  submarine 
model  used  in  ELSHOK. 

One  of  the  conveniences  realized  in  using  a  DDAM 
analysis  to  carry  out  design  checks  is  the  absence  of  a 
submarine  model  in  the  computations.  The  design  shock 
values  are  intended  to  provide  all  inputs  to  the  substruc¬ 
ture  model .  The  penalty  which  must  be  weighed  here  is  the 
lack  of  regard  which  this  places  on  the  pecul larities  of  a 
given  submarine  structure.  On  the  other  hand,  the  ELSHOK 
submarine  model  is  difficult  to  construct  requiring  a  grea 
deal  of  information  and  skill  to  properly  model  a  given 
hull,  but  the  price  is  paid  only  once  for  each  submarine 
class.  All  subsequent  calculations  in  this  thesis  are 
carried  out  utilizino  the  same  shell  model  calculations. 


C.  SUBSTRUCTURE  MODELING  USING  SAP IV 

In  both  DDAM  and  ELSHOK  calculations,  a  substructure 
model  is  -formulated  to  determine  the  modes,  masses,  and 
natural  frequencies  of  the  equipment  system  being  design 
checked.  When  DDAM  was  first  developed,  a  computer  program 
was  available  which  could  solve  the  eigenvalue  problem  for 
up  to  twenty  degrees  of  freedom.  However,  the  scarcity  of 
computer-  resources  dictated  that  most  users  utilize  hand 
computations  for  this  purpose.  Consequently,  most  early 
design  checks  were  limited  in  scope.  DDAM  is  well  suited 
to  tabular  methods  of  computation  but  with  the  common 
availability  of  smal 1  computers  today,  hand  computations 
are  no  longer  required,  and  the  restrictions  on  degrees  of 
freedom  are  largely  removed.  This  is  not  to  say  that  great 
numbers  of  degrees  of  freedom  are  required  in  every  case. 

By  careful  examination  of  the  equipment  installation  being 
checked,  the  major  response  contributions  can  be  captured 
using  a  small  number  of  modes.  In  the  cases  considered  in 
this  thesis,  emphasis  was  placed  on  using  equivalent  models 
for  DDAM  and  ELSHOK  rather  than  closely  model ing  a  real 
component  in  a  submarine  as  it  would  be  in  each  case  in 
actual  practice.  The  interest  here  lies  in  how  the  results 
compare,  for  the  same  models. 

Phase  II  of  ELSHOK  utilizes  the  SAP IV  finite  element 
code  to  model  submar i ne- i nstal 1 ed  equipment.  Because 


ELSHOK  is  written  utilizing  program  modules  which  execute 
independently  and  exchange  data  through  output  files,  the 
SAP IV  code  could  be  used  to  run  the  modal  analysis  of 


equipment  models  for  DDAM  and  ELSHOK  thus  ensuring  equival¬ 
ence  in  input  to  both  methods.  The  models  selected  for 
analysis  were  kept  sma! 1  so  that  compar i sons  could  be  made 
without  excessive  complication. 

Three  equipment  models  have  been  examined  for  this 
thesis.  Case  I  is  a  cantilever  foundation  mode  1  with  a 
1000  lb  weight  attached  to  its  free  end.  Cases  II  and  III 
are  a  simple  beam  foundation  supporting  two  weights  at 
center-span  and  having  attachment  points  on  two  separate, 
discretering  stiffeners. 

C . 1 .  Substructure  Case  I  Model 

The  case  I  model  development  is  depicted  in 
figure  5.  The  model  represents  a  1000  lb  valve  supported  by 
a  cantilever  foundation.  The  foundation  i s  bu i 1 t  up  from 
A I  SC  C  10x30  steel  channels  CRef.  113  and  the  valve  is 
rigidly  fixed  to  the  free  end.  Since  the  intent  was  not  to 
qualify  any  particular  valve  design,  the  entire  valve  has 
been  modeled  as  a  particle  (lumped  mass).  The  beam  mass  is 
represented  by  six  lumped  masses.  The  foundation  is 
designed  to  be  fixed  to  a  discrete-ring  stiffener  in  the 
submarine  in  an  upright  position.  Although  3-dimensional 
motoion  has  been  allowed,  the  maximum  deflection  occurs  in 
the  athwart  ship  direction  in  response  to  a  side -on  shock 


case 


loading.  This 


s  representative  o f  a  re! at i v e 1 y  1 


equipment  installation. 

C .  2 .  Substructure  Case  II  Model 

The  case  II  mode  1  development  is  depicted  in 
•figure  6.  The  model  represents  a  foundation  constructed 
from  AISC  W  27x177  I-beams  CRef.  Ill  supporting  two  10,080 
1b  weights.  The  foundation  structure  is  represented  by  IS 
beam  elements  and  1?  masses.  It  spans  two  discrete 
stiffeners  to  which  its  ends  are  fixed.  This  model  is 
representative  of  a  large  pump  or  turbomachinery 
installation, 

C .  3  .  Substructure  Case  III  Model 

The  case  III  model  is  similar  to  the  case  II 
miodel  except  the  foundation  structure  has  been  changed  to 
AISC  W  27x144  I-beams  CRef.  11]  and  the  supported  weights 
have  been  reduced  to  5,808  lb  apiece.  This  model 
represents  an  intermediate  weight  equipment  installation 
such  as  a  main  feed  pump. 


AISC  C  10X30  Steel 
Proper  t i es 


Foundat i on 


V.  ANALYSIS 


The  goal  of  the  analysis  procedure  used  in  this  inves¬ 
tigation  is  two-told.  The  primary  objective  is  to  analyze 
equivalent  substruc tures  using  DDAM  and  ELSHOK  so  a  com¬ 
parison  o-f  the  maximum  predicted  de-f  lections  can  be  made. 
Secondly,  any  hu 1 1 /substruc ture  interaction  effects  are  to 
be  noted.  The  analysis  is  complicated  by  several  factors. 
All  the  calculations  are  based  on  constant  energy  flux  which 
is  dependent  on  charge  weight  and  standoff.  Since  DDAM 
shock  inputs  ultimately  represent  the  results  of  explosive 
shock  tests,  the  geometry  of  the  “analytical  charge"  is 
invisible  to  the  user  and  assumed  in  the  empirical  shock 
spectrum.  ELSHOK  uses  inputs  of  charge  weight  and  standoff 
to  calculate  shock  loading  by  eq .  <4>  so  the  effects  of 
their  variation  will  change  the  transient  response  of  the 
substructure  even  though  the  shock  intensity  is  constant. 

In  this  analysis,  three  sets  of  charge  weight  and  standoff 
were  used.  A  further  complication  results  due  to  the  fact 
that  DDAM  only  can  be  used  to  calculate  maximum  relative 
deflections  or  forces.  ELSHOK  calculates  transient 
velocities  incurred  by  the  model  which  must  be  converted  to 
maximum  relative  deflections. 


A.  CASE  ANALYSIS  PROCEDURE 

Each  case  analysis  was  begun  by  -first  constructing  the 
equipment  models  to  be  investigated  using  SAPIV.  A  dynamic 
analysis  was  performed  to  find  the  natural  frequencies  and 
mode  shapes  of  the  model .  This  information  was  saved  along 
with  the  mass  matrix  for  subsequent  DDAM  calculations. 

Using  PICRUST,  the  SAP IV  data  was  then  reduced  to  a  suit¬ 
able  form  to  be  merged  with  the  hull  structure/fluid  data 
from  Phase  I  calculations.  Finally,  all  the  resulting  data 
were  integrated  using  the  USLOB  code  to  compute  the  tran¬ 
sient  velocity  profiles  for  each  equipment  installation 
configuration  and  charge  weight/standoff  set. 

The  PICRUST  code  was  used  to  specify  the  installation 
configuration  of  each  model.  Model  coordinate  system  orien¬ 
tation  to  the  hull  system  is  determined  by  the  angle  jn 
fig.  (4).  The  location  of  the  equipment  along  the  longi¬ 
tudinal  axis  is  also  specified.  For  example,  in  the  case  I 
model,  the  angle  is  279  degrees  and  the  base  attachment  is 
at  frame  95  in  the  submarine  model  (within  the  compartment 
model).  For  the  case  II  and  III  models,  two  different 
orientations  were  investigated,  one  with  the  model  mounted 
a thwart ship  and  the  other  with  the  model  mounted  verti¬ 
cally,  i.e.,  shock  input  from  the  side  and  bottom  of  the 


submar i ne . 


The  USLOB  code  allows  input  of  charge  weight  and 
standoff.  To  investigate  how  these  parameters  affect  the 
transient  response  of  the  substructure,  three  sets  ot  val¬ 
ues  were  used  tor  these  parameters  based  on  three  different 
charge  weights.  The  parameters  chosen  are  1 isted  in  the 
■following  table. 

Charge  Weight  Standoff  Distance 

(lbs  TNT)  < Inches/ Fee t > 

5,000  1,414/117.9 

10,000  2,000/166.7 

15,000  2,450/204.2 

TABLE  III  -  Charge  Weights  and  Standouts  -for  Analysis 

After  selection  of  all  factors  affecting  the  geometry 
of  the  problem,  the  time  step  increment  and  integration 
limits  were  specified.  In  each  case,  enough  time  steps 
were  chosen  to  capture  the  peak  response  amplitudes.  This 
number  was  found  by  tr i al  and  error . 

After  calculation  of  the  velocity  profiles  for  a  given 
model,  this  information  was  integrated  using  Simpson's  1/3 
rule  to  obtain  a  deflection  history  and  the  maximum  deflec¬ 
tion  response  of  the  model.  As  alluded  to  earlier,  the 
models  were  constructed  so  that  suitable  differences  taken 
between  the  velocity  profiles  of  designated  points  would 
yield  the  relative  deflection,  at  any  instant  in  time,  of 


the  point  in  q  u  e  s  t  i on . 


This  procedure  unfortunately  re¬ 
lies,  to  a  certain  extent,  on  symmetry  in  the  model .  For 
instance,  in  the  case  I  model,  the  de-flection  o-f  the  weight 
on  the  end  o-f  the  cantilever  in  the  athwartship  direction 
can  be  obtained  by  integrating  the  difference  between  the 
velocity  histories  o-f  the  base  and  the  tip  of  the  beam. 
Figures  <7>  and  <8)  show  schemat i cal  1 y  the  differences 
taken  to  calculate  deflections  for  each  model  configura¬ 
tion. 


^rel  -  ^base , athwar t sh i p  ~  Vt i p , athwar tsh i p 

Figure  7  -  Relative  Velocities  Used  to  Calculate 
Relative  Deflection  Between  Mass  and 
Base  for  Case  I 


Figure  8  -  Case  II  &  III  Relative  Velocity  Calculations  -for 
Athwartship  and  Vertical  Shock  Waves 


After  the  ELSHGK  results  were  obtained  for 


each  mode  1 

DDAM  was  used  to  calculate  the  maximum  deflections  using 
the  information  saved  from  the  SAP  IV  runs.  A  file  was  con 
structed,  using  an  editing  program,  to  input  the  masses, 
natural  frequencies,  and  rriode  shapes  for  the  mode!  to  the 
DDAM  program.  During  execution,  the  program  util ises 
prompts  to  query  the  user  about  the  type  of  installation 
and  shock  input  direction.  The  maximum  relative  deflect'd 
is  calculated  by  the  NRL  formula  [Ref.  123: 

^m  =  abs's  J  j  , maxi!  +  *J£ ( A  i  ’  £  —  J  j  , m a x *  (  5  ) 

where  ; 

abs<djjmaX>  =  the  largest  rriodal  deflection  of  a 
point  on  the  model 

•^i  =  modal  relative  deflection  at  a  point  on 
the  model 

*%■,  =  max  i  mum  relative  deflection  of  a  point  o 
the  model 

This  formula  is  used  to  reflect  the  fact  that  all  modes  do 
not  experience  their  max i mum  de 1 e  c  t i on  s  s i mu  1 t  an  e  ou  s 1 >  . 

Figures  <?>,  <10>,  and  <11>  summarize  the  model 
configurations  and  shock  input  directions  for  cases  I,  II, 
and  III.  Tables  IV,  V,  and  VI  are  summaries  of  the  in¬ 
formation  used  in  the  DDAM  and  ELSHOK  calculations  of  thes 


TABLE  IV  - 


I  Incut s  t  o  0 DAM  and  EL S H 0 K 


Number  of  Degrees  of  Freedom 

Number  o-f  Connected  Degrees  of  Freedom 

Number  o-f  Unconnected  Degrees  o-f  Freedom 


Number  o-f  Int< 


Points 


Number  o-f  Frequencies  In  Analysis 
Substructure  Natural  Frequencies  tops) 


1  445.41 

2  1054.85 

3  1343.07 

4  4837.85 


5  7608.71 

6  9553.85 

7  18454.9 


8  12169.6 

9  12284.8 

1  0  1 80 14.2 


Number  o-f  Time  Steps  In  the  Analysis- 
Size  o-f  Time  Steps 


0 . 5E-6 


TABLE  V  -  Case  II  Inputs  to  DDAM  and  ELSHGK 


Number  ot  Degrees  of  Freedom 

Number  o-f  Connected  Degrees  o-f  Freedom 

Number  o-f  Unconnected  Degrees  of  Freedom 

Number  of  Interface  Points 

Number  of  Frequencies  In  Analysis- 

Substructure  Natural  Frequencies  (cps) 

1  48.7743  4  551.143 

2  392.390  5  583.723 

3  392.508 


Number  of  Time  Steps  In  the  Analysis 
Size  of  Time  Steps 


Number  of  Degrees  of  Freedom 

Number  of  Connected  Degrees  of  Freedom 

Number  of  Unconnected  Degrees  of  Freedom 

Number  of  Interface  Points 

Number  of  Frequencies  In  Analysis 

Substructure  Natural  Frequencies  (cps) 

1  52.7625  5  570.582  8  1213.51 

2  388.0991  6  353.971  9  1433.73 

3  381.061  7  1213.45  10  1492.67 

4  537.903 

Number  of  Time  Steps  In  the  Analysis 
Size  of  Time  Steps  2 . 


VI .  RESULTS 


Although  the  number  of  cases  and  their  permutations 
explored  in  this  investigation  are  not  numerous,  they  have 
been  contrived  to  accentuate  the  similarities  and  differ- 
ences  between  the  DDAM  and  EL3H0K  methods  of  analysis.  As 
mentioned  previously,  ELSHOK  provides  a  great  deal  more 
in-formation  about  the  shock  response  of  submar  i  ne- i  nstal  1  ed 
equipment  than  DDAM.  However,  DDAM  is  not  intended  to  be  a 
theoretical  tool  but  rather  is  simply  a  design  check.  For 
this  reason,  where  ELSHOK  will  predict  different  maximum 
deflections  depending  on  charge  weight  and  standoff,  DDAM 
will  predict  only  one  value  for  all  variations. 

A.  CASE  I  <  1,000  LB  CASE > 

The  predicted  responses  of  DDAM  and  ELSHOK  compare  well 
for  this  case.  The  DDAM  calculation  yields  a  maximum  de¬ 
flection  of  0.0143"  for  the  end  of  the  cantilever  beam  when 
subjected  to  an  athwartship  shock  input,  whereas  the  ELSHOK 
calculation  predicted  a  0.0131"  deflection  for  a  10,000  lb 
charge  and  a  2,000  inch  standoff  as  in  fig  <?b>.  The 
15,000  lb  and  5,000  lb  charges  both  caused  slightly  lower 
deflections.  In  this  case,  DDAM  is  more  conservative 
(about  ?'/.')  than  the  ELSHOK  calculation.  This  illustrates 


that  where  the  equipment  weight  is  only  about  1,000  lb, 


little  interaction  is  to  be  expected  between  the  hull 
response  which  is  relatively  low  -frequency,  and  the  higher 
•frequency  equipment  response.  The  highest  hull  frequence 
is  only  about  100  Hz.  Coupt  ing  effects  between  the 
responses  o-f  the  equipment  mass  and  the  hull  mass  are 
i  ns  i  gn  i  f  i  gan  t  due  to  the  large  ratio  o-f  the  masses.  If.  the 
mass  of  only  segment  five  in  the  full  model  is  considered 
to  be  the  entire  hull  mass,  then  the  effective  mass  of  the 
hull  is  about  4100  slugs.  The  mass  of  the  model  is  less 
than  2.9  slugs.  This  yields  a  mass  ratio  of  about  0.00073. 
DDAM  is  based  upon  the  approximation  that  the  hull  mass  can 
be  considered  infinite  which  implies  a  mass  ratio  of 
8.00000...  providing  no  account  for  coupling  other-  than 
through  the  safety  factor  built  into  its  design  spectrum. 
The  close  agreement  between  ELSHOK  and  DDAM  is  probably  due 
to  the  lack  of.  response  coupling  in  both  methods  of 
anal y  s i s • 

B.  CASE  II  (20,000  LB  CASE) 

Two  equipment  configurations  were  investigated  in  the 
ELSHOK  analyses  for  this  case.  In  the  problem  geometry  in 
which  the  submarine  is  side  loaded,  figure  (9c),  the 
maximum  deflection  is  again  in  good  agreement  with  that 
predicted  by  DDAM;  0.2433  inches  in  the  former  versus 
0.2744  in  the  latter.  However,  the  input  motion  to  the 
structure  is  not  primarily  in  the  direction  of  rrpaximutTi 
deflection  in  this  instance.  In  the  second  configuration 


o-f  this  case,  corresponding  to  -figure  <9d),  a  much  greater 
de-flection  <  1 33V.  greater)  is  predicted  by  ELSHOK.  There 
are  two  p  oss i b 1 e  e  x  p 1  an  a  t i on  s .  Some  hull /s  t ructure 
interaction  is  to  be  expected  because  the  fundamental 
natural  frequency  of  the  substructure  is  43.3  Hz  which 
falls  in  the  middle  of  the  entire  range  of  hull  frequen¬ 
cies.  If,  as  in  the  former  case,  a  hull  effective  mass  of  • 
4100  slugs  is  assumed,  and  the  model  mass  is  67.4  slugs  in 
this  case,  then  the  ratio  of  the  masses  is  increased  to 
0.0164  or  about  1.6  percent.  In  this  example,  the  influ¬ 
ence  of  coupling  on  the  response  can  be  expected  to  amp  1  ify 
the  expected  result  for  an  infinite  mass  hull  system 
considerably.  Also  a  source  of  uncertainty  is  the  test  data 
upon  which  the  shock  input  spectrum  used  in  DDAM  is  based. 


Most  explosive  tests  are  set  up  with  the  charge  geometry  of 
figure  ( 9c >  to  avoid  problems  associated  with  bubble  migra¬ 
tion  to  the  hull  due  to  bouyancy  effects.  In  the  second 


configuration,  figure  <9d),  the  charge  was  situated  as  if 
it  had  detonated  beneath  the  keel  of  the  submarine,  a 
realistic  possibility  when  considering  modern  weapons 
tactics  but  perhaps  not  simulated  well  by  explosive  test¬ 
ing.  This  particular  configuration  subjects  the  structure 
to  greater  shock  input  motions  in  its  more  flexible  direc¬ 
tion  producing  the  larger  deflections,  however  DDAM 
provides  the  same  inputs  for  both  problems. 


--  ^  „*■  o  .*■  .*■  **•  .* 


37 


c. 


CASE  III  (18,000  LB  CASE' 


Two  equipment  configurations  were  investigated  in  the 
ELSHOK  analysis  for  this  case.  The  results  are  very 
similar  to  the  results  for  case  II.  In  the  configuration 
of  Fig.  10(a),  a  maximum  deflection  of  0.2184  inch  was 
predicted  and  in  the  configuration  of  Fig  10(b),  a -maxi  mum 
deflection  of  0.892  inch  was  predicted  by  ELSHOK  whereas 
DOAM  predicted  a  maximum  of  0.2743  inch  for  both  config¬ 
urations.  Again,  the  same  possible  causes  can  be  cited  for 
the  large  <150%>  difference  in  the  second  instance.  In 
this  case  the  fundamental  frequency  of  the  foundation  is 
52.7  Hz.  The  mass  of  the  equipment  model  in  this  case  is 
38.0  slugs  yielding  a  mass  ratio  between  the  hull  and 
substructure  of  about  0.00877. 

A  summary  of  the  calculated  results  is  presented  in 
Table  VII.  With  the  limited  calculations  available  in  this 
work,  no  attempt  could  be  made  to  deduce  systematically  the 
exact  cause  for  discrepancies  between  ELSHOK  and  DOAM  in 


cases  II  and  III. 


LE  VII  SUMMARY  OF  RESULTS 


■v'll  CONCLUSIONS 


Since  its  adoption  by  the  naval  shock  community,  DDAM 
has  provided  an  easy-to-use  and  convenient  means  to  design 
check  equipment  installations  proposed  -for  submarines.  It 
is  independent  of  any  well  de-fined  submarine  structural 
input  and  can  be  applied  early  in  the  overall  design  pro¬ 
cess  before  such  details  are  well  developed.  However,  the 
same  qualities  that  allow  tor  the  -flexibility  in  DDAM  may 
also  contribute  to  possible  inaccuracies  when  this  method 
is  applied  to  submarines  o-f  radically  di-f-ferent  size  or 
design  than  the  ones  reflected  in  the  empirical  Design 
Spectrum . 

The  work  in  this  thesis  is  in  no  way  construed  to  be 
an  all-encompassing  evaluation  of  DDAM  and  ELSHQK.  From 
the  infinite  number  of  possible  equipment  configurations 
and  sizes,  three  simple  models  have  been  selected  which  ex¬ 
amine  only  the  hull  mounted  equipment  problem.  The  intent 
of  this  investigation  has  been  to  examine  the  capability  of 
the  two  methods  to  predict  the  shock  response  of  these 
equipments  and  not  to  rank  one  against  the  other.  ELSHQK 
or  any  variation  of  this  code  is  much  too  complex  and 
computer  resource  intensive  to  use  as  a  design  tool  and  was 
not  intended  for  this  purpose. 


In  the  case  of  equipment  whose  response  is  relatively 
high  -frequency  compared  to  the  hull  response,  DDAM  and 
ELSHOK  deflection  predictions  are  in  good  agreement.  The 
maximum  hull  response  frequency  is  limited  to  about  100  Hz 
for  the  6900  LT  submarine  used  in  this  work.  The  good  cor¬ 
relation  is  attributed  to  the  lack  of  interaction  and 
coupling  effects  between  the  hull  and  substructure.  When 
the  response  frequencies  of  the  substructure  approached 
those  of  the  hull  and  the  mass  of  the  substructure  became  a 
significant  portion  of  the  overall  system,  the  comparisons 
between  ELSHOK  and  DDAM  deflection  predictions  were  not 
good . 

The  indications  of  this  work  are  that  DDAM  does  not 
correctly  reflect  the  hu 1 1 -subs true ture  response  coupling 
amplification  which  becomes  apparent  when  the  substructure 
response  is  tuned  to  the  hull  response.  Additionally,  the 
empirical  database  upon  which  DDAM  is  reliant  for  shock 
input  may  fail  to  represent  the  heavy  equipment  problem 
well. 

More  basic  research  remains  to  be  done  before  many  of 
the  questions  or  criticisms  posed  by  this  thesis  can  be 
considered  conclusive.  Due  to  the  scarcity  of  published 
reports  describing  the  input  data  used  to  generate  the 
design  shock  values  used  by  DDAM,  the  application  of  this 
method  as  a  general  design  qualification  should  be  regarded 


with  caution.  The  basis  of  the  Design  Spectrum  used  in 
DDAM  should  be  available  -for  analysts  and  revision. 

Additional  work  remains  to  be  done  to  determine  where 
the  discrepancies  between  ELSKGK  and  DDAM  first  become  sig¬ 
nificant.  The  cases  analyzed  here  represent  both  ends  of 
the  substructure  .response  spectrum  but  lend  little  infor¬ 
mation  to  what  happens  in  the  transition  region  between  the 
1 ow- i n t er ac t i on  and  h i gh- i n terac t i on  frequency  response 
regimes. 

DDAM  should  be  retained  as  a  design  check  method. 
However,  some  modifications  are  recommended  here.  The 
present  Design  Shock  inputs  used  in  DDAM  should  be  updated 
to  reflect  the  i ncresed  size  of  modern  day  submarines.  If 
the  details  of  the  original  f or mu  1  at i on  were  known  or 
revealed,  the  current  values  could  be  updated  by  generating 
new  analytical  data  based  on  computer  simulations.  This 
new  spectrum  could  then  be  verified  by  explosive  testing. 

It  is  apparent  from  the  hu 1 1 -subs true ture  coupling  phe¬ 
nomenon  noted  here  that  a  modified  shock  spectrum  should  be 
formulated  for  each  new  class  of  submarine  which  would 
reflect  the  peculiarities  of  this  class  versus  the  previous 
ones . 

At  the  present  time,  when  new  submarines  cost  billions 
of  dollars  and  are  by  no  means  numerous,  every  effort  must 
be  made  to  ensure  their  survival  in  an  underwater  explosive 


shock  environment.  This  should  include  taking  a  new  look 
at  DDAM  and  ensuring  it  upholds  the  standards  set  by  nave, 
shock  policy.  There  is  a  great  deal  of  -flexibility  inher¬ 
ent  in  DDAM  and  together  with  modern  numerical  techniques 
it  can  be  updated  to  reflect  current  technology. 
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The  DDAM  Program  contained  herein  provides  an  easy 
means  -to  perform  a  small  to  medium-scale  DDAM  analysis  of 
given  internal  equipment  design.  The  program  is  written 
IBM  PC  BASIC  language  and  is  not  restricted  to  submarine 
analyses  only  but  includes  the  standard  DDAM  case 
p oss i b i 1  i t i e s .  The  solution  a  1 gor i t hm  foil ows  the  me  t h od 
cited  in  Ref .  2 . 

A.  USING  THE  DDAM  PROGRAM 
A . 1 .  Option  Selections 

The  DDAM  program  is  menu-driven  and  written  to  b 
self-explanatory.  Prior  to  starting  the  program,  it  must 
be  av a i 1 ab 1 e  on  a  diskette  along  with  the  BAS I CA . COM 
program . 

To  start  the  program  type: 

BAS I CA< ENTER) 

< F3> DDAM < ENTER > 

<F2> 

The  first  command  1  ine  loads  and  executes  the  BASIC  inter¬ 
preter  and  the  second  line  loads  the  DDAM  program.  The 
third  line  starts  execution  of  the  DDAM  Program. 

Upon  starting,  the  title  of  the  program  <  DDAM >  i 
written  on  the  screen  and  several  seconds  later  a  message 


I 


b 

P 


appears  describing  the  provided  functions.  Follow  the 
d i r ec t i ons . 

Four  execution  options  will  appear  for  source  of 
data  and  program  exit  to  the  operating  system.  They  are: 

[13  Input  Mode  Shapes  and  Mass  Matrix  for  the 
problem  from  Keyboard 

C23  Input  Mode  Shapes  and  Mass  Matrix  from  a 
disk  file 

[33  Input  a  Mass  Matrix  and  Stiffness  Matrix 

for  previously  formulated  equations  of  motion 

[53  Exit  to  DOS 

No  op  t i on  [43  is  p  r  e  sen  t ;  h  owe  v  e  r ,  p  r  ov i s i on  s  have  bee  n 
made  to  install  a  user-defined  subroutine  at  1  ine  25009  in 
the  program. 

Option  [13  is  useful  for  small  problems  and  is  not 
recommended  for  situations  involving  more  than  four  degrees 
of  freedom.  It  allows  all  data  input  to  be  carried  out 
using  the  keyboard.  In  larger  problems,  the  entry  of  mode 
shapes  becomes  tedious  and  error  prone. 

Option  [23  is  recommended  for  most  problems.  It 
allows  the  use  of  a  diskfile  for  input.  The  file  is 
user -spec i f i ed  and  utilizes  a  free-format  data  structure. 
The  information  must  be  present  in  the  file  in  a  set  order 
which  is  described  later. 


6? 


Selection  of  option  C31  causes  a  message  to  be 
displayed  stating  that  the  JACOBI  program  must  -first  be 
run  to  solve  the  eigenvalue  problem  for  mode  shapes  and 


natural  -frequencies.  This  program  listing  -follows  the  DDAti 
listing  in  this  appendix.  Incorporating  the  JACOBI  program 
in  DDAM  would  have  reouced  significantly  the  size  of  pro¬ 
blem  which  could  be  handled.  The  JACOBI  program  is  in¬ 
cluded  as  a  convenience  to  the  reader  and  was  not  integral 
to  the  analysis  of  the  case  studies  done  in  this  thesis. 

Its  operation  is  described  later. 

A . 2 .  QptionCl]  Execution 

Selecting  option  Cl]  allows  input  of  the  masses, 
modes,  and  natural  frequencies  of  the  problem  from  the 
keyboard.  The  number  of  modes  and  masses  are  first  input. 

A  mass  input  subroutine  queries  the  user  interactively  for 
the  diagonal  entries  of  the  mass  matrix.  After  input  is 
completed,  the  diagonal  elements  are  displayed  enabling  the 
user  to  make  corrections  as  necessary.  A  zero  correction 
terminates  mass  matrix  input.  Next,  the  program  asks  for 
input  of  the  natural  frequencies  of  the  system.  It  re¬ 
quires  one  natural  frequency  for  each  mode  shape.  The 
input  format  is  similar  to  the  one  used  to  input  the  mass- 
elements.  The  natural  frequencies  are  required  in  units  of 
radians/second.  Each  mode  shape  is  then  input.  After 
initializing  each  mode,  a  correction  option  is  provided. 


The  program  next  calculates  the  participation  -factors  and 
mods  1  we  i  gh  t  s  -f  or  t  h e  p r  ob  1  em  . 


Following  these  preliminary  steps,  DDAM  requires 
information  to  fix  the  installation  details  for  the  equip¬ 
ment  and  type  of  analysis  to  be  carried  out.  The  standard 
cases  of  Ref.  2  are  allowed.  Once  the  case  is  fixed,  the 
Design  Values  are  calculated  and  displayed  as  DA< 1 ) . . . DA< n ) 
where  n  is  the  number  of  modes  selected  for  the  analysis. 
The  Design  Values  are  the  equivalent  static  accelerations 
calculated  from  empirical  formulas  for  each  mode . 

Several  output  options  are  provided.  Modal  forces 
or-  deflections  or-  both  can  be  selected.  After  selection, 
the  appropriate  calculations  are  made  and  the  results  dis¬ 
played.  Initially  a  print  option  was  separately  provided. 
However,  in  the  present  version,  selection  of  the  display 
to  screen  option  with  printer  on  will  cause  output  to  go  to 
the  screen  and  the  printing  device.  Disabling  the  printer 
will  cause  output  to  go  only  to  the  screen.  Following 
the  output  options,  choices  are  provided  to  change  the 
equipment  installation  type  and  repeat  the  analysis  or  to 
exit  the  program . 

A . 3 .  O&tion  C  2  3  Execution 

Selecting  option  C  2  3  allows  input  of  masses,  nat¬ 
ural  frequencies,  and  mode  shapes  from  a  disk  file.  After 
the  data  is  read  in,  the  program  execution  is  similar  to 
that  in  option  Cl].  Option  C  2  3  prompts  the  user  for  the 


dimension  of  the  mass  matrix,  the  number  of  modes  in  the 
analysis,  and  the  name  of  the  input  file.  The  input  file 
is  a  free-format  ASCII  file  containing  the  masses,  natural 
frequencies,  and  mode  shapes  for  the  problem.  This  file 
may  be  constructed  from  structural  analyzer-  (SAPIU)  output 
or  any  other  problem  formulation  method  using  a  suitable 
editor  program.  (EDLIN,  the  editor  provided  with  the  “IBM 
Disk  Operating  System"  is  not  recommended  for  this  purpose 
due  to  its  line  orientation.)  The  format  for  the  input 
file  foil ows : 

mass( 1 )  ....  mass(m)  (slugs) 

omega(l)  ....  omega <m)  (rad/s) 

p h i < 1  , 1 )  ....  p h i ( 1  , n ) 

ph i (m, 1 )  ....  ph i (m,n) 

where  m  =  number  of  masses  and 

n  =  number  of  modes  in  the  analysis 

In  the  above  format,  the  data  is  read  in  sequentially  from 

left  to  right.  The  data  is  echoed  to  the  screen  and/or 

printer  in  order  to  provide  a  check  for  proper  input 

operation. 


8.  PROGRAM  LISTING  -OOAM 


18  CLS:  KEY  OFF  :  LOCATE  8,1,8 

28  PRINT" ********************************************************************** 

21  PRINT"** 

**• 


38  PRINT"** 

DDDD 

DDDD 

m 

m 

tfff 

**" 

40  PRINT"** 

DO  D 

00  0 

A 

A 

Ml 

m  m 

**• 

58  PRINT"** 

DO  D 

DD  D 

rfWWVWi 

Iff 

M  Iff 

**■ 

48  PRINT"** 

DO  0 

DO  0 

A 

A 

m 

Iff 

**" 

70  PRINT"** 

DDDD 

DDDD 

A 

A 

Iff 

**■ 

88  PRINT"** 

**" 

90  PRINT"** 

DYNAMIC 

DESIGN 

ANALYSIS 

METHOD 

**" 

188  PRINT"** 

**" 

118  PRINT" »* 

**■ 

1 28  PRINT" ********************************************************************** 
***** 

121  DIM  MASS! 48) 

122  DIM  PHI (48, 40) 

123  DIM  W(48) 

124  DIM  MX(48,40),MX2(40,48),SUll(40),SIlM2(48)tP<48) 

125  DIM  UA(48) 

127  DIM  V0T48) ,A8(48) ,DA(48) 

128  DIM  AA(40) ,VA(48) ,F(48,48) ,0EFL<48) ,DELTA<40 ,48) 

138  FOR  I  =  1  TO  5888  :  NEXT  I 

140  CLS:  LOCATE  1,1,0 

158  PRINT"  The  purpose  of  this  progran  is  to  allow  the  user  a  convenient  raeth 
od":PRINT 

140  PRINT'by  which  to  carry  out  a  snail  to  reediun  sized  DOAM  analysis  of  a  give 
n":PRINT 

170  PRINT'piece  of  equipment,  Fanil iarity  with  NORMAL  MODE  ANALYSIS  and  sene  u 
n-":PRM 

188  PRINT'derstanding  of  the  theory  behind  00AM  is  prerequisite  to  confident  us 
e":PRINT 

198  PRINT'of  this  progran.  The  program  is  prompt  driven  but  the  user  should  hav 
e  " :PRINT 

280  PRINT" avail  able  to  him  /  her  either  :':PRINT:PRINT 

218  PRINT"  A.)  Mode  Shapes  and  a  Mass  Matrix  for  the" 

220  PRINT"  problem  of  interest." :PRINT:PRINT 

230  PRINT"  8.)  A  Mass  Matrix  and  Stiffness  Matrix  from" 


»  '  i  *  *  ^  »"*  «  •  , 
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241  PRINT*  previously  formulated  EQUATIONS  of  MOTION* ;PRM 

:PRINT:  PR  I  NT 

251  PRINT*  (  Press  ANY  key  to  continue  )':LQCATE  25,1,8 

268  At  «  INKEY* :  IF  A*  =  **  THEN  269 
279  CLS: 

288  LOCATE  5,1,8 

298  PRINT*  In  order  to  initialize  the  program  input,  choose  one  of  the  follow 
ing* :PRINT 

388  PRlNT'menu  selections  = ==>* ‘.PRINT sPRINT 

318  LOCATE  18,28,8  sPRINT’tlJ  Input  Node  Shapes  and  Mass  Matrix* 

328  LOCATE  11,28,8  sPRINT*  for  the  problem  from  keyboard* 

325  LOCATE  13,28,8  :PR1NT*I21  Input  Mode  Shapes  and  Mass  Matrix* 

328  LOCATE  14,28,8  sPRINT*  fron  a  disk  file  * 

338  LOCATE  16,28,8  sPRINT* C33  Input  a  Mass  Matrix  and  Stiffness  Matrix* 

348  LOCATE  17,28,8  sPRINT*  for  previously  formulated  equations  of* 

358  LOCATE  18,28,8  sPRINT*  motion* 

368  LOCATE  28,28,8  sPRINT* [51  Exit  to  DOS* 

378  A*  =  1NKEYI  :  IF  A*  )=  CHR$(49)  AND  A*  <=  CHR$(53;  GOTO  488 

388  IF  A*  =  **  GOTO  378 

398  BEEP  :  SOTO  378 

488  8  =  VAL(A$)  :  CLS 

418  ON  8  GOSUB  5888,15888,28888,25888,588 

588  SYSTEM  :0fD 

688  'HIHHHHHHHHHHHIHHIiHHIHHiiHimmiHiHHHHHilllilHH 

681'  SUBROUTINE  TO  INPUT  MASS  MATRIX 

683  'HHHHiHHtillHHiHHHlHHtHHiHliiiHliHiHHIiliHHHHiHHi 

685  INPUT  ’Enter  the  dimension  of  the  Mass  Matrix  — >  *;  N 
618  IF  N  <=  28  GOTO  622 

628  PRINT  sPRINT  sPRINT*  PROBLEN  SI2E  IS  LIMITED  TO  28  X  28  MATRICES  mm 
*  *:PRINT:PRINT:  GOTO  688 

621  PRINT  sPRINT sPRINT 

622  INPUT  'Enter  the  number  of  modes  you  wish  to  use  =)  *jM 

638  CLS:  PRINT*  »mm  Enter  the  elements  of  the  DIAGONAL  Mass  Matrix  mm**jPR1 
NT sPRINT 

648  PRINT*  If  you  make  an  error  you  will  have  an  opportun-’ 

658  PRINT*  ity  later  to  correct  it.* sPRINT sPRINT 

678  FOR  I  «  I  TO  N 

688  PRINT  *  M<*jli*,*sl;*>  =  *;:1NPUT*  *,  MASS(I) 

698  IF  <1  =  5)  OR  (I  =  18)  OR  <1  =  15)  THEN  CLS 
788  IF  MASS(I)  0  8  GOTO  738 

718  PRINT*  mm*  MASS  MATRIX  MUST  BE  POSITIVE  DEFINITE  —  INPUT  NON-ZERO  VALU 
E  MM* 

728  PRINT:  GOTO  688 
738  PRINT:  NEXT  I 
748  CLS 

758  FOR  1  =  1  TO  N 

768  PRINT  *N(*:I {■,*;!;*)  =  *;MASS<3), 

778  NEXT  I 

788  PRINTsPRINT sPRINT 

798  PRlNT’Enter  diagonal  element  number  to  correct  mistaken  entries  or  * 

888  PRINT: INPUT  *  Enter  <  ZERO  >  to  continue  =>*,  RESPONSE 


818  IF  RESPONSE  =  8  THEN  RETURN 

828  JF  RESPONSE  <=  9  OR  RESPONSE  >  N  THEN  CLS  :  SOTO  758  ELSE  I  =  RESPONSE 
858  PRINT; PRINT: INPUT  'Enter  correct  element  value  =>  *,  M<1) 

868  PR1NT:PRINT'M(*;I;*,';I ;')  CfWfGED  TO  *ilttSS (1) 

878  FOR  1  =  1  TO  1588  :  NEXT  I  :  CLS  :  GOTO  758 
888  RETURN 

1888  'HeHHiiHfHHIlHHHHMHMHHHifHMlHiHfilHHiHHHHHHHi 

1818  '  SUBROUTINE  TO  INPUT  MODE  SHAPES 

1821  '*********«*************************#************************************ 
1858  CLS:  LOCATE  5,1,8 

1868  PRINT1*****  Enter  OMEGA'S  corresponding  to  Mode  Shapes  ****** :PRINT 
1878  PRINT'  <  smallest  to  largest  )': PRINT sPRINT 

1888  FOR  I  =  1  TO  N 

1898  PRINT  '  onega(';I;*)  =  ' ; •.INPUT''  *,«(!) 

1188  IF  (1=5)  OR  (1=18)  OR  (1=15)  THEN  CLS 
1118  NEXT  1 

1128  CLS  :  FOR  1  =  1  TO  M 

1138  PRINT*0MEGA(';1;‘)  =  " ;W( 1 ) , 

1148  NEXT  1 

1158  PRINT :PRINT :PRINT 

1168  PRINT  'Enter  frequency  number  to  correct  mistaken  entries  or  *:PRINT 
1178  INPUT  '  Enter  <  ZERO  )  to  continue  =>\ RESPONSE 

1188  IF  RESPONSE  =  8  THEN  GOTO  1238 

1198  IF  RESPONSE  (=  8  OR  RESPONSE  )  N  THEN  CLS  :  GOTO  1128  ELSE  I  =  RESPONSE 
1288  PRINT :PRINT: INPUT  'Enter  correct  frequency  value  =>  ',  U(I> 

1218  PRINT •.PRINT'QMEGACjl;*)  changed  to  ,}U<1) 

1228  FOR  I  =  1  TO  2888  :  NEXT  I  :  GOTO  1128 
1238  CLS;  LOCATE  5,1,8 

1248  PRINT'*****  Enter  ElGEWECTQR's  corresponding  to  Mode  Shapes  ****** ’.PRINT 
1258  PRINT*  <  begin  with  Nod#  One  )*:PRINT:PRINT 

1268  FOR  I  =  1  TO  M 
1278  FOR  J  =  1  TO  N 

1288  PRINT*  MODE  SHAPE  ';!;*  PhiCjJj')  =  *  j :  INPUT  *  *,PHI(J,J) 

1298  IF  (J=5)  OR  (J=18)  OR  <J=15>  THEN  CLS  :  LOCATE  5,1,8 

1380  NEXT  J 
1305  CLS  ;  PR1NT:PRINT 

1318  PRINT  ****************  MODE  NUMBER  * ;  I  j*  #**************’ 

1320  PRINT  :PRINT  :FQR  J  =  1  TO  N 

1338  PRINT'  PhiCjJj*)  =  *;PHI(1,J), 

1340  NEXT  J 

1350  PRINT  :PRINT 

1360  PRINT’Enter  element  number  in  mode  to  change  or*:PRINT 

1370  INPUT*  Enter  <  ZERO  )  to  continue  =>',  RESPONSE 

1380  IF  RESPONSE  =  8  THEN  GOTO  1430 

1398  IF  RESPONSE  <=  8  OR  RESPONSE  )  N  THEN  CLS: GOTO  1310  ELSE  J  =  RESPONSE 

1400  PRINT:  INPUT’Enter  correct  element  value  =>*;PHI(I,J):PRINT 
1410  PRINT'Ph i (* ; J; * )  CHANGED  TO  *;PH1<I,J) 

1420  FOR  K  =  1  TO  2808  :  NEXT  K  :  CLS  :  GOTO  1310 
1438  CLSiNEXT  I 
1440  RETURN 

5180  '  ************************************************************************** 


5005  '  SUBROUTINE  TO  DRIVE  NODE/MASS  INPUT 
5820  LOCATE  5,1,0 

5030  PRINT*  The  progran  has  been  initialized  to  allow  user  input  ct  the'iPRINT 
5840  PRINT*Mode  Shapes  and  Mass  Matrix.  H  you  desire  one  ot  the  other  op-* :PR 
INT 

5050  PRINT'tions,  PRESS  <  escape  key  >;  . .otherwise  =>*:PR1NT:PR1NT 

5055  PRINT*  PRESS  <  space  bar  >  to  continue  )* 

5008  A*  =  INKEY*  :  IF  A*  =  CHM(27)  GOTO  270 
5070  IF  A«  <)  CHR*<32>  GOTO  5000 
5075  CIS  :  LOCATE  5,1,8  i  GOSUB  088 
5880  CLS  :  LOCATE  5,1,8  :  GOSUB  1888 
5090  CLS  :  80SU8  0809 
5100  CLS:  LOCATE  25,1,8 
5110  PRINT 
5128  G0SU8  8008 
5130  CLS:  GOSUB  8288 
5149  CLS  :  GOSUB  18909 
0888 

0010  '  SUBROUTINE  TO  FIND  PARTICIPATION 
0820  '  FACTORS 

6839  UiimwmiHIHIHHmWHiHHiHHWIllHHiHIHHitlHiHIHH 
0049  LOCATE  25,1,8 
0050  T9<1)  =  TIME* 

0008  PRINT*  CALCULATING  PARTICIPATION  FACTORS  =  start  :  *  ;TIMEt 

0089  FOR  1  =  1  TO  M 

0890  FOR  J  =  1  TO  N 

7008  MX(I,J)  =  MASS(J)  *  PHI < I , J) 

7810  NX2<1,J)  =  MX< 1 , J>  *  PHI <X , J) 

7020  NEXT  J 

7830  NEXT  I 
7840  FOR  I  =  1  TO  M 
7858  SUNl(I)  =  MX(I,1) 

7000  SUM2<I)  =  NX2<I,1) 

7870  FOR  J  =  2  TO  N 

7880  SUMKI)  =  SUMKI)  <  MX(I,J) 

7090  SUH2<  I  >  =  SUH2U)  ♦  MX2U,J) 

7100  NEXT  J 

7118  NEXT  I 

7120  FOR  I  =  1  TO  M 

7130  PU)  =  SUMKI)  /  SUM2U) 

7140  NEXT  I 

7150  FOR  I  =  1  TO  H 

7100  MBAR(I)  =  P<I)  *  SUMKI) 

7178  next  I 

7180  LOCATE  25,08,0  :  PRINT*  done:  * jTIME* 

7199  oumr  =  M 

7280  IF  OtWY  )  19  THEN  KX  =  M/2  ELSE  KX  =  M 
7218  L*1:J*« 

7215  LOCATE  5,1,8 

7228  PRINT’*****  PARTICIPATION  FACTORS  WO  MODAL  MASSES  ****** :PRINT:PRINT 


7221  LPRINT1*****  P  FACTORS  AND  MODAL  MASSES  *«***• 

7222  LPRINT 

7225  PRINT*  N  MODAL  PARTICIPATION  MODAL  MSS  1 

7226  PM  =■!*  »#.II«»8»AAAA  Ri .lMMMMt**'‘AA* 

7227  LPR1NT*  N  MODAL  PARTICIPATION  MODAL  MSS  1 

7238  FOR  I  =  L  TO  KX 

7248  PRINT  USING  PM  i  I,P(I),MBAR<1) 

7241  LPRINT  USING  PM  ;  I  ,P<  I )  ,MBAR<  I  > 

7258  J  *  J  ♦  1 
7268  NEXT  I 

7278  IF  Dimr  >  18  AND  J  <  11  THEN  A*  =  INKEYf  iGOTQ  7288  ELSE  GOTO  7311 
7288  PRINT :PRINT :PRINT 
7281  LPRINT :LPR1NT 

7298  PRINT1  PRESS  ANY  KEY  TO  CONTINUE.,.1 

7388  IF  A*  =  11  THEN  GOTO  7388 

7318  KX  =  M  :  L  =  18:  CLS  :  GOTO  7215 

7311  PRINT  -.PRINT  :PR1NT"  PRESS  WY  KEY  TO  CONTINUE...1 

7328  A*  =  INKEY*  :  IF  A*  =  11  THEN  GOTO  7328 
7338  RETURN 

8888  'lilHIHHHHHltHHlHIHHlHitHHiHHiHHmiiHHHlHHlHiHH 

8818  '  SUBROUTINE  TO  CALCULATE  MODAL 
8828  '  WEIGHTS 

8838  limHHIHIHHHHHmiMIHHHIltHtHHmiiHHHWHHHHHHi 
8835  CLSsLOCATE  5,1,8 

8848  PRINT1*****  CALCULATED  MODAL  WEIGHTS  ***«• 

8841  LPRINT1*****  CALCULATED  MODAL  WEIGHTS  **»**■ 

8845  UU*  »  ‘88  888.888888* 

8846  PRINT  :  PRINT  1  N  WEIGHT  IN  KIPS1 

8847  LPRINT  1  N  WEIGHT  IN  KIPS1 
8858  PRINT 

8851  LPRINT 
8878  FOR  1  =  1  TO  M 

8888  UA<I)  =386  «  SUMKI)  *  P<I)  /  1888 
8898  PRINT  USING  WU*  ;  I,WA<I) 

8891  LPRINT  USING  WW*  ;  1,WA<1) 

8895  NEXT  I 

8188  LPR1T4T : LPRINT : PRINT  -.PRINT  jPRINT1  Press  any  Ke 

y  to  continue...1 

8118  A*  =  INKEY* :  IF  A*  =  11  THEN  GOTO  8118 
8128  RETURN 

8288  '************************************************************************* 
8218  '  SUBROUTINE  TO  DETERMINE  GEOMETRY 
8228  '  OF  ANALYSIS  AND  A  <8)  ,  Y  (8)  ,  D  <A) 

8238  'x*********************************************************************** 
8248  PRINT1  Sa»e  information  is  now  needed  to  complete  the  details  about  th 
e1 

8258  PRINT* type  of  vessel  the  structure  1  an  analyzing  is  installed  in  and  whet 

her1 

8268  PRINT’or  not  the  analysis  is  to  consider  elastic  or  elastic-plastic  defor 

na-1 

8278  PRINT1 1 i ons . 1 


8289  PRINT: PRINT 

8298  PRINT’ PLEASE  ENTER  ONE  OF  THE  FOLLOWING  =>’ 

8399  LOCATE  9,28,8:  PRINT’ Ill  SUBMARINE  (  hull  mounted  system  )’ 

8318  LOCATE  11 ,28 ,8 :PRINT’I2]  SUBMARINE  <  decs  mounted  system  )’ 

8328  LOCATE  13 ,28 ,8 sPRIKT’ C33  SUBMARINE  <  shell  plating  mounted  system  >’ 

8338  LOCATE  15,28,8:PRINT’[43  SURFACE  SHIP  (  hull  mounted  system  )’ 

8348  LOCATE  1 7 ,28 ,8 :PRINT’ C 5]  SURFACE  SHIP  <  deck  mounted  system  )’ 

8358  LOCATE  19,28,9:PR1NT’163  SURFACE  SHIP  <  shell  plating  mounted  system  )’ 

8368  A*  =  INKEY* :  IF  A$>=CHM<49)  AND  At  <=  CHRt<54)  GOTO  8398 

8378  IF  A*  =  "  THEN  GOTO  8368 

8388  SEEP:  GOTO  8368 

8398  8  =  UAL (AS) : CLS 

8488  PRINT’ENTER  ==>  [13  Tor  ELASTIC  ANALYSIS’ :  PRINT :  PRINT 

8418  PRINT’  [23  Tor  ELASTIC  -  PLASTIC  AMLYS1S  ’ 

8428  AS  =  INKEYS  :  IF  AS  )=  CHRS<49)  AND  AS  <=  CHRS<58)  GOTO  8458 
8438  IF  AS  =  ’’  THEN  GOTO  8429 
8449  BEEP:  GOTO  8429 
8458  C  =  UAL (AS) 

8468  IF  C=2  AND  <8  =  3  OR  B  =  6)  THEN  60T0  3478  ELSE  GOTO  8588 

8478  PRINT :PRIMT :PRINT :PR1NT’*****  ELASTIC  -  PLASTIC  OPTION  NOT  AVAILABLE  FOR  ** 

***’: PRINT 

8488  PRINT'*****  SHELL  PLATING  MOUNTED  SYSTEMS  *«***’ 

8481  LOCATE  18,28,9 

8482  PRINT’  PRESS  ANY  KEY  TO  CONTINUE  .... 

8483  AS  =  INKEYS  :  IF  AS  =  “  THEN  GOTO  8483 
8498  CLS  .-GOTO  8248 

8588  CLS: LOCATE  5,1,8 

8519  PRINT’ENTER  ===>  [1]  Tor  motion  in  UERT1CAL  direction’ -.PRINT 

8528  PRINT’  [23  Tor  motion  in  ATHUARTSH1P  direction’ :PR1NT 

8538  PRINT’  [33  Tor  motion  in  FORE  AND  AFT  direction’ 

8548  AS  =  INKEYS:  IF  AS  >=  CHRS<49)  AND  AS  <=  CHRS(51)  GOTO  8579 
8558  IF  AS  =  ”  THEN  GOTO  8548 
8568  BEEP:  GOTO  8548 
8578  D  =  UAL(AS) 

8571  IF  8  =  1  THEN  LPRINT’  SUBMARINE  (  hull  mounted  system  )*; 

8572  IF  B  =  2  THEN  LPRINT’  SUBMARINE  (  deck  mounted  system  )’; 

8573  IF  8  =  3  THEN  LPRINT’  SU8NARJNE  (  shell  plating  mounted  system  >’,- 

8574  IF  B  =  4  THEN  LPRINT’  SURFACE  SHIP  (  hull  mounted  system  >’; 

8575  IF  B  =  5  THEN  LPRINT’  SURFACE  SHIP  (  deck  mounted  system  )’; 

8576  IF  B  =  6  THEN  LPRINT’  SURFACE  SHIP  (  shell  plating  mounted  system  )’; 

8577  IF  C  =  1  THEN  LPRINT’  Tor  ELASTIC  ANALYSIS' 

8578  IF  C  =  2  THEN  LPRINT’  Tor  ELASTIC  -  PLASTIC  ANALYSIS’ 

8579  IF  D  =  I  THEN  LPRINT: LPRINT  ’  Tor  motion  in  UERTICAL  direction’ -.LPRINT 

8588  IF  D  =  2  TH0T  LPRINT sLPRINT  ’  Tor  motion  in  ATHUARTSHIP  direct  ion’ :LPRINT 

8581  IF  D  =  3  THEN  LPRINT:LPR1NT  ’  Tor  motion  in  FORE  AND  AFT  direct  ion’ sLPRINT 

8598  FOR  I  =  1  TO  N 

8689  IF  8  0  1  THEN  GOTO  8688 

8618  A8(I)  =  18.4*  <(488+UA<I))/<28+UA(I))> 

8629  U8(l)  =  28!  »  «488+UAU»/(180*UA(I»> 

8638  IF  <0=1  OR  0=2)  AND  C=1  THEN  NULTI=1  :NULT2=1  .-  GOTO  9110  ELSE 

8648  IF  <D=3  AND  C=!)  THEN  MULT1=.4:MULT2=,4  :  GOTO  9118  ELSE 


3658  IF  (0=1  OR  0=2)  MO  C=2  THEN  MULT1=1:NULT2=.5  :  GOTO  9118  ELSE 

8660  MULT1=.4  :  MULT2=.2 

8678  GOTO  9118 

8688  IF  8  <>  2  THEN  GOTO  8788 

8698  A8< I >  =  5.2  *  ((488  ♦  WAU))A28  ♦  UAU))) 

8788  U8U)  =  18!  *  <<480+UA(I))/(180H#(I)>) 

8718  IF  (0=1  AND  OD  THEN  MULT1=1  :NULT2=1  :  GOTO  9118  ELSE 

8720  IF  (0=2  fWD  C=l)  THEN  HULT1=2:MULT2=2  :  GOTO  9118  ELSE 

8738  IF  (0=3  AND  C=l)  THEN  HULT1=.8:MULT2=,8  :  GOTO  9110  ELSE 

8748  IF  (0=1  AND  0=2)  THEN  MULT1=1 :MULT2= .5  :  GOTO  9110  ELSE 

8758  IF  (0=2  ANO  C=2)  THEN  MULT1=2:NULT2=! !  :  GOTO  9118  ELSE 

8768  MULTI  =  .8  :  MULT2  =  .4 

8770  GOTO  9118 

8788  IF  8  <>  3  THEN  GOTO  8858 

3798  A0(1)  =  52  *  ( ( 488+UA<  I )  )A28+UAU ) ) ) 

8808  V8U)  =  188  *  <<488+UA(I))A188+UAU))) 

8810  IF  0=1  THEN  MULT1=1  :MULT2=1  :  GOTO  9118  ELSE 

8820  IF  0=2  THEN  MULT1=.2:MULT2=.2  :  GOTO  9110  ELSE 

8838  MULTI  =  .88  :  MULT2  =  ,88 

8840  GOTO  9118 

8858  IF  B  0  4  THEN  GOTO  8958 

8868  A0(I)  =  20  *  <(37.5  ♦  UA<1>)#<12  +  W(I))A6  ♦  UAU)A2)) 

8878  V0(I)  =  68  *  ((12  ♦  UA(l))/'(6  ♦  UA(I)>) 

8880  IF  <0=1  AND  C=l)  THEN  HULT1=1 :MULT2=1  !  GOTO  9118  ELSE 

8898  IF  (0=2  AND  C=l)  THEN  MULT1=.4:MULT2=.4  :  GOTO  9118  ELSE 

8900  IF  <0=3  AND  C=l)  TH04  MULT1=.2:MULT2=.2  ;  GOTO  9118  ELSE 

8910  IF  (0=1  AND  C=2)  THEN  MULTOl  :MULT2=.5  :  GOTO  9118  ELSE 

8920  IF  (0=2  AND  02)  THEN  MULT1=.4:MULT2=.2  :  GOTO  9110  ELSE 

8938  MULTI  =  .2  :  MULT2  =  .1 

8940  GOTO  9110 

8950  IF  B  0  5  THEN  GOTO  9850 

8960  A0(I)  =  10  *  (<37.5  ♦  UA(I))*<12  ♦  UA(l))/(6  ♦  UA(I)A2)) 

8978  >78(1)  =  38  *  <<12  ♦  UAU))A6  +  l#(l»> 

8980  IF  <0=1  AND  C=l)  THEN  MULT1=1 :MULT2=1  :  GOTO  9118  ELSE 

8990  IF  <D=2  AND  C=l)  THEN  MULT 1= . 4 ;MULT2= . 4  !  GOTO  9118  ELSE 

9000  IF  (0=3  AND  C=!)  THEN  MULT1=.4:MULT2=.4  :  GOTO  9110  ELSE 

9018  IF  (0=1  AND  C=2)  THEN  MULT1=1:MULT2=,5  :  GOTO  9118  ELSE 

9028  IF  (0=2  AND  02)  THEN  MULTi=,4:MULT2=.2  :  GOTO  9118  ELSE 

9838  MULTI  =  .4  :  MULT2  =  .2 

9840  GOTO  9110 

9858  A0(I)  =  40  *  ((37.5  *  UAU))*U2+  UA(I))/(6  ♦  UAU))A2) 

9068  08(1)  =  120  *  ((12  ♦  UAU))A6  ♦  UAU))> 

9888  IF  0=1  THEN  MULT1=1 sMULT2=l  :  GOTO  9118  ELSE 

9898  IF  02  THEN  MULT1=.2:MULT2=,2  :  GOTO  9110  ELSE 

9108  MULT!  =  .1  :  MULT2  =  .1 

9118  00(1)  =  MULT 2  *  V8U)  :  ABU)  =  MULTI  *  ABU)  :  NEXT  1 

9128  FOR  I  =  1  TO  M 

9138  AAU)  *  A0(I)  »  386 

9140  UAU)  =  V8U)  «  W(I) 

9158  IF  ABS(AAU))  >  A8S(VAU))  THEN  DAU)  =  VAU)  ELSE  DAU)  =  mi) 

9168  IF  ABS(DAU))  <  2316  THEN  DAU)  =  2316 


9178  NEXT  1 

9171  CLS:  LOCATE  5,1,8 

9172  FDR  I  =  1  TO  M:  PRINT  ■OACjI;1)  =  ‘;DA<l>f 

9173  NEXT  1 

917*1  PRINTsPRINTsPRINT*  Press  any  key  to  continue  ...* 

9175  A*  =  1NKEY*:IF  A*  =  "  THEN  GOTO  9175  ELSE  CLS  :  LOCATE  5,1,8 

9176  FOR  I  =  1  TO  N  :  LPRINT  :  LPR1NT  *DAC ;I  ;*)  =  ■  ;0A<1) :  NEXT  1  :  LPR1MT 
9180  RETURN 

10988 

10010  '  SUBROUTINE  TO  CALCULATE  NODAL 

18020  '  F  0  R  C  E  S  0  R  0  E  F  L  E  C  T  !  0  N  S 

10930  ' ************************************************************************ 

18048  LOCATE  5,1,8 

18050  PRINT*  Almost  done....  Please  ENTER  the  number  oi  the  type  ol  in Format 
ion* 

10060  PRINT :PRINT*you  desire  as  output  ==>* 

10070  LOCATE  10,1,0  iPRINT*  III  NODAL  FORCES* 

10080  LOCATE  12,1,0  sPWNT*  (23  NODAL  DEFLECTIONS* 

19090  LOCATE  14,1,0  sPRINT*  [31  FORCES  &  DEFLECTIONS* 

10095  LOCATE  16,1,8  sPRINT*  14)  CHANGE  AWLYSIS  TYPE* 

10096  LOCATE  18,1,0  sPRINT*  [5)  EXIT  TO  DOS* 

10189  A$  =  INKEY0!  IF  A*  >=  CHR»<49)  AND  A*  <~  CHRN53)  THEN  GOTO  18148 
10118  IF  A$  *  **  THEN  GOTO  18180 

10130  BEEP  :  GOTO  19180 

10140  FLAG1  =  WAL<A1)  :  CLS 

10150  IF  FLAG1  =  2  THEN  GOTO  10228 

10155  IF  FLAG1  =  4  THEN  GOTO  16160 

18156  IF  FLAG1  =  5  THEN  SYSTEM 

10160  REM  ******  CALCULATE  MODAL  FORCES  ******** 

10170  FOR  I  »  1  TO  N 
19180  FOR  J  =  1  TO  N 

10190  F(1,J)  =  PCI)  *  NASS(J)  *  PHI(I,J)  *  DA(I) 

10280  NEXT  J 

18210  NEXT  I 

18220  IF  FLAG1  =  1  THEN  GOTO  12890 

18238  REM  ******  CALCULATE  MODAL  DEFLECTIONS  ********* 

10240  FOR  I  =  1  TO  M 
10258  FOR  J  =  1  TO  N 

10260  DELTA(l.J)  =  PHI(1,J)  *  P<I>  *  DA<I)  /  U(I)*2 

10290  NEXT  J 

18309  NEXT  I 

10310  FOR  J  =  1  TO  N 

19329  BIGGEST  =  8 

10330  FOR  1  *  1  TO  N 

10348  IF  BIGGEST  <  ABS(DELTA<J,J))  THEN  BIGGEST  =  ABS<DELTA(I,J))s  K=I 

10350  NEXT  I 

10360  S  -  8 

10370  FOR  1  *  1  TO  N 

10388  IF  1  *  K  THEN  GOTO  10480 

10390  S  =  S  *  DELTA(I,J)*2 

10489  NEXT  I 


10418  DEFL(J)  =  ABS(D£LTA<K,J)>  ♦  SQR<Sj 
13420  NEXT  J 

12000  REN  *********  OUTPUT  RESULTS  *********** 

12010  CLS  :  LOCATE  5,1,3 

12828  PRINT  *  ENTER  Ill  TO  ECHO  RESULTS  TO  SCREEN' :P 

R1NT:  PRINT 

12830  PRINT  *  [21  TO  ECHO  RESULTS  TO  PRINTER*: 

PR1NT:PR1NT 

12040  At  =  INKEYt 

12058  IF  A*  >=  CHRt<49)  AND  At  <=  CHRK50)  THEN  SOTO  12078 
12048  IF  At  =  "  THEN  GOTO  12048  ELSE  8EEP:  GOTO  12848 
12070  FLAG2  =  UAL<At) 

12080  IF  FLAG2  =  2  THEN  GOTO  12380 
12090  IF  FLAG1  =  2  THEN  GOTO  12198 
12100  FOR  1  =  1  TO  M 

12118  CLS  :  LOCATE  5,1,8  :  PRINT*  ***  M  0  D  A  L  FORCES  FOR  MODE 
<•;!}•)  [lbs]  *»*•  : PRINT 

12111  LPRINT:LPRINT*  ***  M  0  0  A  L  FORCES  FOR  MODE  (  *  ;I  ;*  )  [lb 
si  *«*  ;  LPR1NT 

12112  LPRINT  *  N  *:  LPRINT 

12114  PRINT  *  N  *  :  PRINT 

12115  FFt  =  *0#  ii.HMHHtr-'*''* 

12128  FOR  J  =  1  TO  N 

12138  PRINT  USING  FFt  ;  J,F(J,J> 

12131  LPRINT  USING  FFt  ;  J,FU,J) 

12140  NEXT  J 

12150  PRINT :PRIKT  :PRIKT :PR1NT 

12140  PRINT*  Press  any  key  to  continue  ...* 

12170  At  =  INKEYt:  IF  A*  »  **  THEN  GOTO  12178 

12188  NEXT  I 

12190  IF  FLAG1  =  1  THEN  GOTO  12280 
122? 8  CLS:  LOCATE  5,1,0 

12210  PRINT*  *«  DEFLECTIONS  AT  MASS  ATTACHMENT  PC 
I  N  T  S  ***“ :PRINT:PRINT 

12211  LPRINT  :  LPRINT1**  DEFLECTIONS  AT  MASS  ATTACH  ME 
NT  POINTS  ***:LPRINT 

12215  DDt  =  *01.088888  * 

12220  FOR  I  =  1  TO  N 

12230  PRINT* DELTA  <*;!;*)  =  *; 

12231  LPRINT*  DELTA  <* ;I ;*)  =  *; 

12232  LPRINT  USING  ODt  ;  DEFL< I )  ; 

12233  LPRINT* in.* 

12235  PRINT  USING  DDt  ;  DEFL(I)  ; 

12234  PRINT  *in.* 

12240  NEXT  I 

12250  PRINT  :PRINT:PRINT 

12251  LPRINT  sLPRINT 

12240  PRINT*  Press  any  key  to  continue  ...* 

12270  At  =  INKEYt  :  IF  At  =  **  THEN  GOTO  12270 
12288  CLS  :  LOCATE  5,1,0 
12298  GOTO  10008 
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12388  PRINT  •  OPTION  NOT  AVAILABLE  YET;—)  USE  <  Prtsc  >  KEY  * 

12318  FOR  I  »  1  TO  3888  :  NEXT  1  GOTO  12888 

15808  ' ************************************************************************ 
15818  '  SUBROUTINE  TO  INPUT  MASS/MOOES 
15828  '  FROM  DISK  FILE 

15838  'IHlHHHOiHHHHiHiHHIiHlttHHHHHiHHHHtHHHHiHHiH 

15848  CLS  :  LOCATE  5,1,8 

15858  PRINT*  The  program  has  bun  initialized  to  allow  user  input  of  the*. -PRINT 
15848  PRlNT’Mode  Shapes, Frequencies,  and  Masses  via  an  input  tile.':PRINT 
15878  PRINT  :  PRINT*!!  you  desire  one  ot  the  other  options,* :PRINT:PRINT 

15888  PRINT*  PRESS  (escape  key);  . otherwise  =>*:PRINT  :  PRINT 

15898  PRINT*  PRESS  (space  bar)  to  continue  * 

15188  A*  =  INKEY*  :  IF  A*  =  CHR*<27)  GOTO  278 
15118  IF  A*  ()  CHR*<32)  GOTO  15188 
15121  CLS:LOCATE  5,1,8 

15138  INPUT  'Enter  the  dimension  o!  the  Mass  Matrix  =>  *;N 
15148  IF  N  (=  58  THEN  GOTO  15178 

15158  PRINT : PRINT : PRINT* *****  PROBLEM  S12E  IS  LIMITED  TO  58  DEG  OF  FREEDOM  **** 
** 

15148  FOR  1  =  1  TO  1888  :  NEXT  1  :  GOTO  15138 

15178  PRINT:PRlNT:lNPUT*Enter  the  number  o!  modes  you  wish  to  use  =>  *;M 
15175  PRINT :PRlNT:INPUT*Enter  the  input  tile  name  =>  *,FIN* 

15188  CLS 

15198  'HIHIlHHHHHiHHilimiHHIHKHHIHKHHHiHliltHHiHiKH 

15288  '  INPUT  MASS, FREQ,  ANO  PHI 

15218 

15228  OPEN  FIN*  FOR  INPUT  AS  R1 
15248  FOR  I  «  1  TO  N 
15258  INPUT  81,  MASS(I) 

15248  NEXT  I 

15278  FOR  I  =  1  TO  M 

15288  INPUT  II,  U< I) 

15298  U(I)  =  U(l)  *  4.283185 

15388  NEXT  I 

15318  FOR  J  =  1  TO  N 

15338  FOR  I  =  1  TO  M 

15348  INPUT  #1, PHI (I, J) 

15358  NEXT  I 
15348  NEXT  J 

15378  'HHiHtHHHHHHHiHHHHHHHIlHHHIlHHHiHIHHHHHilHtH 
15388  '  OUTPUT  MSS,  OMEGA,  WD  PHI 

15398  'HHHHKIHIHHHHHiHiHHIiflHilliHHiHIIHIHillHHHHHHIlKI 
15488  LOCATE  5,1,8  :  PRJNT'MSS  MATRIX  DIAGONAL  ELEMENTS  *  :  PRINT 

15481  LPRINT1*****  MASS  D1A60NAL  ELEMENTS  ****** 

15482  LPRINT 

15485  Itl*  =  *88  88.88888'*“ 

15418  FOR  1  =  1  TO  N 
15428  PRINT  USING  tttt  ;  I,MSS(I) 

15421  LPRINT  USING  HI*  ;  I,MASS(I) 

15438  NEXT  I 
15431  PRINT  :PR1NT* 


....PRESS  ANY  KEY  TO  CONTINUE1 


15432  A*  =  INKEYI  :  IF  A$  =  “  THB4  GOTO  15432 
15434  CLS 

15448  PRINT :PRINT :PRINT 

15458  PRINT’ FREQUENCIES  USED  IN  ANALYSIS  (RPS)‘:PR1NT 

15451  LPR1NT 

15452  LPRINT'*****  FREQUENCIES  USED  !N  ANALYSIS  (RPS) 
****** : LPRINT 

15468  FOR  I  =  1  TO  M 
15478  PRINT  I;*  *;U(I) 

15471  LPRJNT  I;’  *{W<1) 

15488  NEXT  I 

15481  PRINT SPRINT*  ....PRESS  AMY  KEY  TO  CONTINUE* 

15482  A*  =  INKEY*  :  IF  A*  =  ’*  THEN  GOTO  15482 
15484  CLS 

15490  PRINT-.PRINT sPRINT 

15491  LPRINT: LPRJNT 

15508  PRINT’MODES  USED  IN  ANALYSIS*  sPRINT 

15501  LPRINT'*****  NODES  USED  IN  ANALYSIS  ****** :  LPRINT 
15510  K  =  M  :  U$  =  *U.»i|»ll**“  * 

15520  IF  K  >  5  THEN  GOTO  15688 
15530  FOR  J  =  1  TO  N 
15540  FOR  1  =  1  TO  N 

15550  PRINT  USING  U»;  PHK I, J); 

15551  LPRINT  USING  Ut;  PHKI.J); 

15560  NEXT  I 

15578  PRINT  *  * 

15571  LPRINT  *  ’ 

15588  NEXT  J 

15590  PRINT  sPRINT 

15591  LPRINT  :  LPRINT 

15592  PRINT:  PRINT*  ....PRESS  ANY  KEY  TO  CONTINUE* 

15593  A$  =  INKEY*  :  IF  A*  =  **  THEN  GOTO  15593 

15594  CLS:  GOTO  15740 

15600  FOR  J  =  1  TO  N 

15610  FOR  I  =  1  TO  5 

15620  PRINT  USING  U$;  PHKI.J); 

15621  LPRINT  USING  U$;  PHKI.J); 

15638  NEXT  I 

15640  PRINT  *  * 

15641  LPRINT  *  * 

15650  NEXT  J 

15660  PRINT :PRINT:LPRJNT  sLPRINT 

15661  PRINT*  ....PRESS  ANY  KEY  TO  CONTINUE* 

15662  A»  =  INKEY*  :  IF  A*  =  **  THEN  GOTO  15662 

15663  CLS:  PRINT :  PRINT 

15678  FOR  J  =  1  TO  N 

15680  FOR  I  =  6  TO  N 

15690  PRINT  USING  U»;  PHKI.J); 

15691  LPRINT  USING  U*;  PHKI.J); 

15708  NEXT  I 

15710  PRINT  *  * 
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15711  LPRINT  '  * 

15728  NEXT  J 

15738  PRINT  :PRINT :LPRIMT : LPRINT 

15731  PRINT*  . . . .PRESS  ANY  KEY  TO  CONTINUE* 

15732  A*  =  INKEYS  :  IF  A*  =  ’’  THEN  GOTO  15732 

15734  CIS:  PRINT  :  PRINT 

15735  CLOSE  #1 

1 5748  ' ************************************************************************ 
15758  '  NORMALIZE  MODE  SHAPES 

15748  '************************************************************************ 
15778  FOR  1  =  1  TO  M 
15788  BIG  =  8! 

15798  FOR  J  =  1  TO  N 

15888  IF  BIG  <  PHIU.J)  THEN  BIG  =  PH1U.J) 

15818  NEXT  J 

15828  FOR  J  =  1  TO  N 

15838  PHI ( I , J)  =  PHI<I,J)/BIG 

15848  NEXT  J 

15858  NEXT  1 

1 5848  ' ************************************************************************ 
15878  '  OUTPUT  NORMALIZED  MODE  SHAPES 

15888  /  ************************************************************************ 
15898  PRINT ’NORMAL  I  ZED  MODES  USED  IN  ANALYSIS’  '.PRINT  :  LPRINT 

15891  LPRINT’*****  NORMALIZED  MODES  FOR  ANALYSIS  ****** 

15892  LPRINT 
15988  K  =  M 

15918  IF  K  >  5  THEN  GOTO  15988 

15921  FOR  J  »  1  TO  N 

15938  FOR  I  =  1  TO  M 

15948  PRINT  USING  Ift;  PHKI.J); 

15941  LPRINT  USING  US;PHI<I,J); 

15958  NEXT  I 

15941  PRINT  *  ’ 

15941  LPRINT  ’  ’ 

15978  NEXT  J 

15971  PRINT:  PRINT’  ....PRESS  WY  KEY  TO  CONTINUE’ 

15972  A$  =  INKEYS  :  IF  AS  =  ”  THEN  GOTO  15972 

15973  CLS 

15988  LPRINT  :  LPRINT  :  GOTO  14128 

15998  FOR  J  =  1  TO  N 

14888  FOR  I  =  1  TO  5 

14818  PRINT  USING  OS;  PH1(I,J); 

14811  LPRINT  USING  US;  PHIU.J); 

14828  NEXT  1 

14838  PRINT  1  ’ 

14831  LPRINT  ’  ’ 

14848  NEXT  J 

14158  PRINT:PRJNT:LPRJMT:LPR1NT 

14851  PRINT’  ....PRESS  WY  KEY  TO  CONTINUE’ 

14852  A$  =  INKEYS  :  IF  AS  =  ’’  THEN  GOTO  15852 

14853  CLS 


16069  FOR  J  =  1  TO  N 

16078  FOR  1  s  6  TO  N 

16888  PRINT  USING  Wj  PHI< I, J); 

16081  LPR1NT  USING  U$;  PHH1.J) ; 

16098  NEXT  I 

16188  PRINT  *  1 

16101  LPR1NT  1  * 

16110  NEXT  J 

16120  PRINT  ;PR!NT:LPR1NT :LPR1NT 

16121  PRINT*  ....PRESS  (M  KEY  TO  CONTINUE* 

16122  A$  =  INKEYi  :  IF  A*  =  **  THEN  GOTO  16122 

16138  CLS  :  G0SU8  6808 

16140  CLS  :  LOCATE  21,1,0 

16158  PRINT 

16155  CLS  :  GOSU8  8888 

16160  CLS  :  GOSU8  8280 

16178  CLS  :  GQSUB  18888 

28808  CLS 
20818  LOCATE  5,1,0 

20828  PRINT****  YOU  MUST  FIRST  USE  THE  PROGRAM  JAC08I  TO  SOLVE  EIGEWALUE  *** 
20038  PRINT****  PROBLEM.  IT  WILL  CALCULATE  NATURAL  FREQUENCIES  MD  MODE  *** 
20848  PRINT'***  SMPES  THAT  CAN  THEN  BE  USED  TO  RUN  THE  DDAM  PROGRAM.  *** 
28858  PRINT:  PRINT:  PRINT: 

20868  PRINT*  .  ENTER  <F3)JAC0BI (RETURN)* :PRINT 

20078  PRINT*  <F2>* 

28888  STOP 

25880  CLS'.PRINT'***  USER  MAY  USE  OPTION  4  TO  INSTALL  HIS/HER  SUBROUTINE  ***• 
25010  PRINT*  ***  BEGINNING  AT  LINE  25008.**** 

25028  FOR  1  *  1  TO  5808  :  NEXT  I  :  GOTO  278 


PROGRAM  JACOBI 


INSTRUCTIONS  FOR  USE 


C . 

The  -following  BASIC  program  may  be  used  as  is  or 
modified  by  the  user  to  solve  the  Eigenvalue  problem  once 
the  equations  of  motion  are  formulated.  It  is  self  explan 
atory  in  execution  with  the  exceptions  that  the  A  matrix 
replaces  the  stiffness  matrix  and  the  B  matrix  replaces  t 
mass  matrix.  The  program  only  requires  that  the  upper  tr 
angle  elements  of  the  stiffness  matrix  be  suppl ied  along 
with  the  diagonal  elements  of  the  mass  matrix. 


0.  PROGRAM  JACOBI  LISTING  Adapted  from  ref.  13) 


4  CIS 

5  REM  ®CYCL1C  JACOBI  METHOD 

18  DIM  A<18,18) ,  8(10) ,U< 10 ,10) 

19  REM  = INPUT  DATAf= 

28  GOSUB  78 

28  REINITIALIZE® 

38  GOSUB  668 
48  GOSUB  988 

45  REM  =0UTPUT  RESULTS® 

58  GOSUB  758 
48  GOTO  3888 
65  REM  =  INPUT  DATAf® 

78  PRINT 

88  PRINT  *  CYCLIC  JACOBI  METHOD1 
168  PRINT 
188  PRINT 

198  PRINT  1  ENTER  THE  SIZE  OF  PROBLEM  (i.e.  no.  of  rows)1 
288  PRINT 


218  INPUT  N  -  ;  N 

228  LPRINT  'PROBLEM  DIMENSION  =  *i  N 

238  INPUT  'NUMBER  OF  SIGNIFICANT  FIGURES  =  *  j  SI 

248  LPRINT 

258  PRINT 

268  PRINT  '  ENTER  THE  UPPER  TRIANGULAR  ELEMENTS  OF  tttTRIX  A  (  by  columns  ): 
278  PRINT 

288  LPRINT  'UPPER  ELEMENTS  OF  MATRIX  A  =' 

298  LPRINT 

320  FOR  J  =  1  TO  N 

338  PRINT 

340  PRINT  '  ENTER  UPPER  PART  OF  COLUW  '  j  J 

350  PRINT 

368  FOR  I  =  1  TO  J 

370  PRINT  '  Al’jljVjJ;')  5  'i 

388  INPUT  A<1 , J) 

381  LPRINT  '  A<';IiY;Jj')  »  AU,J> 

382  LPRINT 
390  NEXT  1 
468  PRINT 
448  PRINT 
450  NEXT  J 

460  FOR  I  =  1  TO  N 
470  II  =  I  -  1 
488  FOR  J  =  1  TO  II 
498  AO,J)  =  A(J,I) 

580  NEXT  J 
510  NEXT  I 
520  PRINT 

530  PRINT  '  ENTER  ELEMENTS  OF  DIAGONAL  MATRIX  B:' 

531  LPRINT  '  ELEMENTS  OF  DIAGONAL  MATRIX  8  * 

532  LPRINT 
548  PRINT 

570  FOR  I  =  1  TO  N 
588  PRINT  '  B<* jl ;•)  * 

590  INPUT  BO) 

591  LPRINT  '  B('ili')  =  "^XIX 
680  NEXT  I 

618  PRINT 
611  LPRINT 
650  RETURN 

655  REM  =1NIT1AL1ZE  TOLERANCE  AND  MAX  NO.  OF  ROTATIONS  = 

668  Z  *  2  *  SI 
678  T1  «  1/08AZ) 

680  PRINT 

690  CLS:  PRINT  'TOLERANCE  =  ';T1 

691  LPRINT  ’TOLERANCE  =  ' ;T1 


788  R  »  5*N*N 
718  R1  =  8 
728  T2  =  .1 
738  N1  =  N  -  1 


748  RETURN 

744  REN  =OUTPUT  E1GENSOLUTION= 

750  PRINT 
748  PRINT 

778  PRINT  ■SOLUTION* 

771  IPRINT 
788  PRINT 

781  IPRINT 

782  IPRINT  ■******»****»**«*  SOLUTION  ***»*********«*«• 

783  IPRINT 

800  PRINT  ■  NUMBER  OF  ROTATIONS  REQUIRED  =  *;R1 
881  LPRINT  ’NUMBER  OF  ROTATIONS  REQUIRED  =  *;R1 

802  LPRINT 

803  LPRINT 
818  PRINT 
820  PRINT 

838  FOR  J  =  1  TO  N 
848  PRINT 

850  PRINT  ’  EIGENVALUE  ’jJj’  IS  ’ ;B<J) 

851  LPRINT  ’EIGENVALUE  ’jJj’  IS  ’ ;B<J) 

852  LPRINT 

853  LPRINT  MTS  EIGENVECTOR  IS  ’ 

848  PRINT 

870  PRINT  ’  ITS  EIGENVECTOR  IS  ’ 

880  PRINT 

890  FOR  I  *  1  TO  N 
980  PRINT  U(1,J) 

901  LPRINT  U(1,J) 

982  LPRINT 
9 10  NEXT  1 

920  PRINT 

921  LPRINT 

930  PRINT  ’  PRESS  SPACE  BAR  TO  CONTINUE  ’ 

948  IF  INKEYS  =  ’  ’  THEN  942 

941  GOTO  940 

942  CLS 
958  PRINT 
940  NEXT  J 
978  RETURN 

975  REM  =»E 1 604PROBLB1  SOLlfTION= 

988  60SUB  1138 

985  REM  ^PERFORM  ONE  CYCLE  OF  ROTATIONS 
991  GOSUB  1290 

995  RBI  =  CHECK  TOLERANCE  = 

1088  IF  XI  <  T1  THEN  1110 

1005  RBI  =  CHECK  NO.  OF  ROTATIONS  = 

1010  IF  R1  >  R  THEN  1840 
1820  T?  =  .1  ♦  XI 
1030  GOTO  990 
1040  PRINT 

1050  PRINT  ’  **»  ERROR  **«’ 
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1851  LPRINT  ‘  it*******************  ERROR  ******************** 
1868  PRINT 

1878  PRINT  *  NO  CONVERGENCE  ATTAINED* 

1871  LPR1NT  *  NO  CONVERGENCE  ATTAINED  * 

1888  PRINT 
1881  LPR1NT 

1898  PRINT  *  UITH  ,jRl;*  ROTATIONS* 

1891  LPRINT  *  UITH  *;R1;*  ROTATIONS* 

1188  END 
1118  G0SU8  1838 
1128  RETURN 
1130  FOR  3  =  1  TO  N 
1148  FOR  J  =  1  TO  N 
1158  U(I,J)  =  8 
1168  NEXT  J 
1170  U<I,1)  =  1 
1188  NEXT  I 
1198  FOR  I  =  1  TO  N 
1208  81  =  SQR<8( I ) ) 

1210  B<I)  =  1/81 

1221  NEXT  I 

1230  FOR  1  =  1  TO  N 

1240  FOR  J  =  1  TO  N 

1251  A<I,J)  =  8< I )  *  A< I , J)  *  8(J) 

1260  NEXT  J 
1270  NEXT  I 
1280  RETURN 
1290  XI  =  8 
1300  FOR  K  =  1  TO  N1 
1310  K1  =  K  ♦  I 
1320  FOR  L  =  K1  TO  N 
1330  A1  =  A(K,K) 

1340  A2  =  A(K,L) 

1358  A3  =  A(L,L) 

1360  X  =  A2  *  A2  /<A1  *  A3) 

1365  REN  =  CHECK  IF  ROTATION  IS  NEEDED  = 

1370  IF  X  >  XI  THEN  1398 
1388  GOTO  1800 
1390  XI  =  X 

1408  IF  X  <  T2  THEN  1888 

1410  R1  =  R1  ♦  1 

1415  RBI  =  COMPUTE  ANGLE  = 

1420  IF  A1  =  A3  THEN  1470 
1438  Z  =  .5  MAI  -  A3)  /A2 
1440  21  =  1  ♦  1/  <Z  »  Z) 

1450  T  =  -Z  *  <1  ♦  SOR(ZD) 


1460  GOTO  1488 
1470  T  =  1 

1480  C  »  1  /  SQR(  1  ♦  T  *  T) 
1490  S  =  C  *  T 
1588  S2  =  S  *  S 


1518  C2  2  C  *  C 
1528  AIK,L)  2  8 

1525  REM  =  TRANSFORM  DIAGONAL  ELEMENTS  = 

1538  A8  2  2  *  A2  *  C  *  S 

1548  A(K,K)  =  A1  *C2  ♦  A8  ♦  A3  *  S2 

1558  A<L,L)  2  A1  *  S2  -  A8  ♦  A3  *  02 

1555  REM  22  TRANSFORM  OFF  DIAGONAL  ELEMENTS  22 

1568  FOR  I  =  1  TO  N 

1578  IF  !  <  K  THEN  1688 

1588  IF  I  >  K  THEN  1641 

1598  GOTO  1748 

1688  A8  2  A< I ,K) 

1618  ACI,K)  2  C  *  A8  ♦  S  *  AU,L> 

1628  A(1,L)  2  -S  *  AB  ♦  C  *  A(1,L) 

1638  GOTO  1748 
1648  IF  I  <  L  THEN  1678 
1658  IF  I  )  L  THEN  1718 
1668  GOTO  1740 
1678  A0  2  A(K,1) 

1688  AIK, I)  2  C  *  A8  ♦  S  *  A1I,L) 

1698  A(I,L)  2  -S  *  A8  ♦  C  *  A(I,L) 

1788  GOTO  1748 
1718  A0  2  AIK, I) 

1728  AIK, I)  2  C  *  A8  ♦  S  «  AIL, I) 

1738  AIL, I)  2  -S  *  A8  ♦  C  *  AIL, I) 

1748  NEXT  I 

1745  REM  22  TRANSFORM  fttTRIX  U  TO  GENERATE  EIGENVECTORS 
1758  FOR  I  2  1  TO  N 
1768  U0  2  11(1  ,K) 

1770  UII ,K)  2  C  *  U8  ♦  S  *  Oil ,L) 

1788  UII ,L)  2  -S  *  U8  ♦  C  *  UII  ,L) 

1790  NEXT  I 
1388  NEXT  L 
1810  NEXT  K 
1828  RETURN 

1825  REM  22  NORMALIZE  EIGENVECTORS  22 

1830  FOR  I  2  1  TO  N 

1848  FOR  J  2  1  TO  N 

1858  U(1,J)  2  UII.J)  *  BID 

1868  NEXT  J 

1870  NEXT  I 

1880  FOR  I  2  1  TO  N 

1890  BID  2  All, I) 

1900  NEXT  I 

1905  RBI  «  ORDER  E1GENSOLUTJON  = 

1910  FOR  I  2  1  TO  N! 

1920  11  2  I  ♦  1 
1938  Z  2  811) 

1940  N  2  J 

1958  FOR  J  2  11  TO  N 


1968  IF  Z  1  8<J)  THEN  1998 


*>  •• 


APPENDIX  8:  PROGRAMS  USED  TO  CONCERT  VELOCITY  HISTORIES 
TO  DEFLECTION  HISTORIES 


20  '  PROGRtfl  **  SIMPS1  »*  USED  FOR  CASE  1 

30  '  INTEGRATES  VELOCITY  HISTORY  URT  TINE  TO  OBTAIN  DEFLECTIONS 

41  - - 

50  CLS  :  01N  A$<3) 

60  F<1)=8  :  F<2)=0  :  F<3)=0  :  F0RM1$=*  8.1*18  *  :  F0RN2$=‘  «*.««««««  * 

70  INPUT'  01TER  THE  WME  OF  (ti  INPUT  FILE...,,,,... .  \FILE1N*  :  PRINT 

88  INPUT*  ENTER  THE  ORDIMTES  TO  INTEGRATE  OVER .  *  .NUMBER  :  PRINT 

90  INPUT’  ENTER  THE  SIZE  OF  THE  INTEGRATION  INTERVAL . ',H 

180  CLS 

118  OPEN  F1LEJN*  FOR  INPUT  AS  II  :  OPEN  ’D1SPR2.0UT’  FOR  OUTPUT  AS  *2 

121  ' - 

130  '  USE  SIMPSON'S  1/3  RULE  TO  CALCULATE  DEFLECTIONS 

140  ' - 

150  INPUT  Ki ,N,TINE(8),VEL1K8),V£L12{0),VEL13(8) ,VEL6K8) ,VEL62(8),VEL63(8) 

160  FOR  I  =  1  TO  2 

178  INPUT  II ,N ,TIHE< I ) ,VEL1 KI) ,VEL1 2< I ) ,VEL13< I ) ,VEL61 < 1 ) ,VEL62( 1 ) ,VEL63< I ) 
180  NEXT  I 

190  VIA  =  VEL6K0)  -  VEL1K0) 

280  V1B  =  VEL6K1)  -  VELIKI) 

210  VIC  *  VEL6K2)  -  VEL1K2) 

220  Q1  *  (VIA  +  4*V1B  ♦  VIC)  *  H  /  3 
230  V2A  =  VEL62(0)  -  VEL12(8) 

240  V2B  *  VEL62U)  -  VEL12(1) 

258  V2C  =  VEL62(2)  -  VEL12(2) 

260  82  =  (V2A  ♦  4*V2B  ♦  V2C)  *  H  /  3 
270  V3A  =  VEL63(0)  -  VEL13(8) 

280  V3B  =  VEL63U)  -  VEL13U) 

290  V3C  =  VEL63(2)  -  VEL13<2) 

300  83  =  (V3A  ♦  4*V38  ♦  V3C)  *  H  /  3 

310  F(l)  =  F(l)  ♦  81 

320  F<2)  =  F<2)  ♦  82 

330  F(3)  =  F(3)  ♦  83 

340  T  =  TIME(l)  »  1880 

350  FOR  I  *  1  TO  3 

360  IF  A8S«MW(D)  >  ABS(F(D)  THEN  GOTO  388 
370  OttXU)  =  F(I)  :  TNW<I)  =  T 

380  NEXT  I 

398  PRINT  82,  '  T  =  PRINT  12,  USING  FORMII-.T; 

480  PRINT  12,  *(151  c  OX  =  ';  :  PRINT  12,  USING  FQRN2$;F<1); 

418  PRINT  12,  *  in.  DY  =  * ; :  PRINT  12,  USING  F0RM2»;F(2); 

428  PRINT  12,  'in.  DZ  =  PRINT  12,  USING  F0RM2*;F<3) ; :PRINT  82,'in.' 

430  VEL1K8)  =  VELIK2)  :  VEL12<0)  =  VEL12(2)  :  VELI3T8)  =  VEL13<2) 


449  VEL6K8)  =  VEL6K2)  :  VEL62<0)  =  VEL62<2>  :  VEL63<8)  =  VEL63(2> 

458  IF  H  ♦  2  <  NUMBER  THEN  GOTO  168 

468  f»<l)  =  *(X)  =  1  :  A4<2)  =  *<Y)  =  *  ;  A*(3>  =  '(2)  =  1 

478  FOR  I  =  1  TO  3 

488  PRINT  82,  *TWX*  ;A$<  1 ) ;  :  PRINT  82,  USING  FORMIC  ;TmX<  1 )  { 

485  PRINT  82,’  nsec  0N4X*;A«I)j:  PRINT  82,  USING  FQ812*,*QMAX(1) 

486  PRINT  82, ‘  in.1 
498  NEXT  I 

588  CLOSE  81  i  CLOSE  *2 


PROGRAM  **  S1MPS2  **  USED  FOR  CASES  II  AND  HI 
INTEGRATES  VELOCITY  HISTORY  URT  TINE  TO  OBTAIN  DEFLECTIONS 


58  CLS  :  DIM  A4<3> 

68  F< 1 )=«  :  F<2)=8  :  F0RM14=*  89.8119  *  :  FORM24=*  88.899881  * 

70  INPUT*  ENTER  THE  NWE  OF  AN  INPUT  FILE .  *  .FILEIN*  :  PRINT 

88  INPUT*  ENTER  THE  ORDINATES  TO  INTEGRATE  OVER .  *, NUMBER  :  PRINT 

90  INPUT*  ENTER  THE  SIZE  OF  THE  INTEGRATION  INTERVAL .  *,H 

190  CLS 

118  OPEN  F1LEJN*  FOR  INPUT  AS  81  :  OPEN  ‘D1SPR2.0UT*  FOR  OUTPUT  AS  82 
115  PRINT  82, *H  *  *;H;*  SEC* 


USE  SIMPSON'S  1/3  RULE  TO  CALCULATE  DEFLECTIONS 


150  INPUT  II ,M,T1ME(0> ,VEL12<0> ,VEL92<0) ,VEL52<0> ,VEL112<0) ,VEL192<0> .VEL152L 0> 

168  FOR  I  =  1  TO  2 

178  INPUT  81  ,M,T1ME< 1 ) ,VEL1 2< I ) ,VEL92( I ) ,VEL52< 1 ) ,VEL1 12(1) ,VEL1 92(1)  ,VEL 152(1) 
188  NEXT  I 

198  VIA  =  <VEL12(8)  ♦  VEL92<8)>/2  -  VEL52I8) 

288  VI B  =  (VEL12U)  ♦  VEL92U))/2  -  VEL52U) 

210  VIC  =  (VEL12I2)  ♦  VEL92<2))/2  -  VEL52<2) 

220  Q1  =  (VIA  ♦  4*V1B  ♦  VIC)  *  H  /  3 

230  V2A  =  (VEL112<8)  ♦  VEL192(8))/2  -  VEL152(0) 

240  V2B  =  (VEL112U)  ♦  VEL192<l))/2  -  VEL152U) 

258  V2C  =  (VEL112<2)  ♦  VEL192(2))/2  -  VEL152I2) 

260  Q2  =  (V2A  ♦  4*V2B  ♦  V2C)  *  H  /  3 

310  F(l)  =  F<1)  ♦  Q1 

328  F(2)  =  F(2)  ♦  02 

348  T  =  TIME(l)  *  1998 

358  FOR  1  =  1  TO  2 

368  IF  A8S(0MAX(D)  >  A8S(F(D)  THEN  GOTO  380 
370  OMAX(I)  =  F(I)  i  TMAX(I)  *  T 

388  NEXT  1 

398  PRINT  82,  *  T  *  *j.-  PRINT  82,  USING  FORMAT; 

409  PRINT  82,  *ns*c  D(5,Y)  =  *:  :  PRINT  12,  USING  F0RM24;F<1); 


418  PRINT  *2,  ‘in.  0<15,Y)  =  *;:  PRINT  82,  USING  F0RM2*iF<2) ; 

428  PRINT  *2,  ‘in." 

438  VEL1 2(8 )  *  VEL12<2)  :  VEL92<8)  =  UEL92<2)  :  VEL52<8)  =  UEL52<2) 

448  VEL112<8)  =  gEL112<2)  :  gEL192<0)  =  VEL192(2)  :  gEL152<8)  =  gEL152(2) 

458  IF  N  +  3  <  NUMBER  THEN  GOTO  148 

448  A*<1)  =  ■  < 5 |Y)  =  *  :  Al<2)  =  "<15, Y)  =  * 

478  FOR  !=  1  TO  2 

488  PRINT  82,  ’TNWjAid);  :  PRINT  82,  USING  FORNUjTNAXU) ; 

485  PRINT  82,'  nsec  UMAX’ ;A$< I ) ; :  PRINT  82,  USING  F0RM2$;QMAX(1); 

484  PRINT  82,'  in.' 

490  NEXT  1 

588  CLOSE  81  :  CLOSE  82 


APPENDIX  C  -  TYPICAL  CASE  I  ANALYSIS 

As  noted  in  Chapter  V,  the  -first  step  in  each  anal  ys 
was  the  construction  of  an  appropriate  SAPIV  model  of  th 
equipment  substructure  under  consideration.  In  this  cas 
the  model  represents  a  1,908  lb  value  attached  to  the  fr 
end  of  a  cantilever  foundation. 

A.  SAP IV  INPUT  DATA  FOR  CASE  I 

The  following  data  file  was  used  to  generate  the  mod 
information  for  the  Case  I  model  [Ref.  10]. 

HEADING  CARD 

CASE  1  —  VALVE  ON  CANTILEVER  FOUNDATION 
MASTER  CONTROL  CARD 

7.1 .8.6.1  ,9,9, 0,9,0 
NODAL  PT  DATA 

0,1  ,1  ,0,8,1  ,1  ,1  ,9.  ,8.  ,9.  ,1  ,9.0 
0,2,1 ,8,0, 0,1 ,1 ,9. ,8. , 5 • ,1 ,0.9 
0,6,1  ,  0 , 0 , 0 , 1  ,1  ,  0  .  ,  0  .  ,  25  .  ,1  ,0.0 
0,7,1 ,1 ,1 ,1 ,1 ,1 ,10. ,0. ,10. ,0.0 

TYPE  2  -  3D  BEAM  ELEMENTS 

2.5.1  ,0,1 

MATERIAL  PROPERTIES 

1 , 30 . 0E6 , 0 . 3 , 7 . 29E-4 ,0.0 
ELEMENT  PROPERTY  CARD 

1,17.64,0.0,0.0,539.95, 333 .95,206.0 


0.0, 0.8,0. 0,0.0 
'0.0, 0.0, 0.0,0. 8 
0.0, 8.0, 0.0, 0.0 
0 . 0 , 0 . 0 , 0 . 0  ,  @  .  0 

BEAM  DATA  CARDS 

1 ,1 ,2,7,1 ,1 ,0,8,0,0,000000,1 
5, 5, 6, 7,1 ,1 ,0,0,0,0,000000,1 

CONCENTRATED  MASS  ON  FREE  END 

6,0 ,2.591  ,2.591  ,2.591  ,0 .8 ,0 .8  ,8  .0 
0,0,0.000,0.008,0.000,0.0,0.0,0.0 

DYNAMIC  ANALYSIS  CARD  (IN  ACCORDANCE  WITH  REF . 4 

0,0, 0,0. 0,0. 0,0, 0.0 

1  ,0 

1  ,1  ,1  ,1  ,1  ,1  ,1 

END  OF  DATA  (DUMMY  TITLE  CARD) 

0,0, 0,0, 0,0, 0,0, 0,0 

The  output  -from  the  SAP IV  code  is  not  included  here  but 
all  the  essential  i  n  -forma  t  i  on  is  included  in  the  USLOB  code 
ou  tpu  t  . 

B.  PI  CRUST  INPUT  DATA  FOR  CASE  I 

The  PI  CRUST  code  is  used  to  process  the  SAP  IV  output 
and  calculate  influence  coefficients  corresponding  to  the 
base  or  support  motions  of  the  substructure.  Additionally, 
connectivity  between  the  hull  and  foundation  is  specified 
and  a  substructure  file  is  produced. 

The  following  input  fi'e  was  used  for  CASE  I: 

LIST  OPTIONS 

1  ,1  ,1  ,1  ,0,1  ,1  ,1  ,0,0 
0,1  ,1  ,1  ,1  ,1  ,0 , 1  ,0 ,0 


MODAL  SPECIFICATION 


k 

!• 

► 

► 

F 


5 

I 


5,0,0 

CONNECTIVITY 

1 ,5,10 
270 

MODES  FROM  SAP IV  INCLUDED  IN  ANALYSIS 
.  1 , 2 , 3 , 4 , 5 

END 

The  output  listing  -from  the  PI  CRUST  code  is  not 
included  here  but  all  essential  output  -from  these 
calculations  is  included  in  the  USLOB  output. 

C.  USLOB  INPUT  DATA  FOR  CASE  I  <1414“  STANDOFF) 

The  -following  input  data  was  suppl  i  ed  to  the  USLOB 
code  to  specify  the  parameters  -for  the  time  integration 
and  associated  output.  A  total  o-f  1201  time  steps  are 
specified  for  the  integration  interval,  each  covering  an 
interval  of  5.0E-6  seconds.  The  XXXX  parameters  denote 
parameters  for  eq.  5  which  are  ommitted  due  to  the 
sensitivity  of  this  equation. 

The  following  input  data  correspond  to  the  USLOB  output 
included  in  this  appendix. 

TIME  STEP  AND  CONTROL  CARD  FOR  INTEGRATION 
1201,1,1,9,11,1,1 

EQUATION  < 5)  PARAMETER  SPECIFICATION 

5.0E-6, XXXX, 1606,0 .0 
5000., XXXX , XXXX , XXXX , XXXX 
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0.  LISTING  -  USL08  OUTPUT  FOR  CASE  I 


CASE  I  --  WLVE  ON  FOUNDATION  VERT  OISPL 


CONTROL  PARAMETERS  FOR  SUBSTRUCTURE  1 


NPOINT  = 

7 

NOQF  = 

18 

NCDOF  * 

3 

NUDOF  = 

IS 

NFREQ  * 

10 

NFACE  = 

1 

NSYMTY  = 

0 

NFORCE  = 

13 

COWECT1QN  DATA  FOR  INTERFACE  POINTS 

SHL 

SUB 

X 

Y 

2 

XX  YY 

22  LSGPT  LF1ND 

DE6CQN 

54 

1 

1 

1 

1 

8  8 

8  5818  27 

2.780000E+82 

0DEGROT 

=  9.808080E+01 

1 EQUATION  NUMBERS  FOR  SUBSTRUCTURE 

1 

POINT 

X 

Y 

z 

XX 

YY 

22 

1 

1 

2 

3 

0 

0 

8 

2 

4 

5 

4 

0 

0 

0 

3 

7 

8 

9 

0 

0 

0 

4 

10 

11 

12 

0 

8 

8 

5 

13 

14 

15 

8 

8 

0 

4 

14 

17 

18 

8 

8 

0 

7 

0 

0 

0 

8 

0 

0 

FIXED-BASE  WTUIWL  FREQUENCIES  FOR  SUBSTRUCTURE  1 

MODE  FREQ  MODE  FREQ  NODE  FREQ  NODE  FREQ 

1  4.45541E+02  2  1.85485E+03  3  1.34307E+03  4  4.03705E403 

5  7.48871E+83  4  9.55805E+83  7  1 .04549E+04  8  1.21496E+04 

9  1.22840E404  10  1.88142E+84 

1 COMPACTED  FIXED-BASE  NODES  FOR  SUBSTRUCTURE  1 


I  2  3 

1  1.9139E-18  -1.2402E-01  9.5490E-10 

2  3.1548E-28  8.9711E-09  1.2482E-01 

3  -1.2482E-81  -4.8508E-20  3.7107E-11 

4  3.8184E-18  -2.4745E-01  1.9053E-09 

5  4.2987E-28  1.7900E-08  2.4745E-01 
4  -2.4745E-01  -9.4785E-20  5.9702E-U 

7  5.7052E-18  -3.4971E-01  2.8465E-B9 

8  9.4105E-20  2.4743E-88  3.4971E-81 


4 

6.2397E-12 
1 .1 J83E-12 
-1.4719E+88 
1.8034E-I1 
1.8084E-12 
-2.3485E+08 
9.988BE-12 
1.781IE-I2 


2.4775E-11 
1.4989E-13 
-2.3785E+80 
3.9845E-11 
2.7333E-13 
-1 . 4497E+80 
3.9373E-11 
2.4992E-13 
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'  1  V 


9  -3.4971E-81  -1.4448E-19  5.8957E-11  -2.3393E+08  1.4839E+88 

li  7.5444E-18  -4.9828E-01  3.7743E-09  5.9893E-12  2.3492E-11 

11  1.2478E-19  3.5459E-88  4.9028E-01  1.8477E-12  1.4183E-13 

12  -4.9028E-01  -1.9173E-19  3.5172E-11  -1.3957E+88  2.3572E+08 

13  9.3888E-18  -4.8834E-81  4.4848E-89  -3.9555E-13  -1.57I8E-12 

14  1.5485E-19  4.4884E-88  4.8834E-01  -5.9949E-14  -1.8775E-14 

15  -4.0834E-01  -2.3794E-19  -2.3534E-12  9.3383E-82  -4.2342E-82 

4  7  8  9  18 

1  1.4719E+88  3.5428E-14  1.2133E-28  1.2095E-12  -2.3785E+88 

2  1.3424E-13  2.2993E-11  -1.4719E+88  6.4535E-13  9.7987E-11 

3  2.4183E-18  -2.3485E+88  -3.3444E-12  -1.4434E+88  -1.8812E-12 

4  2.3485E+00  5.3848E-14  1.5884E-28  -7.4588E-13  -1.4497E+88 

5  2.1925E-13  3.4999E-11  -2.3485E+88  1.8488E-12  5.7223E-11 

4  -1.5814E-18  1.4719E+8I  5.4581E-12  2.3493E+88  -3.7449E-12 

7  2.3393E+88  4.9358E-14  1.1252E-28  -4.4834E-12  1.4839E+88 

8  2.1455E-13  3.6543E-U  -2.3393E+80  1.0514E-12  -4.8350E-11 

9  -1.4849E-18  1.4538E+88  -5.4572E-12  -2.3724E+88  -5.4847E-12 

18  1.3957E+88  2.7551E-14  4.4918E-21  -4.5984E-12  2.3572E+08 

11  1.2928E-13  2.1883E-11  -1.3957E+08  4.2577E-13  -1.8783E-10 

13  2.4274E-18  -2.3754E+88  3.3858E-12  1.4715E+88  -7.4239E-12 

14  -9.3383E-82  -1.9932E-15  -4.7535E-22  -5.2983E-12  -4.2342E-02 

15  -8.4444E-15  -1.4588E-12  9.3383E-82  -4.1994E-14  -7.7491E-12 

14  -2.2918E-12  2.2428E-82  -2.3883E-14  -1.8835E-02  -9.2143E-12 


1 SUBSTRUCTURE  CONTROL  PARAMETERS 


SUB  NFREQ  NDOF  NCOOF  NPTS  LOCO  NEBB  NBLOCK  NROUFC 

1  18  18  3  7  518  15  1  3 

N8LF0R  KOLSAR  NBL8AR  BE6FRO  8EGSSL  BE6BAR  8E8FOR 
1  254  1  1  1  1  1 

CONTROL  PARAMETERS  FOR  SOLUTION 


NT1ME  = 

1281 

NEQS  = 

438 

NSKIP  = 

2 

NRECS  = 

481 

NCHRG  * 

1 

NTHETA  = 

17 

NQUAO  * 

9 

NOELTA  = 

21 

NFINE  = 

11 

N800T  = 

244 

KOUPLE  * 

2 

NSUBS  = 

1 

SOLUTION  PARAMETERS 

TOTTALM  »  3 .541 492E+84 
RHOFLU  =  9.452148E-05 
CELT  »  5.080088E-04 
XLOAO  *  8 . 844880E+82 
RLOAO  =  1  .404008E+03 


AREA  3  3.439844E+04 
CFLU  *  4.080800E+04 
TEND  =  4.088088E-83 
STOFF  »  1  .414888E+03 
ORNRAO  =  8. 
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BATA  FOR  EMPIRICAL  PRESSURE-TIME  HISTORY 


WCHRG  =  5.088080E403  SURCUT  >  XXOOOOOOOOCX 
PZMLT  *  XXXXXXXXXXXX  TH1LT  »  XXXXXX50000(X 
P2EXP  *  ccoooococca  THEXP  =  xxxxxxxxxxxx 

1 VELOCITY  DISTRIBUTION  AT  SHELL  STATION  56  < SEGMENT  5  POINT  18) 


STEP  TIME  UDOT(  8)  «DOT(  98)  WOOT(IBI)  UD0T(278) 


1  8.  8.  8.  8.  8. 

3  1.8888E-85  4.9444E-01  7.5444E-02  1.7439E-81  7.5445E-82 

5  2.8888E-85  1.8855E+08  2.8722E-01  6.5444E-01  2.8728E-81 

7  3.8888E-85  3.8391E+88  5.8428E-81  1 .3805E+88  5.8418E-81 

9  4.8888E-85  6.1484E+8B  9.3465E-01  2.B219E+88  9.3436E-81 

11  5.8888E-05  8.6881E400  1.3193E+88  2.7669E+80  1.3136E408 

13  4.8000E-05  1 .1347E+01  1.7214E+88  3.4983E+88  1.7202E+00 

15  7.8888E-05  1.414X401  2.1416E400  4.227X408  2.1397E+08 

17  8.8I88E-05  1.7083E+01  2.5792E400  4.9559E+B8  2.5763E+88 

19  9.8888E-85  2.8102E+81  3.8245E+08  5.6665E400  3.8224E+88 

21  1.8888E-04  2.3241E+01  3.4895E+80  6.3713E+0B  3.4840E+08 

23  1.1888E-84  2.6547E+81  3.9697E+00  7.8994E+88  3.942X400 

25  1.2888E-04  2.9933E+01  4.4593E400  7.8237E+00  4.4502E+80 

27  1 .3880E-B4  3.3481E+01  4.9651E+88  8.5656E400  4.9539E+08 

2?  1.4888E-84  3.7123E+81  5.4786E+88  9.3141E+88  5.4450E40B 

31  1.588BE-84  4.8922E+81  6.8854E+00  1.8891E+81  5.9891E+88 

33  1.680BE-04  4.4771E+81  4.5359E+88  1.0856E401  6.5167E+88 

35  1.788BE-84  4.874X401  7.0707E+88  1.1647E401  7.B482E+88 

37  1.8888E-I4  5.2721E+01  7.6848E«88  1.24UE+81  7.5788E+80 

39  1.9B88E-84  5.4819E401  8.1394E+88  1.32I8E+81  8.1897E*88 

41  2.B88BE-B4  6.8891E+81  8.6698E+88  1.3957E+81  8.4353E+08 

43  2.1888E-84  4.5836E401  9.189X408  1. 472X401  9.1519E+88 

45  2.2888E-04  6.9161E401  9.7844E+88  1.5459E+81  9.441X408 

47  2.3888E-84  7.3292E+81  1.8289E401  1.618X40]  1 .0161E401 

49  2.4888E-84  7.7437E401  1.8788E401  1.4899E401  1.8648E401 

51  2.5888E-84  8.1493E401  1.U88E401  1.7569E401  1.1122E401 

53  2.6008E-04  8.5686E401  1.1443E401  1.8247E401  1.1580E401 

55  2.7808E-04  8.9672E401  1  .2897E401  1.889X401  1.2027E401 

57  2.8088E-84  9.3669E401  1.2528E401  1.958X40}  1.2451E40! 

59  2.9000E-84  9.7634E401  1.291X401  2.0106E401  1.2835E481 

61  3.88IIE-04  1.B154E+02  1.3381E401  2.846X401  1  .3289E401 

63  3.1080E-84  1.0542E402  1.366X401  2.1219E401  1.3578E401 

65  3.2888E-04  1.0933E402  1.402X401  2.1774E401  1.3918E401 

67  3.3800E-04  1.1317E402  1.437X401  2.2380E401  1.4261E401 

69  3.4H8E-04  1.1781E402  1.4712E401  2.2828E401  1.4585E401 

71  3.5888E-04  1.2884E402  1.5834E401  2.3353E401  1.4896E401 

73  3.6800E-84  1.2445E402  1.5354E401  2.385X401  1.5285E401 

75  3.7088E-84  1.2844E402  1.5655E401  2.4379E401  1.5494E401 


77  3.8088E-04  1.3224E402  1.5948E401  2.4892E40!  1.574X401 


79  3.90B0E-84  1.3599E+02  1.4222E+B1  2.5382E+01  1.4834E401 
81  4.0888E-04  1.3974E+82  1.4490E+81  2.5884E401  1.6291E+01 
83  4.1898E-94  1.4349E+02  1.4740E401  2.4395E+01  1.4524E+01 
83  4.2888E-04  1.47I9E402  1.4987E401  2.4885E+01  1.4759E+81 
87  4.3900E-04  1.5084E+02  1.7222E401  2.7371E+81  1.4979E+01 
89  4.4999E-94  1.5454E+82  1.7432E+01  2.7844E+B1  1.7174E+01 
91  4.5008E-84  1.5817E+82  1.7639E+01  2.8350E+01  1.7345E+01 
93  4.4888E-04  1.4175E+02  1.7843E+01  2.8822E401  1.7553E+81 
95  4.7008E-84  1.4533E+02  1.8817E401  2.9304E401  1.771IE+I1 
97  4.8898E-04  1.4888E+02  1.8188E+B1  2.9778E+81  1.7854E+01 
99  4.9888E-84  1.7237E+02  1 .8342E+B1  3.8233E+01  1.8880E+81 
181  5.8888E-04  1.7588E+02  1.8487E+81  3.8477E+81  1.8127E401 
103  5 . 1088E-84  1.7920E+82  1.84B1E+01  3.1125E+81  1.8223E+01 
105  5.2800E-04  1.8255E+02  1.8714E+01  3.1558E+91  1.8317E+81 
107  3.3800E-04  1.8585E+02  1.8824E+01  3.I964E+81  1.8488E+01 
189  5.4888E~04  1 .8910E+02  1.8900E+01  3.2371E+01  1.8443E+81 
111  5.5000E-04  1.9229E+82  1.8957E+01  3.2765E+01  1.8501E+01 
113  5.4I88E-04  1.9540E+02  1.9013E+01  3.3134E+01  1 .8537E+01 
115  5.7088E-84  1.9B45E+02  1.9043E+01  3.3489E+01  1.8545E+01 
117  5.80B0E-04  2.8147E+02  1.9*7114*1  3.3847E+81  1.8552E+81 
119  5.9000E-04  2.0442E+02  1.9074E+01  3.4182E+01  1.8533E+01 
121  4.8888E-84  2.073BE+02  1 .9*7<SE+*1  3. 4494£*01  1.8512E+01 
123  4.1000E-84  2.1018E+B2  1.9I57E+B1  3.4789E+01  1.8471E+01 
125  4.288BE-04  2.1283E+02  1.8995E+I1  3.5874E+01  1.8384E+01 
127  4.3800E-84  2.1558E+02  1.8932E+01  3.5335E+01  1.8300E+B1 
129  4.4088E-04  2.1811E+82  1.8849E+01  3.5583E+81  1 .8213E«01 
131  4.5BB8E-04  2.2067E+02  1.8790E+81  3.5823E+01  1.8110E401 
133  4.4008E-04  2.2317E+02  1.8481E481  3.4I54E481  1.7978E+01 
135  4.7080E-04  2.2543E+02  1.8572E401  3.42771401  1.7844E401 
137  4.8800E-04  2.2885E+02  1.8443E481  3.4484E401  1.7711E401 
139  4.90B0E-84  2.3041E+B2  1.8340E40]  3.4482E401  1 .7544E+0J 
141  7.0080E-84  2.3271E+82  1.8189E401  3.4873E401  1.7389E401 
143  7.1888E-84  2.3497E+*2  1.8835E401  3.7054E401  1.7211E40! 
145  7.2888E-04  2.3719E+02  1.7882E40I  3.7222E+81  1.7834E401 
147  7.3800E-04  2.3934E402  1 .7723E+01  3.7379E401  1.4851E481 
149  7.4088E-04  2.4147E402  1.7531E401  3.7534E401  1.4435E401 
151  7.5888E-04  2.4354E+02  1.7334E+01  3.7474E401  1.4415E401 
153  7.4080E-04  2.4558E+82  1.7144E401  3.7804E401  1.4201E401 
155  7.7088E-04  2.4757E+92  1.4941E401  3.7924E401  1.5993E401 
157  7.8808E-04  2.4950E+02  1.4734E+01  3.8842E401  1 .5744E+81 
159  7.9080E-04  2.5139E402  1.4504E401  3.8147E401  1.5489E401 
141  8.0080E-84  2.5325E482  1.4279E+81  3.8242E401  1.5248E+81 
143  8.1008E-04  2.5508E+02  1.4842E40!  3.8324E+01  1.5801E401 
145  8.2808E-04  2.5485E+02  1.5889E40]  3.8407E401  1.4725E401 
147  8.3008E-04  2.5857E+02  1.5544E401  3.8481E401  1.4448E401 
149  8.4880E-04  2.4824E402  1.5298E401  3.8542E401  1.4141E401 
171  8.5080E-84  2.4193E+02  1.5839E401  3.8594E401  1.3888E401 
173  8.4800E-84  2.4355E+02  1.4778E401  3.8448E+01  1.3404E+01 
175  8.7888E-84  2.4518E+02  1.4488E+01  3.8488E401  1.3294E401 
177  8.8880E-94  2.4443E+02  1.4200E401  3.8710E+01  1 .2985E+01 


179  8.9888E-04  2.4813E+02  1.3917E401  3.8730E+01  1.2482E401 


188 


181 

9.0808E-84 

2.6961E+02 

1.3643E401 

3.8739E+B1 

1 . 2383E+0 1 

183 

9.1008E-84 

2 . 71 82E+02 

1 .3335E+01 

3.8746E401 

1.2060E401 

18S 

9.2880E-04 

2.7238E+02 

1.3819E401 

3.8746E401 

1 . 1 725E+01 

187 

9.3888E-04 

2.7372E+82 

1.2709E401 

3.8734E401 

1.1397E401 

189 

9.4008E-04 

2.7505E+02 

1.2486E+01 

3.873 3E401 

1 .1077E+01 

191 

9.5000E-04 

2.7633E+02 

1.21I8E401 

3.8686E401 

1.8754E401 

193 

9.4088E-04 

2.7755E+82 

1.1765E+01 

3.8652E481 

1.0403E+01 

195 

9.7000E-84 

2.7873E+02 

1.1429E+81 

3.8608E401 

1.8052E401 

197 

9.8000E-04 

2.7991E+82 

1.1188E40] 

3.8552E+01 

9. 7897E408 

199 

9.9008E-04 

2.8186E+02 

1.8778E401 

3.8489E+01 

9.3739E40B 

281 

1 .8888E-83 

2.8216E+02 

1.8442E401 

3.8421E401 

9.0253E400 

203 

1.8108E-83 

2 . 8328E+82 

1.0089E401 

3.8346E+01 

8.6684E400 

285 

1.0200E-03 

2.8423E+82 

9.7388E+08 

3.8257E+01 

8.2991E+08 

287 

1 .030BE-B3 

2 . 8524E+82 

9.3938E+08 

3.8163E481 

7.9440E+00 

209 

1.0480E-83 

2.8623E+02 

9.8548E+08 

3.8861E+01 

7.5959E+08 

211 

1 .0500E-83 

2.8715E+02 

8.7064E408 

3.7956E+81 

7.2392E+88 

213 

1.0680E-83 

2.8882E+02 

8.3494E408 

3.7849E40J 

6.8749E+00 

215 

1 .0780E-03 

2.8887E+02 

7.9933E+08 

3.7734E401 

6.5125E400 

217 

1.0888E-83 

2.8971E+02 

7.6462E+00 

3.7611E401 

6. 15996400 

219 

1.0900E-03 

2.9054E+02 

7.3036E400 

3.7488E401 

5.8128E408 

221 

1.1088E-03 

2.9138E+02 

6.9413E+00 

3.7345E481 

5.44686+00 

223 

1 . 1 108E-03 

2 .9201E+92 

6.5667E+80 

3.7286E401 

5.8696E+08 

225 

1.120IE-83 

2.9272E+02 

6.1998E+00 

3.7058E+01 

4.70096+00 

227 

1 . 1 388E-03 

2.9342E+82 

S.8414E408 

3.6898E401 

4.3419E+88 

229 

1.1408E-03 

2.9411E+02 

5.4887E400 

3.6729E+01 

3.9896E+80 

231 

1 . 1 500E-83 

2.9473E+02 

5.1168E400 

3.6555E401 

3.6193E+00 

233 

1 .16B0E-03 

2.9531E+02 

4.7374E+88 

3.6375E40I 

3.2426E+8B 

235 

1.1 700E-83 

2.9589E+02 

4.3469E+00 

3.6185E401 

2.8762E+00 

237 

1 . 1800E-03 

2.9646E+02 

4.0038E408 

3.5985E401 

2.51836+08 

239 

1.1900E-03 

2.9703E+02 

3.6475E400 

3.5777E+8! 

2.16846+80 

241 

1.20B0E-03 

2.9752E+02 

3.27UE400 

3 . 5571E+8 1 

1.7998E+00 

243 

1 . 2108E-03 

2.9796E+02 

2.891 1E400 

3.5359E401 

1 .4288E+08 

245 

1.2200E-03 

2.9841E+02 

2.5193E+00 

3.5138E+01 

1.8673E+0B 

247 

1.2308E-83 

2.9886E+02 

2.1528E408 

3.4909E+01 

7.1232E-81 

249 

1.2400E-03 

2.9929E+82 

1 .7896E+00 

3.4676E401 

3.6192E-01 

251 

1.2580E-03 

2.9966E402 

1 .4168E+00 

3.4442E+81 

3.2035E-03 

253 

1.2608E-03 

2.9999E+02 

1.8388E400 

3.4205E+01 

-3.5943E-01 

255 

1 .2780E-B3 

3.0831E402 

6.6715E-01 

3.3959E401 

-7.1455E-01 

257 

1.2880E-03 

3.8063E402 

2.9858E-81 

3.3787E4B1 

-1 .0653E+08 

259 

1.2900E-83 

3.8894E402 

-6.5765E-82 

3.3451E401 

-1 .41076+88 

261 

1.3008E-03 

3.0121E+02 

-4.3242E-81 

3.3190E401 

-1 . 7571E+00 

263 

1.31B0E-I3 

3.0142E402 

-8.0686E-81 

3 . 2925E+0 1 

-2.1 101E+08 

265 

1.3200E-03 

3.0164E+02 

-1.1766E408 

3.2647E401 

-2.4572E+88 

267 

1.338BE-83 

3.8185E402 

-1.5432E408 

3 . 2363E+8 1 

•2.8808E+00 

269 

1 . 3400E-03 

3.8286E402 

-1.9846E408 

3.2070E+01 

-3.1364E+88 

271 

1 .3580E-03 

3.8225E402 

-2.2646E400 

3.1774E401 

-3.4702E+00 

273 

1 . 3600E-03 

3.0236E+02 

-2.4320E400 

3.1478E+01 

-3.8104E+00 

275 

1.3708E-03 

3.8246E40? 

-2.9962E400 

3.1175E401 

-4.1462E+88 

277 

1.38B8E-I3 

3.82S6E402 

-3.3578E+00 

3.8869E401 

-4.4775E+00 

279 

1.3900E-83 

3.8267E402 

-3.7125E408 

3.8553E401 

-4.8023E+80 

281  1.400BE-83  3.0278E+02  -4.0645E+00  3.0233E+01  -5.1225E+00 


283  1.4180E-03  3.0288E+02  -4.4214E+80  2.9917E+81  -5.4443E+8I 

285  M200E-83  3.B288E+82  -4.7754E+88  2.9599E*81  -5.7882E*88 

287  1.4308E-03  3.0279E+02  -5.12A2E+00  2.9278E+81  -8.0813E+00 

289  1.4488E-83  3.8288E+02  -5.4737E+88  2.8944E+01  -6.3929E+80 

291  1.4588E-83  3.0281E+02  -5.8180E+00  2.8808E+01  -6.69B2E+00 

293  1.4888E-83  3.0277E+02  -4.1811E+00  2.8268E401  -7.9821E^09 

295  1.47B8E-83  3.8287E+02  -8.5818E+08  2.7936E+01  -7.3024E+00 

297  1.4888E-03  3.8257E+82  -8.8409E*88  2.7599E+81  -7.5999E*88 

299  1.4988E-83  3.8247E*02  -7.1778E*08  2.7255E+01  -7. 8941 E* 00 

381  1.5888E-83  3.0239E+02  -7.5121E+88  2.8903E+01  -8.1858E+00 

<  This  history  is  abbreviated  tor  illustrative  purposes.) 

! VELOCITY  DISTRIBUTION  ALONG  SUBSTRUCTURE  1 


STEP 

1 

3 

5 

7 

9 

11 

13 

IS 

17 

19 

21 

23 

25 

27 

29 

31 

33 

35 

37 

39 

41 

43 

45 

47 

49 

51 

53 

55 

57 

59 

81 

83 

85 

87 

89 


TINE 

8. 

1.8888E-85 
2.8888E-85 
3.8S8BE-85 
4.0800E-05 
5.8888E-S5 
8.8888E-8S 
7.8888E-85 
8.8888E-85 
9.008BE-85 
1.8888E-84 
1.1088E-84 
1.288BE-04 
1.3888E-04 
1.4088E-84 
1.5000E-04 
1.8800E-04 
1.7808E-84 
1 .8B8BE-84 
1.9800E-84 
2.8888E-04 
2.1000E-04 
2.2888E-84 
2.3888E-84 
2.4888E-84 
2.5888E-84 
2.8888E-84 
2.7808E-84 
2.8888E-84 
2.9888E-84 
3.0080E-04 
3.1800E-84 
3.2888E-04 
3.3880E-04 
3.4888E-84 


VEL<  1,1) 

8. 

-7.8928E-84 

•3.8187E-83 

-8.8848E-83 

-1.8429E-82 

-2.8541E-82 

-3.8728E-82 

•5.2B09E-82 

-8.8819E-82 

-8.8191E-82 

-1.0522E-81 

•1.2514E-81 

•1.4552E-81 

■1.8745E-81 

-1.9827E-81 

-2.1381E-81 

-2.3783E-81 

-2.8205E-81 

-2.8885E-81 

-3.1149E-81 

•3.3881E-81 

-3.8129E-81 

-3.8889E-01 

•4.1829E-81 

-4.3351E-81 

■4.5582E-01 

•4.7787E-81 

-4.9885E-01 

-5.1832E-81 

-S.3828E-81 

•5.5399E-81 

•5.7835E-81 

-5.8558E-01 

■5.9947E-81 

■8.1224E-81 


VEL(  1,2) 

8. 

-4.4088E-03 
-2.0700E-02 
-3.3248E-82 
■3.7845E-82 
-2.8888E-02 
3.8282E-04 
4.4851E-02 
1.8789E-81 
1.883SE-81 
2.9143E-81 
4.1258E-81 
5.5122E-81 
7.1858E-01 
8.8589E-81 
1.8795E+80 
1 .2921E+B8 
1.5173E+88 
1 .7828E+88 
2.9222E+00 
2.2978E+88 
2.5879E+08 
2.8949E+88 
3.2141E+00 
3.5481E+88 
3.8942E+88 
4.2550E+00 
4.8285E+80 
5.8121E+8# 
5.4117E+88 
5 . 8235E+00 
8.2448E+88 
8.8781E+80 
7.1191E480 
7.5894E+00 


VEL< 

8. 

7.5445E-82 
2.8728E-01 
5.8418E-81 
9.3438E-81 
1.3188E+B8 
1.7282E+80 
2.1397E+B8 
2 . 5783E+88 
3.8224E+00 
3.4840E+80 
3.9825E*8B 
4 . 4582E+00 
4.9539E+00 
5.4850E+00 
5.9891E+00 
8.5187E+00 
7.8482E+00 
7.5788E+08 
8.1097E+80 
8.8353E+08 
9 . 151 9E+08 
9.8818E+08 
1 .8181E+01 
1 .8848E+81 
1.11 22E+0 1 
1 . 1 580E+0 1 
1.2027E+81 
1.2451E+81 
1 .2835E+01 
1 . 3209E+0 1 
1 .3570E+81 
1.3918E+81 
1.4281E+01 
1 .4585E+01 


1,3)  VEL<  8,1) 


8.8872E-07 
1.9828E-88 
1 .8589E-08 
-3.5152E-08 
—4 . 3943E-06 
-2.8589E-85 
-1.2187E-84 
-3.7888E-84 
-8.B828E-84 
-1.7841E-B3 
-3.188BE-03 
-5.1171E-83 
-7.7U8E-83 
-1.1817E-82 
-1.5095E-B2 
-1.9998E-82 
-2.5775E-82 
-3.2514E-02 
-4.8329E-82 
-4.9349E-82 
-5.9891E-02 
-7.1427E-82 
-8.4881E-02 
-9.9248E-02 
-1.1538E-81 
-1 .3297E-81 
-1.5208E-81 
-1.7288E-81 
-1 .9473E-01 
-2.1829E-81 
-2.4339E-81 
-2.7882E-81 
-2.9812E-81 
-3.2788E-81 


VEL( 

0. 

-1.9785E-85 
■3.3974E-B5 
1.2365E-05 
4.4999E-B5 
-1.8797E-B4 
-9.1328E-04 
-2.3448E-03 
-4.2851E-03 
-8.0883E-83 
-8.8775E-83 
-5.5941E-03 
-1.2343E-03 
7.0B35E-83 
2.0128E-82 
3.8730E-02 
8.3995E-02 
9.7299E-82 
1 .4B38E-01 
1 .9492E-B1 
2.8285E-01 
3.4499E-81 
4.4315E-81 
5.5824E-81 
8.9154E-01 
8.4451E-B1 
1 .8188E+00 
1 .2181E+00 
1.4379E+88 
1.8854E+00 
1 .9593E*88 
2.2803E+08 
2 . 5898E+00 
2.9481E*00 
3 . 3327E+08 


8,2)  VEL<  8,3) 


-4.3588E-14 
1.4181E-12 
3.8135E-10 
1.3805E-08 
1.8281E-87 
1.5070E-08 
8.8295E-88 
3.7798E-05 
1.3459E-04 
4.0597E-04 
1 .0888E-03 
2.5841E-03 
5.3195E-03 
1.0372E-02 
1.8757E-82 
3.1738E-02 
5.0828E-02 
7.8805E-82 
1 . 1085E-0 1 
1.5337E-01 
2.0588E-01 
2.8573E-01 
3.3527E-01 
4.1357E-81 
5.B881E-01 
5.9853E-01 
7.B155E-01 
8.1594E-B1 
9.3985E-01 
1 , 0733E+80 
1 . 21 83E+00 
1 . 3885E+88 
1.5299E+88 
1 .7801E+00 


."v  r 

rr  '.t  ^  i' 

71 

3.5(((E-(4  -4.244(E-(1 

8 . 0  308E+0B 

1 .4894E+01  -3.5835E-01 

3.7494E+0B 

1.8791E488 

73 

3.488(E-(4  -4.3518E-81 

8.4948E+88 

1.5285E4B1  -3.9824E-81 

4.1949E488 

2.0478E488 

»* 

75 

3.7B08E-84  -4.4494E-81 

8.9484E+08 

1.5494E401  -4.2323E-B1 

4.4753E+88 

2.2439E408 

• 

77 

3.8((8E-(4  -4.5414E-81 

9.4509E+88 

1.5747E+81  -4.5715E-81 

5.1844E+B8 

2.4700E400 

» 

r 

79 

3.9888E-84  -4.4229E-81 

9.9375E+08 

1.4B34E4B1  -4.9189E-81 

5.7234E+00 

2.4857E+00 

81 

4,((((E-(4  -4.7834E-81 

1.B431E+01 

1 .4291E+81  -5.2734E-81 

4.2922E+90 

2.9112E408 

83 

4.1((8E-(4  -4.7754E-81 

1 .(931E+81 

1.4524E+01  -5.4347E-01 

4.8894E+88 

3.1470E408 

•  -  * 

85 

4.2888E-84  -4.8457E-81 

1 .1437E+81 

1.4759E4B1  -4.BB13E-81 

7.5155E+B8 

3.3937E+B0 

«*. «.» 

•  •  • , » 

87 

4.3888E-84  -4.9893E-B1 

1.1948E+B! 

1.4979E4B1  -4.3723E-81 

8.1493E40B 

3.4518E+80 

89 

4.4808E-84  -4.9747E-81 

1 .2445E+81 

li7174E*01  -6«74$4E-8j 

8 .8503E+00 

3.9219E+00 

1 

91 

4.5888E-84  -7.8374E-(1 

1.2989E+B1 

1.7345E+81  -7 . 1217E-81 

9.5574E4B8 

4.2043E+B8 

93 

4.4((8E-(4  -7.1828E-81 

1 ,35I2E+(1 

1 .7553E+B1  -7.4945E-81 

1 .B289E+B1 

4.4993E+88 

95 

4.7888E-94  -7.1422E-81 

1 .4843E+81 

1.7718E+B1  -7.8491E-B1 

1.K44E401 

4.8849E400 

97 

4,8(((E-(4  -7.2251E-81 

1 .458BE+B1 

1.7854E+81  -8.2381E-01 

1 . 1819E+B1 

5.1244E+08 

99 

4.9888E-84  -7.2891E-81 

1 . 51 1 8E+0 1 

1.8800E+01  -8.4822E-B1 

1.2414E+01 

5.4580E400 

*  * 

r 

1 

1(1 

5.8888E-84  -7.3577E-*1 

1.5459E+81 

1.8127E+B1  -8.9482E-B1 

1.3424E4B1 

5.8((2E4(( 

1(3 

5.1888E-84  -7.4298E-(1 

1.4288E4B1 

1.8223E+81  -9.3109E-B1 

1 . 4253E+0 1 

4. 1527E+80 

i. 

IN 

5.2888E-84  -7.5894E-81 

1.474BE401 

1 .8317E+01  -9.4534E-01 

1 .589X4(1 

4.S144E488 

1(7 

5.3(((E-84  -7.5834E-81 

1.7311E+01 

1 .8488E+01  -9.9872E-81 

1.5958E+81 

4.8856E+08 

1(9 

5.4888E-84  -7.4448E-81 

1.7871E401 

1.8443E+01  -1.0311E+08 

1.4814E+01 

7.2452E+8B 

111 

5.5(((E-(4  -7.7554E-81 

1.8433E+01 

1 .85(1E+(1  -1 .(423E+88 

1.7484E481 

7.4533E+B8 

V  1- 

n 

113 

5.489BE-84  -7.8511E-91 

1 .8995E+01 

1.8537E+01  -1.0923E+0( 

1.8542E4B1 

8.B494E408 

•r  :  *  "  * 

V 

115 

5.7(((E-(4  -7.9475E-81 

1 .9559E+01 

1.8545E+81  -1.1209E408 

1.9442E+B1 

8.4538E408 

♦ 

«  _  -  < 

117 

5.8(((E-(4  -8.(49(E-81 

2.B129E+01 

1 .8552E+81  -1 . 1488E+68 

2 . 0321 E+0 1 

8.8457E+80 

119 

5.98B8E-B4  -8.1542E-81 

2.0702E+01 

1 .8533E+81  -1.1735E+08 

2.1199E+81 

9.2848E488 

121 

4.BBBBE-(4  -8.2443E-(1 

2.1275E+01 

1.8512E+81  -1.1974E+M 

2.2074E+01 

9.7185E+0B 

\  * 

123 

4.1BB8E-B4  -8.38B4E-B1 

2.1851E4B1 

1.8471E481  -1.2197E+8B 

2,2944E*81 

1.8142E4B1 

M 

g 

125 

4.2B8BE-B4  -8.5824E-B1 

2.2434E+B1 

1 .8384E+B1  -1 .2483E+0B 

2.3804E+B1 

1.8579E401 

1 

127 

4.3M8E-04  -8.4285E-81 

2.3823E4I1 

1 .83BBE+81  -1.2592E4BB 

2.4459E481 

1.1821E4B1 

. 

129 

4.4(8(E-(4  -B.7428E-B1 

2.3488E+81 

1 . 8213E+81  -1.2744E+BB 

2 .5499E+81 

1.1448E+81 

* 

% 

* 

131 

4.58B8E-B4  -8.87(1E-(1 

2.4194E+01 

1.81ME+81  -1.2919E+88 

2.4325E401 

1.1919E+01 

133 

4.4BB8E-84  -8.9944E-B1 

2 . 4788E+B1 

1.7978E+81  -1.3857E4M 

2.7134E+01 

1 .2375E+81 

*  •/  *. 

135 

4.7(((E-(4  -9.1225E-(1 

2 . 5380E+B 1 

1.7844E481  -1 .31 74E+88 

2.7938E481 

1.2834E«(1 

i 

137 

4.88B8E-84  -9.2519E-B1 

2.5971E4B1 

1.7711E481  -1.3277E+M 

2.8784E+B1 

1 . 3302E+0 1 

139 

4.988IE-84  -9.3799E-81 

2.4542E401 

1.7544E401  -1.334BE4B8 

2.9444E401 

1.3774E401 

-  -  ■■ 

141 

7 .(((BE-84  -9.5B74E-B1 

2.7142E401 

1.7389E+(1  -1 .3425E+B8 

3.8281E401 

1.4258E+(1 

143 

7.1((8E-(4  -9.4374E-I1 

2.7748E+81 

1.7211£4(1  -1.3473E488 

3.89I4E4B1 

1 .4738E481 

s‘ 

145 

7.288BE-B4  -9.7783E-81 

2.8354E+8! 

1.7834E+B1  -1.3584E+0B 

3.1488E481 

1.5215E+01 

147 

7.3888E-84  -9.893(E-81 

2.8958E+01 

1 .485IE*(i  -1.3519E+(( 

3.2277E481 

1 .5783E+81 

g 

149 

7.48B8E-B4  -1.8B21E+BB 

2.9552E+B1 

1 . 4435E+8 1  -1.3518E488 

3.2922E+81 

1 .4193E+81 

151 

7.5(((E-(4  -l.(15(E+(8 

3.(155E4(1 

1 .4415E+B1  -1.3582E4B8 

3.3543E401 

1.4485E401 

*  - 

153 

7 .41f(E-94  -1 .(277E+(( 

3.(755E*01 

1.4201E+01  -1.3472E+BB 

3.4148E+81 

1.7177E+01 

155 

7.7(((E-(4  -1 .8394E+8® 

3. 1 352E+81 

1 .5993E+81  -1 .3429E+08 

3.4713E481 

1.7449E481 

157 

7.8(((E-(4  -1.8519E+B8 

3.1957E+B1 

1.5744E+B1  -1.3371E48B 

3.5243E+01 

1.8159E401 

159 

7.9(((E-(4  -l.(447E+(( 

3.2543E+01 

1.5489E40  -1.33B1E480 

3.5789E4B1 

1.8447E401 

i 

141 

8.(BB8E-(4  -i.(745E4(( 

3.3147E+01 

1.5248E+81  -1.3219E+88 

3 .4292E+81 

1.9131E+B1 

,  _  * 

143 

8.1B88E-B4  -1,(88JE+B( 

3.3744E4(] 

1.5M1E4B1  -1.3125E4BB 

3.4772E4(1 

1.9411E401 

' 

145 

8.2((BE-(4  -1.1((2E+B( 

3.4372E+B1 

1 .4725E+(1  -1.3(21E+(8 

3.7231E+81 

2.B885E401 

\  -* 

S  *  »  ■  _M 

147 

8.3((BE-(4  -1 .1121E+8B 

3.4988E+B1 

1.444(E*(1  -1.2989E4(( 

3.7449E401 

2.9S55E401 

A  « 

149 

8.4BBBE-B4  -1.1235E+(( 

3.5584E+01 

1 .4141E+81  -S.2789E4(( 

3 . 8888E+0 1 

2.K19E4B1 

,.-v: 

i 

171 

8.5B88E-B4  -l,1348E+(( 

3.4189E4B1 

1 .3888E+B1  -1 .2444E+08 

3.8490E401 

2.1477E481 

► 

1(3 

i 

'  *  *  • 

•V  .  .  .  ...  . 

•  •  *  «  1 

173  8.4888E-04  -1. 1448E400  3.4798E*81  1  .3686E+81  -1.2532E+88  3.8875E+81  2.1929E+81 
175  8.7888E-84  -1.1572E+88  3.7488E+81  1.3294E+81  -1.2397E+8B  3.9247E+81  2.2374E+83 
177  8.8888E-84  -1.1479E+88  3.88B8E<81  1.2985E+01  -1 .2257E+88  3.9685E+01  2.28I1E+81 
179  8.9888E-84  -1.1785E+08  3.8413E+81  1.2482E+81  -1.211 5E+88  3.9951E+81  2.3248E+01 
181  9.8888E-84  -1.1892E400  3.9213E+81  1.2388E+31  -1.1971E+08  4.8287E*01  2.3440E«81 
183  9.1888E-84  -1.1993E+00  3.9821E+81  1.2040E+01  -1.1824E+88  4.0615E+01  2.4870E+01 
185  9.2888E-04  -1.2890E+08  4.0438E+81  1.1725E+01  -1.1681E+88  4.0937E+81  2.4471E+B1 
187  9,3088E-04  -1.2193E+08  4.1838E+01  1.1397E+01  -1.1539E+80  4.1257E+81  2.4842E+81 
189  9.4088E-84  -1.2291E408  4.1641E+I1  1.1077E+01  -1.1399E+08  4.1575E+01  2.5245E+B1 
191  9.5808E-84  -1.2383E+00  4.2240E+01  1.0754E+01  -1.1244E+00  4.1894E+01  2.5617E+01 
193  9.4000E-04  -1.2473E+88  4.2849E+01  1.0483E+81  -1.U35E+88  4.2214E+81  2.5980E+81 
195  9.7088E-84  -1.2544E+08  4.3457E+01  1.0052E+01  -1 .1012E+00  4.2542E+01  2.6333E+81 
197  9.8888E-I4  -1.2455E+88  4.4842E+81  9.7897E+08  -1.8895E+08  4.2873E+8!  2.4674E401 
199  9.9888E-84  -1.2737E+88  4.4441E+81  9.3739E+00  -1.0784E+08  4.3213E+01  2.7087E+01 
281  1.8808E-03  -1.2818E+88  4.5263E+01  9.8253E+88  -1 .8486E400  4.3563E+01  2.7325E+81 
213  1.818BE-83  -1.2898E+08  4.5870E+81  8.4404E+88  -1 .0594E+08  4.3926E+81  2.7629E+81 
285  1.8288E-83  -1.2982E«88  4.4477E+81  8.2991E408  -1 .8512E400  4.4383E+01  2.7918E+01 
287  1.8388E-83  -1.3857E+88  4.7079E+81  7.9448E+88  -1.8441E+B8  4.4496E+01  2.8192E+01 
289  1 .8488E-83  -1.3129E+00  4.7474E+81  7.5959E+88  -1 .8381E+08  4.5!B5E<81  2.8448E«81 
211  1.858IE-83  -1.3288E+88  4.8279E+81  7.2392E+80  -1.8334E+80  4.5533E+81  2.8487E+81 
213  1.8688E-03  -1.3272E488  4.8885E+81  6.8749E408  -1.8388E«88  4.5988E*01  2.8988E+81 
215  1.0788E-03  -1.3338E+08  4.9498E+81  6.5125E+80  -1 .0280E+0B  4.6447E+81  2.9U1E+81 
217  1.0880E-83  -1.3399E+00  5.8089E+01  4.1599E+88  -1.8273E«00  4.6934E+0I  2.9294E«81 
219  1.0908E-83  -1.3448E*00  5.0684E+01  5.8128E+00  -1.B280E408  4.7447E+01  2.9458E+01 
221  1.1880E-03  -1.3528E+88  5.1285E+01  5.4448E400  -1.B381E408  4.7982E+B1  2.9402E+01 
223  1.1188E-03  -1.3579E+80  5.1889E+01  5.8494E+00  -1.0335E488  4.8541E+01  2.9727E+81 
225  1 . 1208E-83  -1.3433E+08  5.2498E401  4.7089E408  -1.0384E480  4.9124E401  2.9831E+81 
227  1.1380E-03  -1.3483E+80  5.3085E*01  4.3419E+8B  -1.B444E488  4.9731E+01  2.9914E+01 
229  1 , 1488E-03  -1.3730E+00  5.3474E*81  3.9896E+88  -1.0522E+08  5.0368E+81  2.9980E+81 
231  1.1500E-03  -1.3777E+80  5.4274E+01  3.4193E+08  -1.0613E+08  5.1014E+01  3.0024E+01 
233  1 . 1608E-B3  -1.3828E+08  5.4874E+01  3.2424E+08  -1 .8719E«08  5.1690E+01  3.0048E+81 
235  1.17B8E-83  -1.3868E408  5.5473E+01  2.8762E+08  -1.0838E+00  5.2389E+01  3.0050E+01 
237  1.188IE-03  -1.3900E+88  5.4044E+01  2.5183E+08  -1 .8971E+00  5.3111E*01  3.8831E>81 
239  1.1900E-83  -1.3934E489  5.4454E+01  2.1484E+88  -1.1118E+80  5.3853E+01  2.9992E+01 
241  1.2880E-83  -1.3949E408  5.7248E401  1.7998E+08  -1.1274E+00  5.4417E+01  2.9931E*01 
243  1.2100E-03  -1.4B00E400  5.7844E+01  1.4288E+80  -1.1445E+80  5.5398E+01  2.9858E+01 
245  1.2280E-03  -1.4038E+00  5.8437E401  1.8473E408  -1.1624E408  5.4197E+01  2.9749E+81 
247  1.238BE-03  -1.4858E+88  5.9825E+81  7.1232E-81  -I.1817E+00  5.7012E+81  2.9427E+01 
249  1.2400E-03  -1.4080E408  5.9484E481  3.4192E-81  -1.2017E*88  5.7848E+81  2.9484E+81 
251  1.2508E-03  -1.4180E+08  4.8192E+01  3.2835E-83  -1.2227E+B0  5.8681E+01  2.9328E+01 
253  1.2480E-03  -1.4119E400  4.0781E*81  -3.5943E-01  -1 .2445E+00  5.9533E+01  2.9134E+01 
255  1.2788E-83  -1.4137E408  4.1368E+01  -7.1455E-01  -1.2471E400  4.0395E+01  2.8925E401 
257  1.2880E-03  -1.4151E480  4.1958E401  -1 .8453E+00  -1 .2905E+88  4.1244E+01  2.8492E+01 
259  1 .2988E-03  -1.4140E+B8  4.2524E401  -1.4187E400  -1.3144E408  4.2138E401  2.8435E401 
241  1.3800E-03  *1.4144E408  4.3181E401  -1  .7571E+08  -1 .3387E400  4.3814E401  2.81S3E401 
243  1.3100E-83  -1.4174E480  6.3484E401  -2.1181E400  -1.3433E400  4.3890E+01  2.7845E+01 
245  1.3280E-03  -1.4178E408  4.4248E401  -2.4572E408  -1.3882E408  4.4745E401  2.7512E40I 
247  1.3388E-03  -1.4179E408  4.4844E401  -2.8B00E*00  -1.4131E408  4.5434E+01  2.7154E+01 
249  1.3408E-83  -1.4175E400  4.5419E+81  -3.1344E+08  -1.4388E+08  4.4501E40I  2.4772E401 
271  1.3S80E-I3  -1.4148E408  A.5998E+81  -3.4782E+08  -1.4629E+88  4.734BE401  2.4345E+01 


273  1.3480E-I3  -I.4143E408  4.4548E401  -3.81I4E40I  -1.4875E408  4.82I9E*81  2.5933E40I 


275  1.3788E-03  -1.415AE408  6.7146E+01 
277  1.3880E-03  -1.4142E+00  6.7718E+81 
279  1.3988E-03  -1.4128E400  A.8285E401 
281  1.4888E-03  -1.4189E+00  4.8845E+01 
283  1.4188E-03  -1.4091E+00  4.9415E+81 
285  1.4288E-03  -1.4070E+0B  6.9984E+01 
287  1.43ME-03  -1 .4044E+08  7.0549E+01 
289  1 .44SSE-03  -1.4818E+88  7.1118E+B1 
291  1.4588E-83  -1.3987E+80  7.1664E+81 
293  1 . 4600E-03  -1 . 395AE+00  7 . 2220E+0 1 
295  1.4700E-03  -1.3921E488  7.2788E+81 
297  1.4880E-83  -1.3884E+00  7.3338E+01 
299  1.4900E-03  -1.3844E+08  7.3892E+01 
301  1.5000E-03  -1.3801E+00  7.4440E+81 


-4.14A2E408  -I.5128E408  A.9847E481 
-4.4775E+08  -1.53A0E400  6.9870E401 
-4.8023E400  -1.559AE488  7.8A78E481 
-5.1225E+88  -1 .5825E+80  7.1469E401 
-5.44A3E408  -I.A04AE400  7.2241E+01 
-5.7A62E+00  -1.A259E400  7.2993E401 
-6.8813E+00  -1.6461E408  7.3724E+81 
-A.3920E400  -1.AA52E408  7.4433E+81 
-6.6982E+88  -1.A832E400  7.5118E481 
-7.0021E+I8  -1.A999E+00  7.5779E401 
-7.3824E408  -1.7154E40B  7.6414E481 
-7.5999E408  -1.7295E408  7.7023E+01 
-7.8941E+I0  -! .7421E+08  7.7685E401 
-8.1850E+08  -1.7533E+80  7.8161E401 


2.5478E401 
2.5088E401 
2.4498E401 
2.3973E+01 
2.342AE40! 
2.285AE401 
2.22A4E+01 
2.1652E+01 
2.1018E401 
2.0363E+81 
1.9A88E401 
1.8994E+01 
1 .8280E+81 
1.7546E+81 


(This  history  is  abbreviated  Tor  illustrative  purposes.) 


E. 


INTEGRATION  USING  SIMPS1 


The  -final  step  in  the  ELSHQK  analysis  of  Case  I  is 
use  an  appropriate  integration  algorithm  to  calculate  a 
deflection  history  from  the  velocity  information  provide 
by  the  USLOB  code.  Simpson's  1/3  rule  was  chosen  in  thi 
work  and  in  t.iis  example,  implemented  in  the  form  of 
S I MRS 1 .BAS  which  appears  in  appendix  B. 

The  following  is  an  entire  output  listing  from  this 
program.  In  the  listing,  the  coordinate  directions  are 
defined  as: 


X  =  fore  and  aft 

Y  =  vertical  or  normal  to  the  hull 
2  =  athwar  t  sh i ps 


T  = 

1.0100 

nsec 

DX  = 

0.000000 

in. 

DY  » 

0.000000 

in. 

02  = 

-8.008082 

T  * 

0.0300 

msec 

OX  = 

0.000000 

in. 

OY  = 

0.080881 

in. 

02  = 

-0.000014 

T  a 

0.0500 

nsec 

OX  = 

0.000001 

in. 

OY  * 

8.800001 

in. 

02  * 

-0.000840 

T  = 

0.0700 

nsec 

OX  = 

0.000002 

in. 

OY  * 

0.000000 

in. 

02  » 

-8.080863 

T  = 

0.0900 

nsec 

OX  = 

0.000004 

in. 

OY  * 

-0.000884 

in. 

02  a 

-0.888144 

T  = 

0.1100 

nsec 

OX  = 

0.000006 

in. 

OY  = 

-0.000012 

in. 

02  = 

-0.008223 

T  = 

0.1300 

nsec 

OX  = 

0.000009 

in. 

OY  * 

-0.000026 

in. 

02  a 

-0.000322 

7  = 

0.1500 

nsec 

OX  = 

0.000013 

in. 

OY  » 

-0.080047 

in. 

02  = 

-8.808441 

7  = 

0.1700 

nsec 

OX  = 

0.000018 

in. 

OY  a 

-0.000075 

in. 

02  a 

-0.880581 

7  = 

0.1900 

nsec 

ox  = 

0.000023 

in. 

OY  » 

-0.800112 

in. 

02  a 

-0.000741 

7  = 

0.2100 

nsec 

ox  = 

0.000029 

in. 

OY  a 

-0.800157 

in. 

02  * 

-0.000928 

7  = 

0  "''00 

nsec 

ox  = 

1.000036 

in. 

OY  a 

-0.080210 

in. 

02  = 

-8.001116 

7  = 

0.2500 

nsec 

ox  = 

0.000043 

in. 

OY  a 

-0.000271 

m. 

02  a 

-0.001329 

7  = 

0.2700 

nsec 

ox  = 

0.000050 

in. 

OY  a 

-0.000339 

in. 

02  a 

-0.001555 

7  * 

0.2900 

nsec 

ox  = 

0.800056 

in. 

OY  a 

-0.801414 

in. 

02  = 

-0.081793 

7  = 

0.3100 

nsec 

ox  = 

0.000863 

in. 

OY  = 

-0.808493 

in. 

02  = 

-0.082048 

7  * 

0.3300 

nsec 

ox  = 

0.880069 

in. 

OY  a 

-8.000577 

in. 

02  = 

-0.882295 

7  = 

0.3500 

nsec 

ox  = 

0.000874 

in. 

OY  = 

-0.000662 

in. 

02  = 

-0.802555 

7  = 

0.3700 

nsec 

OX  s 

0.808079 

in. 

DY  a 

-0.800748 

in. 

02  = 

-0.882819 

7  = 

0.3900 

nsec 

ox  * 

0.080882 

in. 

DY  = 

-0.080832 

in. 

02  = 

-0.003086 

7  = 

0.4100 

nsec 

DX  = 

0.808884 

in. 

OY  a 

-0.808913 

in. 

02  = 

-0.803354 

7  = 

0.4300 

nsec 

ox  = 

0.008885 

in. 

DY  a 

-0.000988 

in. 

02  = 
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OX  = 

-8.888185 

in. 

DY  =  -8.880313 

in, 

02  a 

-0.813987 

7  * 

2.9588 

nsec 

DX  = 

-8.888898 

in. 

OY  =  -8.088238 

in. 

02  = 

-8.014897 

7  = 

2.9788 

nsec 

DX  = 

-8.888891 

in. 

OYa  -8.880161 

in, 

02  = 

-8.014166 

* '  . 

7  = 

2.9988 

nsec 

OX  = 

-8.888884 

in. 

DY  a  -0.888885 

in. 

02  = 

-8.814196 

l  - 

7* 

3.8188 

nsec 

OX  = 

-8.808877 

in. 

OYa  -0.888810 

in. 

D2  = 

-0.014185 

7  = 

3.8388 

nsec 

DX  = 

-0.888878 

in. 

DYa  8.808861 

in. 

D2  = 

-0.014135 

7  = 

3.8588 

nsec 

OX  = 

-8.888863 

in. 

DY=  0.800126 

in. 

02  = 

-8.014046 

7  = 

3.8788 

nsec 

DX  = 

-0.880056 

in. 

DY  a  8.886184 

in. 

02  = 

-0.813917 

7  = 

3.8988 

nsec 

DX  = 

-8.888858 

in. 

DY  a  8.888233 

in. 

02  = 

-8.013751 

- 

7  = 

3.1188 

nsec 

OX  = 

-0.088844 

in. 

OY  a  0.088272 

in. 

02  a 

-8.813548 

i~-c_ 

7  = 

3.1388 

nsec 

OX  = 

-8.808839 

in. 

DYa  8.886299 

in. 

02  a 

-8.813309 

Lr,~ 

7  = 

3.1588 

nsec 

OX  = 

-0.888835 

in. 

DYa  0.888315 

in. 

02  = 

-0.813834 

V*  /' 

7  = 

3.1788 

nsec 

OX  = 

-8.880031 

in. 

OYa  0.000328 

in. 

02  = 

-0.812726 

*  „*«.  <  * 

7  = 

3.1988 

nsec 

OX  a 

-8.808028 

in. 

DY  a  0.808312 

in. 

02  a 

-0.812385 

V 

7  = 

3.2188 

nsec 

DX  = 

-0.888027 

in. 

DY  =  0.080294 

in. 

02  a 

-0.812011 

7  = 

3.2388 

nsec 

DX  = 

-0.808826 

in. 

DYa  0,008264 

in. 

02  = 

-0.811688 

***** 

7  = 

3.2588 

nsec 

DX  = 

-8.888826 

in. 

OYa  8.888225 

in. 

02  = 

-0.011176 

7  = 

3.2788 

nsec 

DX  = 

-8.888827 

in. 

DY  =  0.008177 

in. 

02  = 

-8.018715 

7  = 

3.2988 

nsec 

DX  = 

-8.088838 

in. 

OY  =  0.006123 

in. 

02  = 

-0.810231 

7  = 

3.3188 

nsec 

DX  = 

-8.808833 

in. 

DYa  0,000863 

in. 

02  = 

-0.889723 

\  • 

7  = 

3.3388 

nsec 

OX  a 

-0.008837 

in. 

OY  a  8.088888 

in. 

02  = 

-8.089192 

• 

7  = 

3.3588 

nsec 

DX  a 

-8.888842 

in. 

OY  =  -0.880064 

in. 

D2  a 

-8.808642 

7  = 

3.3788 

nsec 

DX  a 

-8.888847 

in. 

DY  a  -8.880127 

in. 

02  = 

-0.808876 

- r. 

7  = 

3.3988 

nsec 

OX  a 

-8.088854 

in. 

DY  »  -0.800188 

in. 

02  = 

-0.807492 

7  = 

3.4188 

nsec 

OX  = 

-0.888860 

in. 

OY  =  -B. 888244 

in. 

02  = 

-0.886897 

•!*;  / 

7  = 

3.4388 

nsec 

DX  = 

-0.808868 

in. 

DY  =  -0.800294 

in. 

02  = 

-6.006290 

7  = 

3.4588 

nsec 

DX  = 

-8.000075 

in. 

OY  =  -0.088336 

in. 

02  = 

-0.885673 

7  = 

3.4788 

nsec 

OX  = 

-8.088883 

in. 

DY  a  -8.086368 

in. 

02  = 

-6.085856 

7  = 

3.4988 

nsec 

OX  = 

-0.800891 

in. 

OY*  -0.008398 

in. 

02  = 

-8.884422 

- 

7  = 

3.5188 

nsec 

OX  = 

-0.088898 

in. 

DY  =  -8.888480 

in. 

02  = 

-8.083791 

7  a 

3.5388 

nsec 

DX  » 

-0.808105 

in. 

DY»  -8.008399 

in. 

D2  » 

-0.083159 

7  * 

3.5588 

nsec 

DX  = 

-0.888112 

in. 

DY  »  -0.008387 

in. 

02  = 

-8.662538 

■  . 

7  » 

3.5788 

nsec 

DX  = 

-8.808119 

in. 

OY  »  -0.888362 

in. 

02  a 

-6.061965 

■V  *. 
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i 

. 

I 

T  * 

3. 5988 

nsec 

OX  = 

•0.080125 

in. 

OY  =  -0.880327 

in. 

02  = 

-0.001286 

• 

T  = 

3.6100 

nsec 

OX  = 

-6.008138 

in. 

DY=  -0.888281 

in. 

02  = 

-0.888677 

.  T  = 

3.6300 

nsec 

ox  = 

-8.808134 

in. 

DY=  -8.880227 

in. 

02  = 

-8.800079 

T  = 

3.6500 

nsec 

ox  = 

-0.688138 

in. 

OY  =  -8.888165 

in. 

02  * 

8.888508 

»  .  . 

T  = 

3.6700 

nsec 

ox  = 

-0.888148 

in. 

DY  =  -0.008897 

in. 

02  = 

0.881078 

T  * 

3.6900 

nsec 

ox  = 

-0.888142 

in. 

OY  =  -8.888024 

in. 

02  * 

0.881632 

T  = 

3.7100 

nsec 

ox  = 

-8.808142 

in. 

OY  =  8.880058 

in. 

02  = 

0.882167 

’..vv 

T  * 

3.7300 

nsec 

OX  * 

-8.888142 

in. 

DY  =  0.888125 

in. 

02  = 

0.882681 

i  T* 

3.7500 

nsec 

ox  = 

-8.088140 

in. 

DY  =  8.880197 

in. 

02  = 

0.803174 

T- 

3.7700 

nsec 

ox  = 

-0.880137 

in. 

DY  =  0.868266 

in. 

02  = 

8.883643 

1  .  . 

T  =? 

3.7900 

nsec 

ox  = 

-8.088134 

in. 

OY  =  8.868338 

in. 

02  * 

8.884886 

T  = 

3.8100 

nsec 

ox  * 

-0.808138 

in. 

OY  =  8.686387 

in. 

02  = 

8.084581 

-  ' 

T  = 

3.8300 

nsec 

ox  = 

-6.888124 

in. 

OY  =  0.868434 

in. 

02  = 

8.804889 

T  = 

3.8508 

nsec 

ox  = 

-0.888118 

in. 

OY  =  0.088472 

in. 

02  * 

0.085246 

T  = 

3.8700 

nsec 

ox  = 

-6.888112 

in, 

OY  =  0.880499 

in. 

02  = 

8.085571 

T“ 

3.8900 

nsec 

ox  = 

-8.008105 

in. 

OY  =  0.868514 

in. 

02  = 

8.865865 

0 

T  = 

3.9100 

nsec 

ox  = 

-0.008897 

in. 

DY  =  0.680518 

in. 

02  = 

8.606125 

T  = 

3.9300 

nsec 

ox  = 

-8.800889 

in. 

OY  =  0.688589 

in. 

02  = 

8.886351 

•  -  . 

:•  t  * 

3.9500 

nsec 

ox  = 

-0.088881 

in. 

DY  =  0.860489 

in. 

02  = 

0.686543 

T  = 

3.9700 

nsec 

ox  * 

•6.608873 

in. 

DY»  0.860458 

in. 

02  * 

0.886699 

T  = 

3.9900 

nsec 

ox  = 

-6.808065 

in. 

OY  3  6.000418 

in. 

02  * 

8.806819 

1  T: 

4.0100 

nsec 

ox  = 

-0.068858 

in. 

OY  =  8.888369 

in. 

02  = 

6.086965 

1“ 

T* 

4.0300 

nsec 

ox  = 

-0.808851 

in. 

OY  =  8.806313 

in. 

02  = 

8.006953 

T  = 

4.0500 

nsec 

ox  = 

-6.888844 

in. 

OY  3  8.088252 

in. 

02  = 

6.806966 

T  * 

4.0700 

nsec 

ox  = 

-0.000838 

in. 

DY  3  0.000187 

in. 

02  = 

6.686944 

T  = 

4.0900 

nsec 

ox  * 

-8.888832 

in. 

DY=  0.080121 

in. 

02  = 

6.686884 

T  * 

4.1100 

nsec 

ox  = 

-0.680028 

in. 

OY  3  0.006056 

in. 

02  * 

0.006789 

*  .  *  .  -  ‘ . 

— 

|  T  = 

4.1300 

nsec 

ox  * 

-8.886824 

in. 

OY*  -8.888007 

in. 

02  * 

0.086659 

T  = 

4.1500 

nsec 

ox  * 

-0.888021 

in. 

OY*  -0.008065 

in. 

02  * 

0.606494 

T  * 

4.1700 

nsec 

ox  * 

-6.800020 

in. 

DY*  -0.888118 

in. 

02  * 

0.806294 

T  = 

4.1900 

nsec 

ox  = 

-0.060819 

in. 

DY  *  -8.888162 

in. 

02  = 

0.006862 

T  = 

4.2100 

nsec 

ox  = 

-0.008819 

in. 

OY*  -0,000197 

in. 

02  = 

6.885796 

T  = 

4.2300 

nsec 

ox  = 

-6.008820 

in. 

OY  *  -0.088222 

in. 

02  * 

8.085560 

T  = 

4.2500 

nsec 

ox  = 

-6.060822 

in. 

OY*  -0.088236 

in. 

02  * 

0.805174 

T  * 

4.2700 

nsec 

ox  = 

-0.008025 

in. 

DY*  -6.808238 

in. 

02  » 

0.884818 

\  T  = 

4.2900 

nsec 

ox  = 

-8.080029 

in. 

OY*  -8.080228 

in. 

02  = 

6.804436 

T  = 

4.3100 

nsec 

ox  = 

-0.080834 

in. 

OY  =  -0.800207 

in. 

02  = 

8.884027 

T  = 

4.3300 

nsec 

ox  = 

-8.068839 

in. 

OY  *  -0.088175 

in. 

02  = 

0.883593 

* „ 

T  = 
T  * 

4.3500 

nsec 

ox  * 

-8.808845 

in. 

OY  =  -0.800132 

in. 

02  * 

0.803135 

4.3700 

nsec 

ox  = 

•0.008052 

in. 

DY*  -6.606881 

in. 

02  * 

8.882657 

• 

T  = 

4.3900 

nsec 

ox  * 

-0.880858 

in. 

OY  =  -0.088021 

in. 

02  * 

0.002158 

T  * 

4.4100 

nsec 

ox  = 

-6.088065 

in. 

OY  *  6.000844 

in. 

02  = 

0.001641 

.  ■ 

T  = 

4.4300 

nsec 

ox  = 

-8.008872 

in. 

DY  =  0.806114 

in. 

02  * 

8.801109 

T  = 

4.4500 

nsec 

OX  = 

-8.088679 

in. 

OY  *  6.880187 

in. 

02  = 

0.886562 

T  = 

4.4700 

nsec 

ox  = 

-8.068686 

in. 

OY  *  6.606268 

in. 

02  = 

8.668885 

*  .*  '  * 

T  * 

4.4900 

nsec 

ox  = 

-6.608893 

in. 

DY  *  6.866331 

in. 

02  * 

-8.080563 

T  * 

4.5100 

nsec 

ox  = 

-0.008899 

in. 

DY  *  0.608399 

in. 

02  * 

-6.081138 

T  * 

4.5300 

nsec 

ox  * 

-0.088104 

in. 

DY  >  6.860462 

in. 

02  = 

-0.801719 

•  vv 

T  = 

4.5500 

nsec 

ox  = 

-0.088189 

in. 

OY  *  8.880517 

in. 

02  = 

-0.062383 

•-V-V 

T  * 

4.5700 

nsec 

ox  * 

-0.666113 

in, 

OY  *  8.060565 

in. 

02  « 

-0.802888 

,v\-. 

T  = 

4.5900 

nsec 

ox  = 

-8.880116 

in. 

OY  *  0.660602 

in. 

02  * 

-8.003472 

*-V-V 

118 

•**•*.**  *. 
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\  •••V.v.-V-VwV-V 

s’  V" 

'•-V-'; 

•  *«.-'•  \  *,  ■>,  »,.**'*  »  * . 

T 

4.6186 

nsec 

OX  = 

-8.888119 

in. 

OY  * 

8.868629 

in. 

DZ  = 

-6.684853 

T 

= 

4.6388 

nsec 

OX  = 

-8.888119 

in. 

OY  = 

8.888644 

in. 

OZ  = 

-8.884629 

7 

s 

4.6588 

nsec 

ox  = 

-8.888128 

in. 

OY  = 

8.886648 

in. 

OZ  = 

-8.885196 

T 

= 

4.6788 

nsec 

ox  = 

-8.888119 

in. 

OY  = 

8.868648 

in. 

OZ  = 

-8.885754 

T 

3 

4.6988 

nsec 

ox  = 

-8.888117 

in. 

OY  = 

8.688628 

in. 

OZ  = 

-8.886299 

T 

3 

4.7188 

nsec 

ox  = 

-8.888114 

in. 

OY  = 

6.868596 

in. 

OZ  = 

-8.886831 

T 

S 

4.7388 

nsec 

ox  = 

-8.888118 

in. 

OY  = 

8.888556 

in. 

DZ  = 

-8.887345 

T 

= 

4.7588 

nsec 

ox  * 

•8.888185 

in. 

OY  = 

6.868561 

in. 

DZ  = 

-8.087841 

T 

s 

4.7788 

nsec 

ox  = 

-8.888188 

in. 

OY  = 

8.888445 

in. 

DZ  = 

-8.888318 

T 

= 

4.7988 

nsec 

ox  = 

•6.888893 

in. 

OY  = 

8.881383 

in. 

DZ  = 

-8.888772 

T 

3 

4.8188 

nsec 

ox  = 

-8.886886 

in. 

OY  = 

8.888318 

in. 

DZ  = 

-0.869282 

T 

3 

4.8388 

nsec 

ox  = 

-8.888879 

in. 

OY  * 

8.888251 

in. 

O 

r*a 

II 

-8.809688 

T 

= 

4.8588 

nsec 

ox  = 

-8.888871 

in. 

OY  = 

8.888184 

in. 

DZ  = 

-8.809987 

T 

3 

4.8788 

nsec 

ox  = 

-8.888862 

in. 

OY  = 

8.888128 

in. 

DZ  = 

-0.818337 

T 

= 

4.8988 

nsec 

DX  = 

-8.888854 

in. 

OY  = 

8.888868 

in. 

DZ  = 

-8.618658 

T 

= 

4.9188 

nsec 

DX  = 

-6.888845 

in. 

OY  * 

8.886885 

in. 

DZ  = 

-0,818948 

T 

= 

4.9388 

nsec 

ox  = 

-8.888837 

in. 

OY  = 

-8.888842 

in. 

OZ  = 

-8.811285 

T 

= 

4.9588 

nsec 

DX  = 

•8.888828 

in. 

OY  = 

-8.886879 

in. 

DZ  = 

-8.811429 

T 

3 

4.9788 

nsec 

ox  = 

-8.888828 

in. 

OY  = 

-8.888187 

in. 

OZ  = 

-0.811618 

T 

3 

4.9988 

nsec 

DX  = 

-8.888813 

in. 

OY  = 

-8.888124 

in. 

OZ  = 

-0.811771 

T 

3 

5.8188 

nsec 

DX  = 

-8.888886 

in. 

OY  = 

-8.888129 

in. 

DZ  = 

-8.811888 

T 

3 

5.8388 

nsec 

OX  = 

8.888688 

in. 

OY  = 

-8.688123 

in. 

02  = 

-0.811969 

T 

3 

5.8588 

nsec 

DX  = 

8.888886 

in. 

OY  = 

-6.688185 

in. 

OZ  = 

-8.012812 

T 

3 

5.8788 

nsec 

OX  = 

8.888811 

in. 

DY  = 

-8.888876 

in. 

DZ  = 

-0.812017 

T 

3 

5.8988 

nsec 

DX  = 

8.888814 

in. 

OY  = 

-8.888837 

in. 

OZ  = 

-0.011985 

7 

3 

5.1188 

nsec 

OX  = 

8.888817 

in. 

OY  = 

8.868812 

in. 

DZ  = 

-0.011914 

7 

3 

5.1388 

nsec 

DX  = 

8.881819 

in. 

OY  = 

8.881498 

in. 

DZ  = 

-0.011886 

7 

3 

5.1588 

nsec 

OX  * 

8.888819 

in. 

OY  * 

8.881561 

in. 

DZ  = 

-8.811659 

7 

= 

5.1788 

nsec 

OX  = 

8.688819 

in. 

OY  = 

8.881629 

in. 

DZ  = 

-8.011474 

7 

3 

5.1988 

nsec 

OX  = 

8.688818 

in. 

DY  = 

8.861788 

in. 

OZ  = 

-8.811252 

7 

= 

5.2188 

nsec 

DX  = 

8.888816 

in. 

DY  = 

8.881771 

in. 

DZ  = 

-0.818992 

7 

3 

5.2388 

nsec 

OX  = 

8.888813 

in. 

OY  = 

8.881842 

in. 

OZ  = 

-0.018696 

7 

3 

5.2588 

nsec 

OX  = 

8.888889 

in. 

OY  = 

8.881989 

in. 

DZ  = 

-0.010363 

T 

1 

3 

5.2788 

nsec 

OX  = 

8.888884 

in. 

DY  = 

8.861971 

in. 

DZ  = 

-0.009996 

7 

3 

5.2986 

nsec 

OX  = 

-8.888881 

in. 

OY  = 

8.882827 

in. 

OZ  = 

-0.889594 

7 

3 

5.3188 

nsec 

DX  = 

-6.888887 

in. 

DY  = 

8.882874 

in. 

DZ  = 

-0.009160 

7 

= 

5.3388 

nsec 

DX  = 

-6.888813 

in. 

OY  = 

8.682112 

in. 

OZ  * 

-0.688694 

7 

3 

5.3588 

nsec 

OX  = 

-8.888819 

in. 

OY  = 

8.882139 

in. 

DZ  = 

-8.088198 

7 

3 

5.3788 

nsec 

DX  = 

-8.888826 

in. 

OY  = 

8.882155 

in. 

OZ  = 

-8.887674 

7 

3 

5.3988 

nsec 

DX  = 

-8.888832 

in. 

DY  = 

8.882166 

in. 

DZ  = 

-0.887122 

7 

3 

5.4188 

nsec 

DX  = 

-8.888838 

in. 

OY  = 

8.882152 

in. 

OZ  = 

-8.806544 

7 

= 

5.4388 

nsec 

DX  = 

•8.888644 

in. 

DY  = 

8.882133 

in. 

DZ  * 

-0.805943 

7 

= 

5.4588 

nsec 

OX  = 

-8.888858 

in. 

OY  = 

8.882184 

in. 

DZ  = 

-8.085318 

7 

3 

5.4788 

nsec 

OX  = 

-8.668855 

in. 

OY  * 

8.862864 

in. 

DZ  = 

-0.884674 

7 

= 

5.4988 

nsec 

OX  * 

-8.888859 

in. 

OY  = 

8.882815 

in. 

DZ  = 

-0.884818 

7 

3 

5.5188 

nsec 

OX  = 

-8.868863 

in. 

OY  * 

8.881959 

in. 

OZ  = 

-8.883329 

7 

3 

5.5388 

nsec 

OX  * 

-8.888866 

in. 

DY  = 

8.861898 

in. 

DZ  = 

-0.882634 

7 

3 

5.5588 

nsec 

OX  = 

-6.888868 

in, 

OY  = 

8.881832 

in. 

OZ  = 

-8.881927 

7 

3 

5.5788 

nsec 

DX  = 

-8.888869 

in. 

OY  = 

8.881764 

in. 

OZ  = 

-0.881209 

7 

3 

5.5968 

nsec 

DX  « 

-8.868869 

in. 

OY  = 

8.881696 

in. 

OZ  = 

-6.680483 

7 

3 

5.6188 

nsec 

OX  = 

-8.888868 

in. 

OY  = 

8.881638 

in. 

DZ  = 

0.888248 

Ill 


T  = 

5.4308 

nsec 

OX  = 

-0.000044 

in. 

OY  = 

0.801547 

in. 

02  = 

0.800983 

T  = 

5.4500 

nsec 

OX  = 

-0.800043 

in. 

OY  a 

0.881510 

in. 

02  = 

8.881728 

T  = 

5.4700 

nsec 

OX  = 

-0.088059 

in. 

OY  = 

0.881441 

in. 

02  = 

8.002454 

T  = 

5.4900 

nsec 

DX  = 

-0.080054 

in. 

OY  a 

0.001428 

in. 

02  = 

8.003189 

T  = 

5.7108 

nsec 

OX  = 

-8.800849 

in, 

OY  = 

8.001389 

in. 

02  = 

0.803917 

T  = 

5.7308 

nsec 

OX  = 

-0.088842 

in. 

OY  a 

8.081349 

in. 

02  = 

8.884438 

T  = 

5.7500 

nsec 

OX  = 

-0.000035 

in. 

OY  a 

8.001348 

in. 

02  = 

0.005350 

T  = 

5.7700 

nsec 

OX  = 

-0.000027 

in. 

OY  = 

0.881343 

in. 

02  a 

8.804851 

T  a 

5.7900 

nsec 

OX  = 

-0.808019 

in. 

OY  = 

8.081377 

in. 

02  = 

8.804737 

T  a 

5.8100 

nsec 

OX  » 

-0.008010 

in. 

OY  = 

8.001413 

in. 

02  = 

8.087408 

T  = 

5.8300 

nsec 

OX  = 

-8.000881 

in. 

OY  = 

8.881439 

in. 

02  = 

8.008041 

T  = 

5.8500 

nsec 

OX  = 

0.008088 

in. 

OY  a 

8.881485 

in. 

02  = 

8.808494 

T  a 

5.8700 

nsec 

OX  = 

8.008017 

in. 

OY  = 

8.881539 

in. 

02  = 

0.089305 

T  = 

5.8900 

nsec 

OX  = 

0.001125 

in. 

OY  a 

8.801599 

in. 

02  = 

0.809892 

T  a 

5.9100 

nsec 

OX  = 

0.808834 

in. 

OY  = 

8.001445 

in. 

02  = 

0.818454 

T  = 

5.9308 

nsec 

OX  a 

8.000842 

in. 

OY  = 

0.081734 

in. 

02  a 

9.018998 

T  = 

5.9500 

nsec 

DX  = 

8.000049 

in. 

OY  = 

0.801884 

in. 

02  = 

8.011497 

T  = 

5.9700 

nsec 

OX  = 

8.088854 

in. 

OY  = 

0.881873 

in. 

O 

m 

II 

0.811974 

T  = 

5.9900 

nsec 

OX  = 

0.088842 

in. 

OY  = 

8.801939 

in. 

02  = 

0.012420 

TrnX(X)  =  3. 7100  nsec  DNH(X)  «  -0,000142  in. 
7MAX<Y)  a  5.3900  nsec  DMAX(Y)  =  0.002140  in. 
TMAX<2)  =  1.8300  nsec  Dmx<2)  =  0.014337  in. 
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